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ЮБИЛЕЙ УЧЕНОГО 
АММТУЕВЕАВУ ОЕ ТНЕ 5СТЕМТТ5Т 


85 лет Чернышеву Юрию Олеговичу — инженеру, педагогу и ученому 


11 июля 2022 г. исполняется 85 лет со дня рождения заслуженному деятелю науки 
РФ, почетному профессору ДГТУ, профессору кафедры «Автоматизация 
производственных процессов» Донского — государственного — технического 
университета (ДГТУ), главному научному сотруднику научно-исследовательского 
центра Краснодарского высшего военного училища имени генерала армии 
С. М. Штеменко Юрию Олеговичу Чернышеву. 63 года своей жизни он посвятил 
научной, педагогической и научно-организационной деятельности. 

Юрий Олегович Чернышев окончил в 1959 году Таганрогский радиотехнический 
институт, в 1967 году защитил кандидатскую, а в 1984 году — докторскую 
диссертации. 

С 1959 по 1980 годы работал в г. Ростове-на-Дону и занимал должность доцента 
' Ростовского института сельскохозяйственного машиностроения (РИСХМ), а также 
старшего инженера — начальника группы Ростовского высшего военного командно-инженерного училища 


ракетных войск имени Главного маршала артиллерии М. И. Неделина (РВВКИУ РВ). С 1981 по 1983 гг. работал 
директором Северо-Кавказского филиала Российского отделения Государственного проектно-технологического 
института по механизации учета и вычислительных работ ЦСУ СССР. 

С 1983г. Юрий Олегович занял должность заведующего кафедрой «Электротехника и автоматизация 
производственных процессов» Ростовской государственной академии сельхозмашиностроения (РГАСХМ). В 1988 г. он 
основал в этом вузе кафедру «Прикладная математика и вычислительная техника» и возглавил ее. После присоединения в 
2009 году РГАСХМа к ДГТУ и по настоящее время работает профессором кафедры «Автоматизация производственных 
процессов» Донского государственного технического университета и главным научным сотрудником научно- 
исследовательского центра Краснодарского высшего военного училища имени генерала армии С. М. Штеменко. 
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Аннотация 

Введение. Статья посвящена проблеме исследования связанных нестационарных термоэлектроупругих полей в 
пьезокерамических конструкциях. Кратко излагаются основные подходы, связанные с построением общего 
решения при исследовании исходных несамосопряженных уравнений, описывающих рассматриваемый 
процесс. Целью работы является построение нового замкнутого решения осесимметричной задачи 
термоэлектроупругости для длинного пьезокерамического цилиндра. 

Материалы и методы. Рассматривается длинный полый цилиндр, электродированные поверхности которого 
подключены к измерительному прибору с большим входным сопротивлением. На цилиндрических 
поверхностях пластины задана температура, изменяющаяся по времени. Используется гиперболическая теория 
термоэлектроупругости Лорда-—Шульмана. Замкнутое решение построено с помощью обобщенного метода 
конечных интегральных преобразований. 

Результаты исследования. Разработанный алгоритм расчета позволяет определить напряженно- 
деформированное состояние цилиндра, его температурное и электрическое поля. Кроме того, появляется 
возможность исследовать связанность полей в пьезокерамическом цилиндре, а также проанализировать 
влияние релаксации теплового потока на рассматриваемые поля. 

Обсуждение и заключения. Использование допущений о равенстве компонент тензора температурных 
напряжений и отсутствие влияния температуры на электрическое поле позволяет сформулировать 
самосопряженную исходную систему уравнений и построить замкнутое решение. 


Ключевые слова: термоэлектроупругость, гиперболическая теория, нестационарная связанная задача, длин- 
НЫЙ пьезокерамический цилиндр, конечные интегральные преобразования. 
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Введение. В последнее время широкое распространение получили технические устройства различного 
назначения, изготовленные из пьезокерамического материала. Здесь особое место занимают приборы, работа 
которых основана на эффекте связанности упругих, электрических и температурных полей [1]. Для описания их 
работы с учетом  связанности полей в настоящее время разработаны различные теории 
термоэлектроупругости [2—4]. При этом для более качественного описания и оценки нестационарных процессов 
в конструкциях возникает необходимость построения аналитических решений. Однако математическая 
формулировка рассматриваемых задач включает систему несамосопряженных дифференциальных уравнений в 
частных производных, интегрирование которых связано с большими математическими трудностями. 

Для решения данной проблемы, как правило, исследуют уравнения в несвязанном виде [5, 6], 
анализируются бесконечно длинные тела [7-11] или рассматриваются задачи термоэлектроупругости в 
квазистатической постановке [12, 13]. 

В настоящей работе рассматривается связанная динамическая задача термоэлектроупругости для 
бесконечно длинного полого пьезокерамического цилиндра. В результате преобразования исходных расчетных 
соотношенийудается сформировать самосопряженную систему уравнений, интегрирование которой 
осуществляется методом неполного разделения переменных в виде обобщенного конечного интегрального 
преобразования [13]. 

Материалы и методы. Пусть полый, длинный, незакрепленный в радиальной плоскости, 


пьезокерамический цилиндр занимает в цилиндрической системе координат (+.,0,=) область ©: 
{а<к <Ь,0<9<2л, © << <5}. На цилиндрических поверхностях задана температура в виде следующих 


нестационарных функций (граничные условия 1-рода) — в(1.)(к=а), @>(&.) (1. =) . Внутренняя 


электродированная поверхность заземлена, а внешняя подключена к измерительному прибору с большим 
входным сопротивлением (режим электрического холостого хода). 

Математическая формулировка рассматриваемой осесимметричной задачи в безразмерной форме включает 
дифференциальные уравнения движения, электростатики и теплового баланса на основании гиперболической 
теории Лорда-Шульмана, а также краевые условия [2, 7, 15]: 


92 
ей В Ой ей уни 8 Рыб. (1) 
д’ ОГ г ог г 2 
у28 и ке СИ га, У© =0, 
дг ог г ог 
09 [2 2? 
—-| —+ О+а.УИ)=0; 
дк [2 С 20 
ви а а 
г=К,1 2; + а -- 7 9=0,9,,=®,0,„, =5,, (2) 
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1 0 т ры 
00-00 е (3) 
01 [41 


где (бр = [тра ь =, {вое} = |0", Т,,@% пер ый ы м _ в. 


33 33 


р® 


2 2 
1 71 езз _ @@зз _ &з6@зз Ь |Сз —_ Бу _ Св 
а ‚ Ч ‚ 4 =К ‚ а ‚ а 
С;зЕзз 23з37зз А\ р Л С. р С» 

И (к.,ь), 9 (..ь), 9..5) — соответственно радиальная составляющая вектора перемещений, потенциал 
электрического поля и приращение температуры в размерной форме (©^"(#.,1.) =Т(&;,1,)-Т,(.) — текущая 
температура и температура первоначального состояния тела); С.,р,е„,=. — модули упругости, плотность, 
пьезомодули и коэффициент диэлектрической проницаемости электроупругого анизотропного материала 


(", $5ЕЪ 3) ;У7з— компоненты тензора температурных напряжений (у, =Са,, у. =Сцза ); 


Л, К, а — коэффициенты теплопроводности, объемной теплоемкости и линейного температурного расширения 


материала; 8,—щ компонента тензора пирокоэффициентов; В.,—щ время релаксации; ®,— известная в 


е1 

в 9 1 
начальныи момент скорость изменения температуры; У = д НР 
гг 


В случае заземления внутренней поверхности пьезокерамического элемента, электрическое напряжение 


У(.) определяется потенциалом на его внешней поверхности: 


У(1.) =9(1,1.). (4) 
При построении общего решения на первом этапе исследования в результате интегрирования уравнения 
электростатики определяется радиальная компонента вектора напряженности электрического поля: 


ми 0 р 


т г. (5) 


где р, — постоянная интегрирования. 
Подстановка (5) в (1)-(3) позволяет сформулировать новую задачу относительно функций И (",1),© (1,1) ; 


При этом условие отсутствия радиальной составляющей вектора индукции электрического поля на внешней 
цилиндрической поверхности элемента (последнее равенство (2) выполняется в случае 2 =0, а условие 


заземления внутренней поверхности ( ф,_, =0) удовлетворяется в результате интегрирования (5). 


На следующем этапе решения неоднородные краевые условия (2) приводятся к виду, позволяющему в 
дальнейшем использовать процедуру неполного разделения переменных методом конечных интегральных 


преобразований. Для этого вводятся новые функции и(г,!), №(г,!) связанные сИ (г,!),©(г,!): 
О (т.г) = Н, (г.р) -+и(т,г),©(г,г)= Н, (г,1)+М (т, г), (6) 
где Н, (и,г) = Л (г) А(1,1) + (К) А(К,1)+ № (г) <, (1) + Л, (г) (1), 
Н, (гг) = № (г) (1) + , (г) < (1), 
12)... 1,(г) — дважды дифференцируемые функции, А(т.г)_„, = (1+а, -а; а, И (т.г). 
Подстановка (6) в расчетные соотношения (1)-(3) относительно функций И (",!),©(т,г) при выполнении условий: 
(1+а, )УН, +(а, -ЮН, =А(,#) Нк =@, Ну =4, (1) 


г=ЮЛ ? 
позволяет получить новую краевую задачу относительно функции и (, г) р М№(", 1) : 


ди и ОМ № 1 2би м 


У р +В +В Е 8 
бе ‘о бд Ак (1+а,) 0’ ' - 
ОМ 1 5? 
У -[бьр-2 аумчаун) = 5 
г= Ю,1 Уи=0 , М=0: (9) 


=би=-Н, (.0),М=-Н, (,0), (10) 
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би _ СН, (г,1) М 6 @Н, (г,г)_ 
61 в м ^ $ 
О а. 
где Е, =-У = [Сев ен, +а,УН,), 
2 
Е ей е +Ь, р +Ь, -ь - — 
дг г дг г (1+а,) м 
ь _ а, +а.а, ь _ @&-1 в -ь+ 26 
1+а, Г 1+а, . 1+а, 


Здесь следует отметить, что А(х Г) является функцией перемещений цилиндрических поверхностей цилиндра. 


Первоначально А(т, 1) приравнивается к нулю с последующим ее определением и уточнением Н,,Ё,Ё. 
Дальнейшие преобразования расчетных соотношений (8)-(10) связаны с использованием следующих 
допущений: В, =1,Б, =0 и введения термоупругого потенциала 
№=УВ. (11) 
Условие В =1 можно принять без большой погрешности, поскольку для пьезокерамических материалов 
Ь =0,94--0,98, а зависимость В. =О выполняется в случае равенства компонент тензора температурных 


напряжений (у, =У., ) и отсутствие влияния температуры на электрическое поле ( з. =0). 


В результате формируется следующая задача относительно и 1 В В (. Г) : 


д д 1 2и 
— Уи-+Ь УВ =Ё, 12 
ОГ > бк (1+а, ) ОР - ие 
д 12) 27 
був в. вчнии)= К; 
г=А,1 Уи =УВ =0 ; (13) 
=Ои=-Н, (,0) ‚УВ=-Н, (*,0}, (14) 
би _ СН, (г,1) 72) ув-@,—2? (^,1) | 
Ге Ге Ге Ге 
ОН 12) 5? 
где Е, =- = +2 ое, +а,Н,),УН, = № (г)@ (1) № (г) ®, (1). 


Начально-краевую задачу (12)—(14) решаем, используя структурный алгоритм обобщенного конечного 
интегрального преобразования (КИП) [14]. При этом для данной задачи удается использовать 


однокомпонентное неизвестное ядро преобразований К (^, г) } 


{С , С, (1) = ии) вк} рьтак, (15) 


{и(т,г),В(к.г)} = >16 (д.0, са, вук (лк, (16) 


: 
К, | — [ка, к)? ; 
К 


где Я, — собственные значения, образующие счетное множество. 


В результате использования алгоритма КИП [14] получаем задачи относительно ядра преобразований К (1, з г) : 


"КА К (Я, 
4 а ( -?-( ых Ка,=0. (17) 
аг г @ г 
Г ЕЛ УК(А,—)=0, (18) 
и трансформант С (4,,1),С.(4,,Ю: 
Я: 1 4С, 
С, -Е - “ =Ён › (19) 


(1+а,) “ (1+а,) а 
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2 А 4? | 
Я, С, - а (а.б,+а,С,)=Ен , 


АС", ас, 
=” бы =би, —“=Си; (20) 


=0С,=С 
ИИ 


10 › Я 


где{Ен, Ен} = [{Е ‚Е } К(Аьгтаг, {Си 6} =-|{Н, (т,0),Н, (›,0)} КС „гта, 


4 в Г ©Н, (г.г) | ЯН, (+) 
{Си» бью} «о аа 5 К(А,,гугаг. 


Общее решение задачи (17), (18) имеет вид: 
К(А,") =\ (4). Л, (Яг)- Ло (4) (4). 20 
Здесь собственные значения 7, определяются с помощью следующего трансцендентного уравнения: 
% (а) (Ав) - (4) (Ав) =0. 
Система дифференциальных уравнений (19) приводится к следующему разрешающему уравнению 4-го 


порядка относительно С,(,,№): 


И а, м 
[2 +, аг +Ь, аг +В Р +, |6; =Рн, (22) 


ог 


1 Я: 
ира : | ! ыы 
В а’ ар В а, а. В 


Поскольку характеристическое уравнение, соответствующее (22), 
КНЬ К +Ь, К? +В К-Ь,, =0, 


(1+ ь 
где Е, Ав, +(1+а,) к,|- 2 ый С к», 
1 


будет действительное, то оно из условия осциллирующего решения для С, (1,1) имеет два действительных 


корня (А, К., ) и два комплексно сопряженных корня 
(К; ЕС, К, =9-№). 
В этом случае общее решение уравнения (22) имеет вид: 


С, (4,1) = Бы ехр(ки)+ Бу, ехр(к.1)+Р., ехр(К;/) + Бу, ехр(К. 1) (23) 


+, (г) [ехр [^, (1 -т)] —ехр[ А, (:— г)] Ат [в (г)ехр[ ©(1 -- г) | (6 зт(а/—@т)-Ь, с0$(@# + @т)} Ат, 


где, = [(, К ) (®, -в)} +07] в = (+5, а ы к. Бо; = КиК (К, +) 0+0? —@°, 


—К 
2 


К, 41 


В: ==): @= 


Функция С,(Й,,Г) определяется из первого уравнения системы (19). Подстановка полученных выражений 
для трансформант в граничные условия (20) позволяет определить постоянные интегрирования Д,,...О.,. 
Подстановка С1(Я,,1) ‚С. (Я,,В)в (16), (11), (6) позволяет получить окончательные выражения для функций 


О (+, 1) ‚@(х»,2): 
(и) = н, (+ а,Ок (АКГ. 64) 


@(,.:) = Н, (.1)+У ба ,бУк (ак 


11 
На заключительном этапе исследования функции Н, (+,г),Н, (1,1) определяются при решении следующих 


дифференциальных уравнений: 
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ие 5; (25) 
дг Г ог 


что позволяет существенно упростить правые части (Е, Е) расчетных соотношений (8). 


Подстановка выражений для Н,,Ньв (25) позволяет сформировать системы уравнений относительно 
функций 1 (1)... (г) ‚ которые определяются при удовлетворении условий (7). 

Потенциал электрического поля пьезокерамического цилиндра определяется в результате интегрирования 
равенства (5) и удовлетворение предпоследнего граничного условия (2): 


ОН, (г,Е Н! (г 
#- | ц м г. ) на.н, бы (26) 
+в „ов“ +4, У. С, (А.В, (4) [к +10% 4) 
1 1=1 
Н Н о — 
пая ии) вождю 
Г Г = 
Р _ 2К(А, К(4, 
+, »`6, (4. ВВ (ЯК, | й ‚ В, (4-1 и и и 2. 
1=1 |\=А 


В, (А) = [УК (4. 


Полученные расчетные соотношений (24), (26) удовлетворяют дифференциальные уравнения (1) и краевые 
условия (2), (3), т.е. являются замкнутым решением рассматриваемой задачи. 

Результаты исследования. В качестве примера рассматривается радиально поляризованный 
пьезокерамический цилиндр (Р=0,02м, К =0,8) состава РИТ-4, имеющего следующие физические 


характеристики [10]: р = 7500 кг/м °, 


{С»С„Сз} = 113,9; 11,5; 7,43} ж10® нм? {е,»е,} = [-5,2; 15,1} Клум”, 
{уз} = 4,6; 3,9} <10°НИм? °С), =, = 5,6210 фм, а, =2х10 *Кл/м? °С), 
К =3х10° Дж/(м * °С), Л =1,6 Вт/(м °С), В. =5х10с. 


На внутренней поверхности ( 7, =а ) пьезокерамического цилиндра действует температурная нагрузка: 


а) «57 ес. кн.) ‚а (1,)=0, 


тах 


где Н (Г) — единичная функция Хэвисайда ( Н (7) =1 при Г>0, Н(Р)=0 при #<0), Ти =Ти, 1, Тик 
— максимальное значение внешнего температурного воздействия и соответствующее ему время в размерной 
форме (Т", = 373 К (100 °С), 7, =293 К (20°С), #'„ =1 0). 

На рис. 1 представлены графики изменения функций ©” (г,г), ( (т,г), $(г,г) по радиальной координате г в 


различные моменты времени 1. Цифрами 1-3 соответственно обозначены результаты, полученные при 


Л * 
151 


тпах тах (ах — Кр? ах ). 


следующих значениях времени: #=1 41 


тах ? 


Анализ результатов расчета позволяет сделать следующие выводы: 

— достаточно большая величина коэффициента линейного температурного расширения @, пьезокерамического 
материала приводит к быстрому прогреву цилиндра; 

— радиальные перемещения на внутренней цилиндрической поверхности ( 7 = К) на первом этапе исследования 


(1 ЕР ) принимают наибольшие значения с последующим уменьшением в течение времени. Обратная картина 


тах 


наблюдается относительно перемещений при г=1; 
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9° (+,:), °С 


80 


60 


40 


Рис. 1. Графики изменения ©" (+, 1 й и (+, УФ ("1 ) по радиальной координате г в различные моменты времени # 


(2-4 щ›3—154 15 ): а— ©" (и,Р) к; 6—0 (и) -+г;в— ф(и,4)+к 


Степень связанности термоэлектроупругих полей наиболее удобно проанализировать с ПОМОЩЬЮ 


Я? аз й 
коэффициента Б,, = т 1+а, +— | равенства (21). Здесь а, определяет связанность электроупругих полей, а 
Я 


а. 
— влияние скорости изменения объема тела на его температурное поле. 
а 

7 


На рис. 2 представлен график изменения перемещений И (1, 1) по времени { с учетом (сплошная линия) и 


без учета (пунктирная линия) наведенного электрического поля. 
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Следует отметить, что предварительная поляризация пьезокерамики приводит к образованию более 
«жесткого» материала (а, =0,353) и соответственно к уменьшению перемещений при деформировании 
цилиндра. 

Связанностью температурного и электроупругого полей в пьезокерамическом цилиндре можно пренебречь 
за счет малого значения 


48 18х10 < 1. 
а 


(г, Г) х10’ 
10 


2х10° 3х10° 


Рис. 2. График изменения И (1, 1 ) по времени 1 


(сплошная линия — с учетом поляризации, пунктирная — без учета электрического поля) 


На рис. 3 показаны графики изменения электрического напряжения У (1) по времени с учетом (сплошная 


линия) и без учета (пунктирная линия, р, =0) релаксации теплового потока. 


У х!0’ 


0 100 200 300 


Рис. 3. График изменения У (1 ) по времени # (1. „ =0,001 с,Ё„ =196) 


Результаты расчета показывают, что для рассматриваемой задачи уточненную гиперболическую теорию 
Лорда-Шульмана необходимо использовать при большой скорости изменения температурной нагрузки 


ав (1.) 
Ги 


(г ‹ 20,001 с, >5,56х10° к ), а при меныших скоростях — классическую теорию 


* |" =0 
термоэлектроупругости (Д., =0). 

Обсуждение и заключения. Построенное новое замкнутое решение связанной динамической задачи при 
удовлетворении граничных условий теплопроводности 1|-го рода позволяет определить все компоненты 
термоэлектроупругих полей в длинном пьезокерамическом цилиндре. Преимущество представленного 
алгоритма расчета заключается в том, что в отличие от несвязанной постановки задачи, отпадает 
необходимость аппроксимации функции температуры при исследовании уравнения движения. При этом 
действительно влиянием скорости изменения объема пьезокерамического тела на его температурное поле 
можно пренебречь. 
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Расчет силовых и кинематических параметров передаточного 
механизма на основе цепи скручивающихся рычагов 


И.Р. Бондаренко 4, А. А. Волошкив ©, В. С. Перевузник ©, Л.А. Ковалев 
Белгородский государственный технологический университет им. В. Г. Шухова, Российская Федерация, 
г. Белгород, ул. Костюкова, 46 

Ч 1уапбоп485 @гатЫег.га 


Аннотация 

Введение. В работе представлены краткий кинематический анализ, а также применение принципа Даламбера к 
нахождению зависимости между силовыми параметрами в передаточном механизме роботизированного 
манипулятора, построенном из цепи скручивающихся рычагов. Применение данного передаточного механизма 
позволит повысить ресурс привода манипулятора, точность его позиционирования, увеличить рабочую 
нагрузку по сравнению с использованием приводов на гибких связях (скручивающихся нитях). Целью работы 
являлось получение зависимостей между перемещениями элементов цепи, а также их силовыми параметрами, 
необходимых для контроля работы системы управления данными устройствами. 

Материалы и методы. В ходе решения задачи рассмотрен элементарный сегмент (элемент) приводной цепи. 
Для нахождения зависимости между нагрузками в элементе использован принцип возможных перемещений. 
При нахождении кинематических связей между перемещениями проведен краткий геометрический анализ 
элемента передаточной цепи. Для проведения сравнительного анализа полученных аналитических 
зависимостей применён метод имитационного моделирования, реализованный на графической модели звена с 
использованием программного пакета МХ. 

Результаты исследования. В ходе исследования получены зависимости для определения величины момента, 
развиваемого на входном звене в зависимости от внешней рабочей нагрузки и угла его поворота, а также для 
определения линейного перемещения выходного звена. Построена имитационная модель приводного элемента, 
которую возможно применить в динамическом исследовании механизма привода с учетом инертности звеньев. 
Обсуждение и заключения. Из полученных аналитических зависимостей установлено значение угла поворота 
входного звена элемента механизма, при котором достигается максимальное значение величины крутящего 
момента для фиксированной рабочей нагрузки на выходной элемент, а также максимальное линейное 
перемещение выходного звена. Расчетные значения находятся в хорошем соответствии с аналогичными 
величинами, полученными по результатам имитационного эксперимента, что говорит о возможности 
применения аналитических зависимостей при формировании системы управления роботов. Кроме того, данные 
зависимости позволяют обеспечить подбор приводных устройств с необходимыми силовыми показателями. 
Ключевые слова: манипуляторы, гибкие связи, приводной элементарный сегмент, жесткие рычаги, принцип 
возможных перемещений, крутящий момент, имитационный эксперимент. 
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Введение. В последнее время начали приобретать популярность устройства, основанные на закручивании 
кабелей или нитей. Данные изобретения универсальны в использовании и востребованы в различных 
направлениях робототехники. Например, при изготовлении манипуляторов, захватов и экзоскелетов' [1-4]. 
Приводные конструкции, основанные на скручивании гибких элементов, имеют ряд преимуществ. К ним 
относятся: размеры устройства, простота сборки и установки, универсальность. Так, например, мехатронный 
привод на скрученных нитях, представляющий собой электромеханическое устройство с числовым 
программным управлением, способен самостоятельно отслеживать состояние рабочего органа на основе 
данных, полученных с датчика положения вала и датчика силы. 

В представленной работе авторами проводится кинематический и силовой анализ передаточного 
механизма привода, построенного на основе цепи скручивающихся рычагов. Данная конструкция 
предназначена для замены приводных устройств, применяющих гибкие передаточные связи, а именно 
скручивающиеся кабели и нити. 

Экспериментальное устройство, показанное на рис. 1, состоит из электродвигателя и нескольких нитей, 
ориентированных вдоль оси вращения двигателя. Нити связывают выходной вал двигателя с моделируемой 


* Пат. 2020144036 Российская Федерация, МПК Е 16 Н 35/00. Мехатронный привод на скрученных нитях / И. Ю. Гапонов, С. И. Неделчев; 
заявитель и патентообладатель АНОВОУ Иннополис; заявл. 2020.12.30; опубл. 2021.06.24, Бюл. № 12. 
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нагрузкой, создаваемой другим электродвигателем. Управление нагрузкой (линейным двигателем) позволяет 
рассматривать данную систему как пружинно-амортизирующее устройство с регулируемыми параметрами. 


Система 
контроля 


Линейный 
двигатель 


Вращательный Датчик 
двигатель положения 


СИЛЫ 
Рис. 1. Экспериментальная установка [4] 


Благодаря приводу на основе скручивания нитей удалось разместить все приводы и двигатели в протезе 
руки при условии, что обеспечены все основные степени подвижности руки [5-8]. При этом наиболее слабым 
местом системы является гибкий кабель или нить. Износ и растяжение скручивающейся нити являются 
ключевыми недостатками данного типа приводов. 

Альтернативным решением данной проблемы может стать замена скручивающихся нитей или кабелей на 
рычажную систему в виде цепи. Заменим фрагмент скручивающихся нитей или кабелей элементом (сегментом) 
из жестких рычагов (стержней), которые шарнирно соединяют между собой две окружности радиусом К и 
центрами Оги О) (рис. 2 а). Входное звено / удерживается подпятником б и имеет возможность вращательного 
движения относительно центра О’1, находящегося на оси Хо. При этом нижняя окружность 4 (выходное звено) 
может совершать поступательное движение вдоль оси Хо ввиду наличия направляющей 5. В результате 
вращения верхней окружности рычаги 3 и 4 совершают сферическое движение относительно звена /, изменяя 


тем самым расстояние / и перемещая рабочую нагрузку Р (рис. 2 6, в). 


420 


ИЕ 


6) в) 


Рис. 2. Схема сегмента приводной цепи: а — в начальном положении; 
б—в промежуточном положении; в — в предельном положении 


Для практической реализации данного конструктивного решения необходимо наличие системы управления, 


основанной на математической модели приводного устройства, которая отражает кинематические и силовые 
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связи между элементами приводной цепи. 

Целью данной работы было нахождение кинематических и силовых связей в элементарном сегменте 
приводной цепи. 

Материалы и методы. Решение поставленной задачи предлагается найти с помощью принципа возможных 


перемещений [9], позволяющего исключить из рассмотрения внутренние силы в звеньях приводного элемента. 93 


Бир://уези-Чопза.ги 


Бондаренко И.Р. и др. Расчет силовых и кинематических параметров передаточного механизма на основе цепи 


Таким образом, в рассмотрении остаются крутящий момент, приложенный к входному звену, а также 
приложенная к выходному звену рабочая нагрузка Р. Учитывая схему, представленную на рис. 4 6, в качестве 
обобщенной координаты для входного звена принимаем угол поворота ф, а для выходного — линейное 
перемещение < [10]. 

В соответствии с принципом возможных перемещений [11] составим уравнение работ: 

М69+Рдз =0, (1) 
где дф — вращательное возможное перемещение входного звена, дх — линейное возможное перемещение 
выходного звена. 

Расстояние между Н звеньями / и 2 связано с линейным перемещением с звена 2, с учетом длины рычага АВ 


= 1, зависимостью Н =/-— <. Тогда с учетом схемы, предоставленной на рисунке 4 б, можем записать: 
2 


й- Ан) = (1-=)'. 


Величина 7 в зависимости от угла поворота ф определится как: 


г - 
5=1- Ё- оз] © ь (2) 
2 
В соответствии с [11] связь между возможными перемещениями можно установить как: 
4 
9&= 47) ф 90. 
4$ 
Тогда, взяв производную Рот (2) и выполнив преобразования, получим зависимость между возможными 
перемещениями: 
К? 
б:= ыы 66. 
2 
Р-| 2 эт (6) 
2 
Подставив д: в уравнение работ (1), получим: 
К° з1 
М5ф+Р рй 56=0, 


\ _ вы) 


откуда сократив на дф выразим момент М как функцию угла поворота ф в виде: 


РЕ? зтф 
2 
р- вы") 
2 


Для проведения оценки полученных аналитических зависимостей момента (3), развиваемого на 


М = 


(3) 


входном звене и линейного перемещения выходного звена от угла поворота (2), был проведён 
имитационный эксперимент [12, 13], реализованный в программной среде МХ Мазгап [14]. С этой целью 
была построена трехмерная твердотельная модель участка приводной цепи (рис. 3 а), состоящего из двух 
последовательно соединенных элементарных сегментов исследуемого механизма. Один из сегментов 
приводного механизма (рис. 3 6) был использован для проведения расчета со следующими параметрами К 
= 12 мм, [ = 54 мм. При построении модели сферические шарниры были представлены в виде торовых 
шарниров, что позволило выполнить элементы цепи в виде изогнутых колец, соединенных жесткими 
рычагами (рис. 3). Данное конструктивное исполнение позволяет рычагам совершать сферическое 
движение. Так как данная связь обеспечивает рычагам требуемую степень свободы, можно считать еб 


соответствующей сферическому шарниру (рис. 3 в). 
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Рис. 3. Трехмерные твердотельные модели: а — участка цепи из двух сегментов; б — сегмента цепи; 
в — проекции сложно изогнутого кольца 


В разработанной трехмерной модели были назначены кинематические связи, определены центры масс звеньев, 
а также задана рабочая нагрузка Р = 100 Н, после чего была проведена симуляция движения элементарного 
сегмента приводной цепи в МХ Мазап [14, 15]. Параметры модели сегмента приведены в таблице 1. 


Таблица 1 
Параметры моделирования 
Момент инерции Начальная координата 
№ связи в относительно каждой из Масса, центра масс, мм 
симуляции осей, кг ` мм? КГ р. 
х ГУ 17 х У 7, 
Г.001 0 1,63 209 Изогнутое 
[004 5,51 4,18 4,18 0,04 0 1,63 139,5 кольцо сложной 
[007 0 1,63 69,8 формы 
[002 —0,03 —12,9 173,4 
-_ 0,03 12,9 | 173,4 0,02 —— — —_ Рычаг 
Г.005 12,9 1,5 103,7 
Г006 —12,9 1,5 103,7 


Результаты исследования. В соответствии с полученными аналитическими зависимостями для величины 
линейного перемещения 7, а также крутящего момента М, был выполнен расчет, произведенный в среде 
МаСАРО, результаты которого представлены на рис. 4. Расчет был выполнен при следующих исходных 
данных: рабочая нагрузка Р = 100 Н, радиус входного и выходного звена А = 12 мм, длина рычага К = 12 мм, 
диапазон изменения угла поворота ф входного звена в пределах 0-180°. 


<, мм М, Н:мм 
5,6 281,3 
4,7 234,4 
3,75 - 187,5 
2,82 140,7 / 
1,9 93,8 \ 
0,94 46,9 / - 
и 


0 30 60 90 120 150 180 9’ 0 30 60 90 120 150 180 ф,у 


а) 6) 
Рис. 4. Зависимости исследуемых величин от угла поворота: а — линейного перемещения выходного звена; 
б — крутящего момента М на входном звене 
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Из представленных графических зависимостей можно установить, что максимальная величина перемещения 
достигается при угле поворота, приближающимся к 180° и составляет 5,6 мм, а величина крутящего момента 
принимает максимальное значение 281,3 Н-мм при угле поворота ф = 93°. 

В результате имитационного моделирования получен график перемещения выходного звена Г.004 в 
зависимости от поворота приводного звена Г.001 (рис. 5). Перемещение элементарного сегмента рычажного 
механизма происходит до тех пор, пока не произойдет касание правого и левого рычагов в крайнем положении. 
Так как рычаги имеют объем, полный поворот выходного звена невозможен из-за касания боковой 
поверхности, что также отображено на графике (рис. 5). Для достижения полного поворота необходимо 
изменить строение рычага на более сложную форму. На рис. 6 представлена зависимость крутящего момента 
при повороте сегмента цепи расчетной модели на 180°. 


< мм 


144.9 : " | 


143,8 


142,6 


141,6 


140,5 


139.5 
0 60 120 180 Фф/’ 


Рис. 5. График перемещения выходного звена в относительной системе координат 
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Рис. 6. График изменения крутящего момента в процессе проведения имитационного эксперимента 


Значения исследуемых параметров, полученные при помощи симуляции (имитационного эксперимента), 
отличаются от значений полученных теоретически не более, чем на 3 %. 

Обсуждение и заключения. В рамках представленной работы проведены кинематический и силовой анализы 
элементарного сегмента приводной цепи манипуляторного устройства. Полученные аналитические зависимости 
позволили определить максимальные значения линейного перемещения и развиваемого крутящего момента в 
приводном сегменте. Результаты имитационного эксперимента, реализованного в среде МХ Мазбап, и 
последующего сравнительного анализа показали хорошее соответствие полученных аналитических зависимостей 
по отношению к трехмерной динамической модели объекта. Это позволяет сделать вывод о том, что полученные 
аналитические зависимости могут быть использованы при разработке системы управления приводов 
вышерассмотренных манипуляторных устройств. Проведенные исследования показали, что цепь из сегментов 
рассматриваемого рычажного механизма обладает хорошими передаточными свойствами и при достаточной 
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длине может заменить передачу на основе скручивания нити. Однако устройство исследуемой передачи более 
сложное. Несмотря на это, передача на основе скручивания цепи сегментов обладает таким преимуществом, как 
возможность обеспечения большего числа рабочих циклов манипуляторного устройства с сохранением требуемой 
точности позиционирования и повторяемости за счет меньших деформаций элементов цепи. 
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Исследование влияния режимов резания на выходные 
параметры при точении фрезерованием быстрорежущей стали 


Г. В. Матлыгин 4, А. В. Савилов ©, А. С. Патых ©. С.А. Тимофеев ®@ 
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Аннотация 


Введение. Рассматривается вопрос повышения эффективности процесса лезвийной обработки изделий из 
порошковой быстрорежущей стали. Предлагается использовать точение фрезерованием как альтернативу 
классической токарной обработке. Показаны преимущества данного метода. Дан обзор исследований, 
посвященных повышению качества поверхности деталей при точении фрезерованием. Целью работы явилось 
определение влияния режимов резания на шероховатость поверхности при ортогональном методе точения 
фрезерованием порошковой быстрорежущей стали монолитной фрезой. 

Материалы и методы. Применялись методы статистического анализа, основанные на создании 
математической модели прогнозирования микрогеометрических отклонений обрабатываемой поверхности. 
Метод экспериментального исследования использовался для проверки адекватности математической модели. 
Планирование эксперимента производилось по некомпозиционному плану, предложенному Боксом и 
Бенкиным. Эксперимент осуществлялся на токарном обрабатывающем центре с приводным инструментом. В 
качестве материала образцов для проведения эксперимента была применена порошковая быстрорежущая сталь 
ВОНГЕВ $390 М!СВОСТЕАКМ. В качестве режущего инструмента применена монолитная твердосплавная 
фреза. В ходе эксперимента варьировались скорость резания, ширина фрезерования и подача на зуб. Измерение 
шероховатости обработанной поверхности производилось контактным профилометром. 

Результаты исследования. Разработана математическая модель формирования шероховатости поверхности в 
зависимости от параметров обработки. В ходе экспериментов определено влияние скорости резания, подачи 
инструмента и радиальной глубины резания на шероховатость обработанной поверхности. Установлено, что 
зависимость шероховатости от подачи имеет линейный характер на всём исследуемом диапазоне режимов 
резания. В свою очередь зависимость шероховатости от скорости резания и ширины резания имеет 
параболический характер. Полученные результаты позволяют достигать шероховатости обработанной 
поверхности Ка = 1,85 без снижения производительности обработки. 

Обсуждение и заключения. Разработанная математическая модель отражает влияние режимов резания на 
шероховатость поверхности при точении фрезерованием быстрорежущей стали монолитной фрезой. 
Результаты проведённых исследований могут быть использованы для определения оптимальных режимов 
резания, обеспечивающих заданное качество поверхности при изготовлении реальных деталей в 
производственных условиях. Продолжение исследований рекомендуется проводить с контролем 
дополнительных выходных параметров, таких как температура и вибрация. Уменьшение влияния 
регенеративных автоколебаний на шероховатость обработанной поверхности можно достичь за счёт назначения 
режимов резания на базе результатов модального анализа технологической системы. 


Ключевые слова: точение фрезерованием, твердосплавная фреза, шероховатость, быстрорежущая сталь, 
моделирование процесса резания. 
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Оё5си5;опт апа Сопсш51оп5. 'ТБе деуеоре4 таетаНса| то4е! геЙес{5 Фе ппрасЁ оЁ сит» то4ез оп Фе 
зигРасе гоизбпез$ \Веп игилие 6126-5рее4 ${ее] \1 а тшопой1с сиНег. ТБе геза$ оЁР фе сопдаисе4 гезеагсв 
сап Бе изе4 10 деегише Фе орНта! сие то4ез Фа ргоу14е а х1уеп зи! Расе дааШу ш Фе тапа{астаге оф геа1 
рамз ипаег Ше ргодисНоп сопд1опз. [ 15 гесоттеп4е4 {0 сопНпие Ше гезеагсВ \И Ше сопёго] оЁ аа опа! 
ошёриЕ рагатеегз, засН аз {етрегааге ап4 уганоп. Кедисш» Фе еНесе о гезепегануе зе-озсШаНопз оп Фе 
гоиорпе$$ оЁ Ше незме4 злиРасе сап Бе геасве4 ШгоиэН аз$1еие (Бе сие то4ез Базе оп Фе гези $ оРа 
пода! апа1у$15 ОР Ше ргосе5$ зузет. 


Кеумога5$: шги-тИШпо, саг14е сацег, гоиоВпезз, Ш2В-зрее4 5ее, сит зпоаНоп. 


Кипаше иогтайоп: Фе аиФогз \ош@ ПКе 10 ФапК Фе еАКогз апа гемлемег$ Рог Феи: саге айепНоп 0 Ше агйс]е 
ап Гог Ше сопитег!$, \уШсВ паде # роз е ю паргоуе 1$ апаШу. 

Ког сНайоп: С. У. Майуеть, А. У. ЗауПоу, А. 5. Руа(уКВ, 5. А. ТипоЁееу. Заду оЁ Ше еЁес+ оЁ сит то4ез 
оп оц риё рагатеег$ ипдег 611 -5рее4 $ее1 аги-шИШп2. А4уапсеа Епз1теегто ВезеагсВ, 2022, уо1. 22, по. 2, 
рр. 99-106. В&рз://401.0г2/10.23947/2687-1653-2022-22-2-99-106 


Введение. Альтернативой классического процесса резания является точение фрезерованием. Первые 
попытки описания данного вида обработки и его практического применения имели место в первой половине 
ХХв. В 1948 г. Г. И. Грановский описал кинематику обработки цилиндрической вращающейся детали с 
помощью фрезы [1]. В 1964г. А. О. Этин описал различные методы точения фрезерованием, в т. ч. 
фрезерование инструментом с осью, параллельной оси детали, и фрезерование цилиндрической фрезой с 
круговой и тангенциальной подачами [2]. Принято считать основоположниками практического применения 
данного метода Шульца и Шпура. Они разделили операции точения фрезерованием на две группы: 
ортогональные и тангенциальные [3]. Значительный вклад в изучение токарного фрезерования внесли 
Карагузел, Баккал, Будак, которые исследовали механику и термические аспекты токарно-фрезерной 
обработки. Они также разработали модели для проверки геометрии процесса, кинематики и механики, качества 
обрабатываемых деталей в ортогональных и тангенциальных токарно-фрезерных операциях [4]. 
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Актуальной проблемой при внедрении операции точения фрезерованием является достижение высокого 
качества поверхности, в первую очередь — уменьшения шероховатости. Анализ публикаций позволяет 
выделить основные направления в исследовании путей достижения высокого качества поверхности: 

— изучение влияния эксцентриситета [5]; 

— определение профиля режущего инструмента со сменными многогранными пластинами с круглой и 
скругленной режущими кромками [6]; 

— изучение влияния температуры в зоне резания [7]; 

— прогнозирование износостойкости режущей кромки при применении операции  точение 
фрезерованием [8]. 

Следует отметить, что в указанных работах в качестве исследуемых материалов применялись традиционные 
конструкционные стали. Исследования быстрорежущей стали, в том числе полученной методом порошковой 
металлургии, не являются распространенными. Это ограничивает применение преимуществ точения 
фрезерованием для изготовления вращающихся режущих инструментов [9]. Данное направление 
механообработки можно отнести к области, где точение фрезерованием может дать наибольший 
экономический эффект. 

На качество поверхности влияет используемый режущий инструмент. Например, применение концевых 
фрез со сменными многогранными пластинами не всегда позволяет достичь требуемой шероховатости 
поверхности [10]. Важно отметить, что рассматриваемая экспериментальная часть исследований проводится на 
многоцелевых станках. Однако точение фрезерованием также может применяться и на токарных 
обрабатывающих центрах с приводным инструментом [11]. 

Целью данной работы является определение влияния режимов резания на шероховатость поверхности при 
ортогональном методе точения фрезерованием порошковой быстрорежущей стали монолитной фрезой. 

Материалы и методы. Эксперимент проводился на токарном обрабатывающем центре с приводным 
инструментом ОМС МЕРАО0 (рис. 1). 


Рис. 1. Зона обработки при эксперименте 


В качестве исследуемого материала использована порошковая быстрорежущая сталь ВОНГЕВ $390 
МСВОСГЕАМ'. Заготовка представляла собой цилиндр размером 88х300 мм и закреплялась в трёхкулачковом 
гидравлическом патроне с поджатием центром, установленным в задней бабке станка (рис. 2). В качестве 
режущего инструмента использована монолитная твердосплавная фреза У@-1 СМ999 12(В3)х12х32х75, #=5. 


'М390-ВОНГЕВ Виза // уму БоМегаи.ги : [сайт]. ОВГ: Вар:/Ам\у БоШегип.ги/ги/ргодис/т1390/ (дата обращения: 14.12.2021). 
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Обработка проводилась без применения смазочно-охлаждающей жидкости. Измерения шероховатости 
обработанной поверхности производились на стационарном профилометре Тау|ог&Нобзоп Еогт Та1узиг{ 200. 

Химический состав стали ВОНГЕВ $390 МСВОСГЕАМ: 

1,64 % С; 0,60 % 51; 0,30 % Мп; 4,80 % Ст; 2,00 % Мо; 4,80 % У; 10,40 $ У; 8,00 % Со. 

Технологические свойства стали ВОНГЕВ $390 МСВОСТЕАМ, представленные в балльной системе из 5 
баллов, соответствующих максимальным значениям свойств данной группы сталей: 

— прочность на сжатие — 4; 

— шлифуемость — 3; 

— прокаливаемость — 4; 

— вязкость — 4; 
износостойкость — 4. 


7! 


Рис. 2. Схема ортогональной обработки монолитной фрезой УС-1 СМС19910 


Математическая модель отклонения шероховатости поверхности Ка может быть представлена 

уравнением [11]: 
шРа = шр + с! ' ШХ, + с11 №72 + с ШИ, + со ' ШИ? + сз - Шаг + сзз  ша2, 

где /, — подача, мм/зуб; У, — скорость резания, м/мин; а, — ширина резания, мм. 

Переименуем это выражение следующим образом: 

У = 90 ` Хо + 91 ` Хи + 911 ' №2 + 92 Хо + 922 ` ХЗ 4 93 Хз + 9зз' 2, 

где у, = ШшКа; хо — фиктивная переменная; хь, х›, хз — кодированные значения факторов; 9% 9ь 9> 93 — 
постоянные коэффициенты. 

Кодирование переменных осуществляется по следующим соотношениям: 


_ 2+ тих) 


= 1; 
ь Ш в — № тт ь 
2(1И + ШУ, 
= тах 1; 
-. № а — м И тт ь 


о 2(тае + Ша.) а 
° Ш Четах — №1 Чета | 


'Афуапсеа Епэтеетия Кезеагсй 2022. Т. 22, № 2. С. 99-106. 155М№ 2687—1653 


где тах и 


а 
Стах’ “@тах 


— верхние уровни факторов; №/;„, И 


‚Че: 
Стт’ “@тт 


— нижние уровни факторов; /»„, И; 


основные уровни факторов. Значения этих параметров представлены в таблице 1. 


Значения уровней факторов эксперимента 


о’ @ед 


Таблица 1 


Фтор Код Нижний уровень (- | Основной уровень | Верхний уровень 
О (0) (+0 
Ё, мм/зуб Хх 0,04 0,08 0,12 
У, м/мин х 80 90 100 
е, ММ хз 1 2 3 


Эксперимент планируется по некомпозиционному плану, предложенному Боксом и Бенкиным. Этот план 
представляет собой выборки строк из полного факторного эксперимента типа 3^. Он включает в себя 15 опытов. 
Матрица плана приведена в таблице 2. 


Таблица 2 
Матрица результатов планирования эксперимента 
№ Выходной 
пит а параметр 
хо х1 № Хз да Эа | ежу | [№ | № Ва, мкм 
1 1 -1 -1 0 1 0 0 1 110 0,9487 
2 1 1 -1 0 -1 0 0 1 110 1,5583 
3 1 -1 1 0 -1 0 0 1 110 1,2183 
4 1 1 1 0 1 0 0 1 110 1,7555 
5 1 -1 0 -1 0 1 0 11011 1,5674 
6 1 1 0 -1 0 -1 0 11011 1,9084 
7 1 -1 0 1 0 -1 0 11011 1,5198 
8 1 1 0 1 0 1 0 11011 1,9505 
9 1 0 -1 -1 0 0 1 01 1 1,8912 
10 1 0 1 -1 0 0 -1 01 1 1,0914 
и 1 0 -1 1 0 0 -1 01 1 1,6132 
12 1 0 1 1 0 0 1 01 1 1,7664 
13 1 0 0 0 0 0 0 0010 1,5668 
14 1 0 0 0 0 0 0 0010 1,5081 
15 1 0 0 0 0 0 0 0010 1,5046 
Вектор коэффициентов определяется по формуле: 
В = (ХТХ)-1ХТУ. 
Дисперсии 5 2(Ь;) коэффициентов регрессии находятся из выражения: 
52(Ъ,) = си5у, 
где с; — диагональные элементы матрицы ХХ"'; 5?у — дисперсия эмпирической регрессии. 


Среднее значение и дисперсия определяется вектором опытов в центре плана. 
Среднее значение: 
— Хоу 
к 
где У, — вектор значений функции отклика в центре плана. 


Дисперсия: 


) 


52 — общ)? 
У 2 - 
По результатам опытов, проведенных по плану (табл. 2) можно определить постоянные коэффициенты 
уравнения: 
У = 40 Хо + 91 *Х1 4 92 ' Х2 4 9з ' Хз + 912 * Х1 ' Х2 + 913 'Х1 * Хз 4 923 `Х2 ' Хз 
+911 ` 2 + 922 ` 2 + 9зз` ХЗ; 
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по 
1 
Зо = я Ус; 
и=1 
М 
9: =А у ху); 
р 
М 
9и =Б У хихиу;; 
1 
М к М по 
2 2 
9з =вУ жи > р.п Ху 
и=1 1=1 ]=1 ы и=1 
где по — число опытов в центре плана; и — номер параллельного опыта в центре плана; уд, — значение 
функции отклика в 1-ом опыте; М — число опытов в матрице планирования; ] — номер опыта в матрице 
планирования; 1, [ — номера факторов; хь, х; — кодированные значение 1-го и 1-го факторов в ]-м опыте; у; — 


значение функции отклика в ]-м опыте. 
Из этого следует: до = 1,388; 4. = 0,24; 42 = —0,022; 93 = 0,268; 911 = 2,413; 922 = 2,267; 933 = 3,072. 
Таким образом, модель зависимости отклонения шероховатости поверхности от параметров обработки и 
ширины резания будет иметь вид: 
Ва = 1,388 + 2,413 . |2 — 0,24. р. 
Результаты исследования. Зависимости шероховатости поверхности от режимов резания 
представлены на рис. 3. 


а) 6) в) 


Рис. 3. Зависимости отклонения параметра шероховатости К,: 
а — от подачи на зуб {. и от скорости резания У’; 6 — от скорости резания У. и ширины фрезерования а.; 
в — от ширины фрезерования а, и подачи на зуб {}. 


Анализ полученных результатов показывает, что зависимость шероховатости от подачи имеет линейный 
характер на всём исследуемом диапазоне режимов резания (рис. 3 а, в). При этом с увеличением подачи 
шероховатость пропорционально увеличивается, что согласуется с основными положениями теории резания. 
Максимальное значение параметра шероховатости А, составляет 1,85 мкм, что соответствует наиболее 
распространённым требованиям при чистовой механообработке в машиностроении. 

Зависимость шероховатости от скорости резания У, (рис. 3 а, 6) имеет параболический характер с 
максимальными значениями шероховатости в диапазоне скоростей резания 90-95 м/мин. Данный факт можно 
объяснить тем, что в указанном диапазоне скоростей резания имеет место эффект нароста материала на 
режущую кромку. Для подтверждения указанной гипотезы необходимо проведение дополнительных 
исследований с контролем температуры в зоне резания и использованием высокоскоростной видеокамеры для 
наблюдения за процессом стружкообразования и возможного формирования нароста. 

Влияние ширины фрезерования а, на шероховатость поверхности имеет также параболический характер 
(рис. 3 6, в), но при этом максимальные значения шероховатости достигаются при а. = 0,25 мм и а, = 3 мм, а 
минимальное значение — при а, = 1,8-2,2 мм. Можно полагать, что в указанном диапазоне ширины резания на 
процесс обработки в меньшей степени влияют регенеративные колебания, зависящие также от конструктивных 
и геометрических параметров режущего инструмента. 
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В игоге можно констатировать, что при установленных режимах резания удалось достичь шероховатости 
поверхности не выше К, = 1,85 без потери производительности процесса. 

Обсуждение и заключения. Разработана математическая модель, отражающая влияние режимов резания на 
шероховатость поверхности при точении фрезерованием быстрорежущей стали монолитной фрезой. 
Результаты проведённых исследований могут быть использованы для определения оптимальных режимов 
резания, обеспечивающих заданное качество поверхности при изготовлении реальных деталей в 
производственных условиях. 

Рекомендуется проведение дополнительных исследований с контролем температуры в зоне резания и 
вибраций. Назначение режимов резания следует проводить на базе результатов модального анализа 
технологической системы для уменьшения влияния регенеративных автоколебаний на шероховатость 
обработанной поверхности. 
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Аннотация 


Введение. За последнее десятилетие спрос в мире на резервуары высокого давления существенно возрос, 
особенно в таких областях промышленности как авиационная, космическая, химическая и нефтегазовая. 
Находясь под постоянным воздействием высокого внутреннего давления, стенки резервуаров испытывают 
повышенное напряжение, что может стать причиной их внезапного разрушения. Для устранения такой 
возможности и улучшения прочностных характеристик резервуары изготавливают в форме металлических 
цилиндров с внутренним покрытием из композитного материала, состоящего из смолы, армированной 
углеродными волокнами. Цель настоящей статьи заключалась в изучении влияния угла наклона углеродного 
волокна на цилиндрические резервуары и определении величины максимального разрушающего давления с 
использованием метода конечных элементов программы АМЗУ5. 

Материалы и методы. С помощью программы АМЗУЗ создана конечно-элементная модель резервуара, 
состоящего из центральной части в виде металлического цилиндра с внутренним покрытием из полимера, 
армированного углеродными волокнами. По торцам резервуара размещены полусферы со спиральной 
намоткой. Для моделирования слоистого композитного материала использована командная оболочка 
НЕТ. 99, для определения критерия разрушения резервуара высокого давления — теория Тза1-\Уи. 
Результаты исследования. Модель цилиндрического резервуара рассчитывалась для двух видов намотки 
волокон: кольцевой и спиральной при различных углах их наклона. Результаты анализа величины давления для 
различных углов наклона волокон показывают, что, начиная со значения угла 0° и до 45° оно увеличивается, а 
затем до значения 65° — уменьшается. Наиболышее давление, которое может выдержать резервуар, 
армированный углеволокном, составляет 207 МПа при угле наклона волокон + 45°. 

Обсуждение и заключения. Анализ исследований показал, что при угле наклона волокон + 45° максимальное 
напряжение оказалась наименьшим, а максимально возможное разрушающее давление при том же угле 
составило 207 МПа. Из этого следует, что оптимальный угол ориентации волокон для обеспечения безопасной 
работы резервуара высокого давления составляет - 45°, а резервуар из углепластика при том же угле намотки 
волокон имеет максимальную прочность. 


Ключевые слова: резервуар высокого давления, компьютерная модель, угол намотки, композиционное 
покрытие, углеродное волокно, полимерное связующее. 
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Введение. Использование резервуаров высокого давления из углеродных волокон нашло широкое 
распространение в различных областях промышленности из-за своих отличительных свойств, таких как малый 
вес и высокая прочность. Поэтому в последнее время спрос на такие резервуары существенно возрос в тех 
случаях, когда вес играет большую роль [1, 2]. 

Важнейшими областями применения резервуаров такого типа являются аэрокосмическая, авиационная и 
химическая техника. Кроме того, резервуары, армированные волокном, широко используются для увеличения 
напора при транспортировке нефти и газа. Во многих случаях резервуары такого назначения подвергаются 
воздействию высокого внутреннего давления, следствием чего может стать скачкообразное повышение 
напряжения на стенках сосуда и внезапное их разрушение, что наносит большой ущерб материальным и 
людским ресурсам [3, 4]. 

В [5] проведен ряд исследований, направленных только на изучение напорного резервуара из многослой- 
ного композита, в ходе которых определена ожидаемая прочность на смятие. В [6] изучалось поведение враща- 
ющегося составного напорного резервуара под воздействием внутреннего давления и осевой нагрузки; в [7] 
проведены исследования по изучению влияния тепловых нагрузок, действующих на многослойный композит- 
ный напорный бак. В работе [8] исследовалось поведение многоугольного композитного напорного бака пяти 
различных форм под действием внутреннего давления различного режима, а в [9] — изучалась работоспособ- 
ность составного напорного бака при воздействии поперечных нагрузок. 

Проектирование резервуара из композитного материала является сложной проблемой, поэтому для проведе- 
ния полного и точного анализа необходимо выбрать приоритетные факторы. Авторами изучался резервуар вы- 
сокого давления, армированный несколькими слоями углеродных волокон и подвергаемый нагрузке от внут- 
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реннего давления. Цель работы — изучение влияния угла наклона углеродного волокна на цилиндрические 
резервуары и определение величины максимального разрушающего давления с использованием метода 
конечных элементов программы АМЗУ5. 

Материалы и методы. Теоретическое исследование. Резервуар высокого давления из композитных 
материалов (СЕВКР) состоит из центральной части, представляющей собой металлический цилиндр с 
внутренним покрытием, состоящим из полимера, армированного углеродными волокнами. Торцевые 
поверхности резервуара имеют форму полусфер со спиральной намоткой (рис. 1). 


- жнаыыызнннннннннни 
рака 
Рис. 1. Продольное сечение резервуара 

Резервуары, работающие при высоких давлениях и армированные углеродным волокном, изготавливаются 
методом филаментной намотки (рис. 2). Для получения необходимой стабильности армирования волокна 
отправляются на движущуюся тележку при тщательном выборе их координации, а затем наматываются на 
цилиндрическую поверхность. На стабильность координации волокон влияет несколько факторов: температура, 
форма и обработка поверхности, а также степень сцепления волокон с матрицей. Угол наклона намотки 
контролируется скоростью перемещения тележки и скоростью вращения цилиндрического барабана. С целью 


получения высоких эксплуатационных и прочностных свойств внутренняя цилиндрическая поверхность 
резервуара покрывается несколькими слоями волокон [10]. 


Угол наклона волокна 


Цилиндрический я и 
| || 
барабан ши ооо 
Ш «— Полимер | 
Движущаяся тележка ——— —=—=— 
2111--: — Волокна 


Рис. 2. Схема намотки волокон 


Угол наклона волокон оказывает большое влияние на свойства резервуаров, поэтому важным вопросом 
является поиск соответствующего угла для каждой части резервуара. Угол ориентации волокон определяется 
необходимой величиной трения между волокнами и слоем композитного материала, как показано в 


соотношении: 
—_ РВ п 
а&(К)=зт" К +5 в (1) 
К К -№ 
где К — расстояние между центром и точкой слоя; К — радиус центральной оси; К, — радиус по 


касательной к поверхности цилиндра при д = 0 [8]. 

Критерии разрушения Тба!-У’а. При исследовании и моделировании резервуара высокого давления, 
выполненного из композитного материала, использовались критерии разрушения по теории Тза1-Уа [11]. 
Реализация уравнения (1) необходима для исследования ожидаемого разрушения ортотропных материалов по 
теории Тза1-У/и: 


2 5 — 
Вой +0» + 5о› +В 10 + Еуть +2500, =1. 
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Упругие свойства определяются четырьмя независимыми константами: Ел, Е?2, С1о, Ут», 
представленными в таблице 1. 


Таблица 1 
Характеристики композитного материала Г.У5052/Т300 
Свойства Углеродное волокно Т300 Эпоксид [./У5052 
Модуль упругости 230 ГПа 3,0 ГПа 
Прочность на растяжение 3,5 ГПа 71,0 ГПа 
Плотность 1760 г/см? 1,14 Г/см? 
Силы рассчитываются по следующим уравнениям: 
1 1 11 1 
р. о: 
1 1 1 
Е уу: 75" 12 по ЕЕ», 
где Х, — растягивающая сила в продольном направлении; У — растягивающая сила в поперечном 
направлении; Х, —Щ сила давления в продольном направлении; У — сила давления в поперечном 


направлении; 5 — сила сдвига. 
Максимальное разрушающее напряжение достигается при выполнении одного из следующих 

соотношений [12]: 

[е 


[@} 
Е 


Х 

Свойства материалов и моделирование методом конечных элементов. Композитный материал, 
использованный для исследуемого резервуара, представляет собой полимер, армированный углеродным 
волокном типа Т300, а в качестве полимерного материала применяется эпоксид Г/У5052. Композитные 
материалы имеют ортотропную природу, поэтому процесс их моделирования конечными элементами сложнее, 
чем изотропных материалов, таких как алюминий и сталь. 

На рис.3 представлена конечно-элементная модель резервуара высокого давления из углепластика, 
внутренний слой которого состоит из алюминиевого сплава, армированного восемью слоями композитного 
материала Сагфоп/Ероху Т300/.У5052. 

Модель имеет следующие размеры: 

— длина резервуара —1200 мм; 


диаметр резервуара по центру — 300 мм; 
— общая толщина — 64 мм; 
— толщина одного слоя — 6,5 мм; 


толщина футеровки — 0,12 мм. 
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Рис. 3. Конечно-элементная модель резервуара СЕКР 
Для исследования резервуаров высокого давления из композитных материалов очень важно выбрать 
110 подходящий тип конечных элементов. Программа АМЗУ5 содержит конечные элементы ЗНЕГ, и ОГ, 


необходимые для моделирования слоистых композитных материалов. Для исследования использовали 
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программу ЗНЕМ, 99, которая ускоряет расчёт конструкции, имеющей до 250 слоев. Этот элемент представляет 
собой многоуровневую линейную структуру с восемью узлами и шестью степенями свободы и позволяет 
пользователю определять гибкость, наклон слоев и плотность каждого слоя. 

При изучении композитных материалов формирование слоев является одним из основных вопросов, 
поскольку каждый слой имеет свой угол наклона, а волокна в каждом слое имеют различные углы наклона, 
поэтому свойства каждого слоя необходимо определять отдельно. 

Для формирования слоев требуются характеристики материала, количество слоев, угол наклона волокон, 
толщина слоя и количество точек интегрирования в каждом слое. Приведем сравнительные сведения о 
свойствах применяемых материалов: СЕБР и алюминия 6061. Плотность: СЕВР — 1570 кг/м°, алюминия — 
2750 кг/м. В табл. 2 показаны механические характеристики применяемых материалов: модуль упругости Е; 
предел прочности при сдвиге С; предел прочности при растяжении У. 

Таблица 2 


Механические свойства СЕКР и алюминия 6061 


Параметры в направлениях СЕВР Алюминий 6061 
координатных осей, ГПа 

Ех 128 7070 

Еу 10,5 70 

Сху 5 70 

Су 5. 

Стх Э 

Уху 0,27 0,3 

У 0,4 0,3 

Ух 0,02 0,3 


Создана модель с использованием девяти слоёв КОМПОЗИТНОГО материала С различными углами наклона волокон: 
[+25° /—25° |, [+30° /—30° |, [+35° /-35° ]‚ [ +40° /—40? |, [+45° /—45° |, [+50° /-50° |, [+55°/—55° |. 


На рис. 4 показана последовательность слоев с наклоном волокон под углом +45°, реализованная в 
программе АМ5У5. 


ГАУЕВ 5ТАСКИМС 


ЕГЕМ - 0 
ЕСТ - 1 
ТАУЕВ$ * 
ТОТАЕ .- 0 
$НОММ 1 ТО 
ЕВОМ 0 
Г ауегз Ма{епа!$ 
1 1 
Тпеге 5 : 
45 3 2 
-45 4 | 2 
45 5 
4 ИЖ: 
45 Я 2 
-45 А 8 2 
45 9 


-45 
45 


Рис. 4. Последовательность слоев с наклоном волокон под углом + 45° 


Машиностроение и машиноведение 


111 


Бир://уези-Чопза.ги 


112 


'А4уапсеа Епэтеетия Кезеагсй 2022. Т. 22, № 2. С. 107—115. 155М 2687—1653 


Цифровое моделирование. Расчет напряжений и перемещений. Работа резервуара высокого 
давления проанализирована с использованием критериев разрушения по теории Тза1-У\У‘иа. Для расчета 
максимальных напряжений и перемещений применялось внутреннее рабочее давление 35 МПа. 
Наибольшее напряжение резервуара наблюдалось при угле наклона кольцевых волокон — 0° и различных 
углах наклона спиральных волокон: 


[+25° /—25° |, [+30°/-30° |, [ +35°/-35° |, [+40° / 40? |, [+45° /—45° |, [+50° /—50° |, [ +55° /—55° |. 


На рис. 5 показано распределение эквивалентных напряжений и перемещений в резервуаре при угле 
наклона волокон + 45°. В таблицах 3, 4 приведены максимальные и минимальные значения перемещений и 
эквивалентных напряжений в направлениях координатных осей Х, У, 7. 


ч%7С | 4 ч 
МОБАЕ $ОШТЮМ ЛМ _ с | МОРАЕ ЗОЦ/ЮМ лм _ 
| 
ЗТЕР=1 313 2014 | ЗТЕР=1 332014 
ЗИВ=1 13:01:07 $иВ=1 13:00:15 
ИМЕ=1 ТИМЕ=1 
и$им ЗЕОУ — (АМС) 
(А\б) ОМХ=6.473 
8$У$=0 $ММ=11.886 
ОМХ=6.473 $МХ=463.756 
$МХ=6.473 
о 1.439 2.877 4,316 5754 11.886 112.302 212.717 313.133 413.549 
719265 2.158 3.596 5.035 6.473 62.094 162.51 262.925 363.341 463.756 
а) 6) 


Рис. 5. Эквивалентное распределение в резервуаре: & — напряжений; б — перемещений 


Таблица 3 
Максимальные и минимальные значения перемещений 
в направлениях осей Х, У, 7 
Направление деформации Минимум, мм Максимум, мм 
суммарная векторная 0 6,473 
по оси Х —3,914 3,914 
по оси У — 6,473 0 
по оси 7 —3,906 3,909 
Таблица 4 
Максимальные и минимальные значения напряжений в направлении осей Х, У, 7 
Направление деформации Минимум, МПа Максимум, МПа 
суммарная векторная 11,886 463,756 
по оси Х — 81,852 528,485 
по оси У —155,587 341,607 
по оси 7 — 80,307 530,07 


Определение ожидаемого разрушающего давления и оптимального угла наклона волокон. Нагружение 
резервуара осуществлялось постепенно увеличивающимся внутренним давлением, начиная с его рабочего 
значения — 35 МПа. Затем полученное расчетное максимальное напряжение внутреннего давления сравнивали 


с конечным допускаемым напряжением по условию: о, < 0, 


и ? 


где: о, —Щ_ расчетное максимальное 


тах 
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напряжение; [в —_ конечное допускаемое напряжение. Допускаемое напряжение ДЛЯ резервуаров высокого 
давления из углепластика составляет 1210 МПа. 
Модель цилиндрического резервуара рассчитывалась ДЛЯ двух ВИДОВ путей намотки волокон: кольцевого и 
спирального, а также при различных углах наклона волокон: 
[+25° /—25° |, [+30°/-30° |, [ +35°/-35° |, [+40° /—40° |, [+45° /—45° |, [+50° /—50° |, [ +55° /—55° |. 


При кольцевом пути намотки угол наклона волокон к оси цилиндра составляет 0°. На рис. 6 показана 
конечно-элементная модель цилиндрического резервуара высокого давления. 


ЕЪЕМЕМТ 1 


Рис. 6. Конечно-элементная модель цилиндрического резервуара высокого давления 


Результаты исследования. На рис. 7 показано изменение максимального напряжения, начиная со значения 
угла от 0°и до 45° оно уменышается, а затем до значения 65° возрастает. Установлено, что максимальное 
напряжение в цилиндрическом резервуаре высокого давления, армированного углеродными волокнами, было 
наименьшим при угле наклона волокон + 45°. 
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Угол наклона волокна 
Рис. 7. Зависимость максимального напряжения от угла наклона волокон 


На рис. 8 представлены результаты анализа величины давления для различных углов наклона волокон, 
которые показывают, что, начиная со значения угла 0° и до 45° давление увеличивается, а затем до значения 
угла 65° уменьшается. Наибольшее давление, которое может выдержать резервуар, армированный 
углеволокном без разрушения, составляет 207 Мпа. Это давление имеет место при угле наклона волокон + 45°, 
т. о. оно и является ожидаемым значением разрушающего давления для резервуара. 
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Угол наклона волокна 
Рис. 8. Зависимость давления от угла наклона волокон 
На рис.9 представлено распределение напряжений в резервуаре с углом наклона волокон + 45°, причем 


максимальное напряжение при приложении давления составляет 1213 МПа и оно выше, чем допускаемое 
напряжение для резервуаров высокого давления из углепластика, которое составляет 1210 МПа. 


—— 
МОБАЕ $ОШТЮМ лм _ 


ЗТЕР=1 а 
$ИВ=1 

ИМЕ=1 

$ТМи! (АМС) 

ТОР 

В3У5=0 
ОМХ=25.255 
$ММ=20.52 
$МХ=1213 


30.62 666.86 818.616 1081 
518 687.183 949.847 1213 


Рис. 9. Распределение напряжений в резервуаре при нагнетании давления и угле наклона волокон + 45° 


Обсуждение и заключения. Конечно-элементная модель резервуара, армированного углеродным волокном, 
спроектирована и исследована с использованием программы АМЗУ$ с применением в процессе моделирования 
элемента ЗНЕГТ, 99. Создано несколько моделей с различными углами наклона волокон. Максимальное 
напряжение и разрушающее давление для резервуара высокого давления рассчитывались с использованием 
критериев Тза1-\и. Установлено, что максимальное напряжение оказалась наименьшим при угле наклона 
волокон + 45°, а максимально возможное разрушающее давление в этих же условиях составило 207 МПа. Это 
означает, что оптимальный угол наклона волокон для обеспечения безопасной работы резервуара высокого 
давления составляет + 45°, а резервуар из углепластика, рассчитанный с тем же углом намотки волокон, имеет 


максимальную прочность. 
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Аннотация 

Введение. Проанализированы отечественные и зарубежные работы о проблемах фрезерования 
сложнопрофильных поверхностей сфероцилиндрическим концевым инструментом. Рассмотрены методы 
контроля качества поверхности и пути обеспечения амплитудных параметров шероховатости, которые 
базируются на данных исследований и натурных экспериментов. Представлены теоретические положения по 
определению сил резания и результаты виброакустической диагностики. 

Материалы и методы. Использовались методы корреляционного анализа, сравнения и обобщения 
результатов. Расчетные данные находились при различных углах наклона инструмента с учетом мгновенных 
сил резания и фиксировались в диапазоне значений переменной подачи на зуб (К) и угла наклона 
поверхности (7). Нашли экспериментальное подтверждение сведения виброакустической диагностики и 
теоретические данные представленной модели при различных углах наклона инструмента. Следовательно, 
такие методы можно использовать для прогнозирования параметров шероховатости поверхности. 

Результаты исследования. Выявлена взаимосвязь между силами резания, углом наклона инструмента и 
данными виброакустической диагностики. Сформулирована модель силы резания и перемещений инструмента с 
учетом наклона поверхности. Определен оптимальный диапазон углов наклона инструмента к обрабатываемой 
поверхности, который обеспечивает минимальные значения амплитудных параметров шероховатости. 
Полученные эмпирическим путем звуковые колебания, представленные в спектральном и волновом виде, хорошо 
согласуются с данными из других источников, что позволило сделать вывод о целесообразности прогнозирования 
и мониторинга параметров шероховатости в реальном времени за счет акустики. 

Обсуждение и заключения. Установлено, что рост сил в направлении а.(Х) и (У) наблюдается при у > 40°. Это 
объясняется распределением составляющих силы резания вдоль режущей кромки и зависит от наклона 
поверхности. Уменьшение амплитудных параметров происходит при увеличении угла от 10 до 40 градусов. 
Найденные взаимосвязи силового анализа, направления обработки и виброакустической диагностики 
подтверждают корректность использования виброакустической диагностики для прогнозирования 
шероховатости поверхности. Акустическая диагностика вне зависимости от компоновки технологического 
оборудования позволяет быстро наладить звуковое устройство и оценить влияние режимов резания на 
параметры шероховатости. 


Ключевые слова: шероховатость, сфероцилиндрическая фреза, виброакустическая диагностика, фрезерование, 
наклон инструмента, вибрация, силы резания. 
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Введение. Специалистам недостаточно программных инструментов, которые помогли бы предсказать 
значение силы резания, амплитуды вибраций и т. п. при программировании станков с числовым программным 
управлением (ЧПУ). В результате приходится подбирать скорость подачи на основе опыта, устанавливая 
постоянную скорость подачи для всей обрабатываемой поверхности сложной формы. Это может снизить 
эффективность обработки и качество поверхностей. Для решения данной проблемы необходимо построить 
высокоточную и надежную модель прогнозирования шероховатости поверхности. 

В настоящее время при таком прогнозировании используют два основных подхода. В первом случае 
фотоэлектрические датчики транслируют данные о шероховатости онлайн. Эффективность этого способа 
достаточно низкая т. к. в процессе резания образуется стружка, а в область регистрируемого микрорельефа 


подается смазочно-охлаждающая жидкость (СОЖ). Ко второй группе можно отнести методы контроля с 
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применением искусственного интеллекта, а также теоретическое моделирование и моделирование 
эмпирического регрессионного анализа. В практике наиболее распространены методы из второй группы. 

Авторы работы [1] исследовали фрезерование сфероцилиндрическим инструментом аэрокосмических 
компонентов из алюминиевого сплава [Мб и разработали модель с использованием структур АКМАХ (англ. 
ашотеотезуе тоушс ауегазе \И ехозепои$ шриё5 — авторегрессионная скользящая средняя с экзогенными 
входными данными) для прогнозирования шероховатости поверхности при различных режимах обработки. 
В [2] представлена математическая модель траектории движения сфероцилиндрического инструмента для 
обработки грибовидных лопаток паровой турбины с целью повышения эффективности обработки и качества 
поверхности. Авторы [3] создали аналитическую модель процесса фрезерования сфероцилиндрическим 
инструментом и рассмотрели ситуацию возникновения вибрации для прогнозирования пределов устойчивости 
инструмента при разных углах наклона и скорости вращения шпинделя. В [4] исследуется обработка 
крыльчатки на пятиосевом станке с ЧПУ при использовании сфероцилиндрического инструмента. Траекторию 
движения инструмента и оптимальные режимы резания пришлось разрабатывать из-за сложной криволинейной 
поверхности и близкого расположения лопастей колеса, так как в этом случае невозможно задействовать 
стандартные системы САМ (англ. сотрщег-а14е тапиасвагте — автоматизированная система или ее модуль 
для подготовки управляющих производственных программ). В [5] описана система оптимизации скорости 
подачи для трехосевого фрезерования сфероцилиндрическим инструментом. Использовали возможности 
геометрического моделирования В-гер АСТ: резание под углом преобразовали в ортогональную модель для 
различных инструментов и заготовок. Авторы [6] создали математическую модель постоянной силы резания [7] 
для обработки сложной криволинейной поверхности сфероцилиндрическим инструментом на станке с ЧПУ за 
счет поддержания постоянной скорости резания при различных углах наклона инструмента. В [8] представлена 
динамическая модель процесса фрезерования сфероцилиндрическим инструментом, включающая режимы 
резания и износ кромки инструмента. Авторы [9] разработали систему мониторинга в реальном времени с 
компенсацией погрешностей для повышения точности при производстве деталей свободной формы. В [10] 
исследуется энергопотребление при различных стратегиях резания с целью оценки эффективности каждой 
стратегии резания. Авторы этой работы установили, что оптимальная стратегия позволяет обеспечить заданные 
параметры качества поверхности. 

Многие ученые использовали метод Тагути (Тагучи) для оптимизации параметров обработки [11-13]. 
Некоторые изыскания носят узконаправленный характер и для масштабирования результатов необходима 
дополнительная работа [14—18]. 

Анализ литературных источников выявил отсутствие рекомендаций по использованию метода 
прогнозирования и мониторинга шероховатости в реальном времени. Это тем более касается случаев, когда 
учитываются звуковые колебания в воздушной среде. 

Материалы и методы. Для определения сил резания применяется модель из [7, 19]: 

АР; = Кс, + Кс 4$ 
АЕ, = Ко, + Кс 45}, 


(1) 
ЧЕ; = Ка, + Ки. 


Здесь К,», К» К — краевые удельные коэффициенты, Ку, К», К — удельные коэффициенты сдвига, 4; — 
бесконечно малая длина режущей кромки /[; [20], 5.; — площадь поперечного сечения срезаемого сегмента. 

В этой модели для задания результирующей силы, действующей на 1-й бесконечно малый отрезок режущей 
кромки, используется множество криволинейных систем координат, перпендикулярных касательным к 
сферической поверхности. 


Геометрия сфероцилиндрической концевой фрезы схематично показана на рис. 1. 
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Рис. 1. Схема для определения: а — сил резания сфероцилиндрической фрезы; 
6 — срезаемого сегмента; в — площади поперечного сечения приведенной фигуры [16]. 


На схеме указаны угол наклона А, фрезы со сфероцилиндрической исходной поверхностью инструмента 
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вдоль режущей кромки. Здесь же видим: 

— тангенсальную 4Ё,, радиальную АЁ,, осевую АЁ.; силы резания, действующие на /-Й зуб, и координаты 
инструмента для сфероцилиндрической концевой фрезы (рис. 1 а); 

— хорду [Г2‹и, угол сегмента Ф!, длину дуги („с (рис. 1 6); 

— площадь поперечного сечения (рис. 1 6, в). 
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где Г2.„ — длина второй хорды, мм; Й’„„‚ — средняя толщина стружки, мм. 


г 3 
с = 2 и 
№, = №: ЗА. (4) 
Длина первой хорды Геи, длина дуги [„. и угол сегмента круга Фу: 
Е 9 
= 2.5 : -9< . р, (5) 
[ел 26 
Ф, . к 
р = ь 
. (6) 
Ф, = 2агссо$ ( =. 5. 
28 (т) 
Угол А (рис. 1 в) для расчета средней толщины стружки: 
Н 
А = агсап |" 
Исп /2 (8) 
Радиус-вектор г текущей точки исходной инструментальной поверхности: 
Кзшоф со$@ 
КЮзшф зтО 
г = (9) 
Ксо5ф 


1 


Уравнение режущей кромки решается как уравнение кривой, расположенной на исходной 
инструментальной поверхности и пересекающей меридианы под углом /. 90°. Параметрическое уравнение 
вида © =О (ф) описывает некоторую линию на сфере. Косинусы углов между осями системы координат 
ХУ/ сфероцилиндрической фрезы и касательной прямой к кривой О=оО (9) на исходной 


инструментальной поверхности: 
ах с05ф с05@ 4ф — зшф зшО аф 


с0$@ = — = о, 
45 \/ах” + зт’ф ао’ 


АУ созф зшО 4ф — зшф с0о$О аф 


созВ т т ИК 
45 \/а" + яп’ф 407 (10) 
ай — зшф аф 

с0$0 = = - 
45 \/а" + зт’ф ао” 


где 45 — дифференциал дуги режущей кромки. 
Косинусы углов между осями координат и касательной прямой к меридиану (при © = сопз0) на исходной 


инструментальной поверхности представляются в виде зависимостей С. П. Радзевича [21]: 


ах 
с05% = —— = с05Ф с050, 
45 
АУ 
созВ = — = созф шо, (11) 
45 
ай . 
с0$0 = — = - 51$. 
45 
Поставим условие: 6 = /.. Тогда уравнения режущей кромки в параметрической форме: 
ю КсозА, 
спа (^. + С)’ 
_ Кэш, (12) 
спа (А+ С) 


й = ва + С), 


где сй — гиперболический косинус, #1 — гиперболический тангенс. 
Для определения бесконечно малой длины режущей кромки можно воспользоваться выражением, 


предложенным в [7, 16, 19]: 
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и = |" (+) + (7 (У) + (2 (4) а, (13) 
где 
г (у) = Л - (усоел - 1) 
—В (цу сое ^, - 1) сов. (14) 


п’ (у) = 
\л — (у соёХ - Г 


=’ (у) = Х? + У? + (в- 2)? со 


Основная часть 

Детали эксперимента. На подготовительном этапе изучена литература [22—27]. Данные при планировании 
эксперимента описаны ниже. Механообработка проводилась в попутном и встречном направлениях (рис. 2), на 
заготовках со свойствами стали 45 (зарубежные аналоги — С45, 1045). Использовали сфероцилиндрическую 
концевую твердосплавную фрезу диаметром О = 8 мм, с двумя зубьями фирмы бапахмк Соготапе (5 = 2). 
Соотношение вылета инструмента, установленного в патроне, приняли как [/), равным 4. Фрезерный 
инструмент из мелкозернистого карбида вольфрама имел противоизносное покрытие Т1А!М и параметр /. = 30°. 
Подача на зуб К = 0,2 мм/зуб, величина припуска для всех образцов а, = Н = 0,2 мм, боковой шаг а, = 0,2 мм. 
Для обеспечения равной скорости резания У, при различных углах наклона (у) инструмента к обрабатываемой 
поверхности частота вращения (п) варьировалась в диапазоне от 1478 пут ' до 8000 пит". 


у 
а, 
>. я 
а) 


6) 


Рис. 2. Виды условной поверхности контакта заготовки и инструмента при различных 


значениях углов наклона инструмента в направлении: а — попутном; 6 — встречном 


Угол у выбирался из диапазона от 10° до 50°, поэтому расчет эффективного диаметра О.и, ([л.„) проводился 
по формулам: 


Бсар = 20соз® = 2. В - 00$ (90° - т) 


у > 0° зили (15) 


ЗЕ. 
ед” ао ЕЯ В 


В экспериментах использовали станок с ЧПУ РМ@а РМО 50 (управляющая система НедепВалт). 
Шероховатость поверхности после обработки измеряли при помощи профилометра-профилографа- 
контурографа З$шЁот 18000. В соответствии с действующими стандартами' проводили измерения по 
различным трассам. Для нахождения максимальных значений параметров шероховатости поверхности 
использовали 50 % Гауссовой фильтрации. Виброакустическая диагностика выполнялась согласно ГОСТ Р 
ИСО 7919-3-99° с применением анализатора спектра 7еЧМаб 017-02. В качестве выходной оценки 
эффективности обработки использовали амплитуду А вибраций шпиндельного узла (мм), изменяющуюся со 
временем (1, с). 


1 ГОСТ 27964-88. Измерение параметров шероховатости. Термины и определения. М., 2018. 13 с. 
2? ГОСТР ИСО 7919-3-99. Контроль состояния машин по результатам измерений вибраций на вращающихся валах. М., 2000. 12 с. 
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Результаты исследования. Эксперимент показал, что реальная шероховатость поверхности существенно 
отличается от данных теоретических моделей == результатов кинематико-геометрической проекции 
инструмента на заготовку. В настоящее время растут востребованность в деталях пространственно-сложной 
конфигурации И запросы К качеству поверхности. В этой СВЯЗИ СТОИТ сфокусироваться на требованиях 
предсказуемости результатов механической обработки, выраженных в том числе показателями шероховатости. 

В таблице 1 И 2 представлены результаты эксперимента амплитудных параметров шероховатости, 
описанных в ГОСТР ИСО 42883. 

Таблица 1 
Параметры шероховатости в зависимости от направления обработки сфероцилиндрическим 
инструментом при попутном фрезеровании 


Амплитудные параметры шероховатости, ит 
Угол, ° Ка Ка В Ч Кр Ку 
10 0,457 0,569 2,912 5,000 1,671 1,240 
20 0,597 0,735 3,936 5,159 2,335 1,600 
30 0,449 0,561 3,007 4,960 1,808 1,200 
40 0,479 0,590 3,232 4,239 1,823 1,407 
50 0,787 1,027 6,615 12,68 4,119 2,496 


Параметры шероховатости в зависимости от направления обработки сфероцилиндрическим 
инструментом при встречном фрезеровании 


Таблица 2 


Амплитудные параметры шероховатости, пит 
Угол, ° Ка Ка К № Кр Ку 
10 0,806 1,006 4,776 5,960 2,776 2,000 
20 0,752 0,935 5,207 6,440 2,847 2,359 
30 0,509 0,647 3,727 5,239 2,223 1,503 
40 0,628 0,805 4,832 7,400 2,912 1,919 
50 0,730 0,915 5,111 7,639 3,039 2,072 


Анализируя данные таблиц 1 и 2, отметим, что превалирующее влияние на «высоту» шероховатости 
оказывает направление резания. Количественное значение амплитудных параметров «ниже» при попутном 
фрезеровании. Далее будем говорить только о таком фрезеровании. 

На основании данных таблицы 1 провели корреляционный анализ и составили коррелограмму (рис. 3). 
Выявленные зависимости соотносимы с результатами работ [28, 29]. 


К Ку 


К! 


Ка-————Ка 


— 0,9...1,00) 0,8...0,89 .--- 0,7...0,79 


Рис. 3. Коррелограмма амплитудных параметров шероховатости после фрезерования в попутном направлении 


Уравнения, полученные в ходе корреляционного анализа, помогают отойти от прямого нормирования 
параметров шероховатости. При этом аналитические зависимости способствуют минимизации числа 
контролируемых факторов. Это снижает трудоемкость перехода к мониторингу звуковых колебаний в 


реальном времени. 


3 ГОСТ Р ИСО 4287-2014. Геометрические характеристики изделий (ОРЗ). Структура поверхности. Профильный метод. Термины, 
определения и параметры структуры поверхности. М., 2015. 20 с. 
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В представленном исследовании рассчитаны силы резания (Ех, Ку, Ез) (рис. 4). Видно, что угол наклона 


поверхности у существенно влияет на силы резания. 


Е, Н Е Н 

100 100 
—5 

50 157 50 
Ех 

0 0 

—50 50 

—100 100 ре 
0 0.005 0.01 0.015 0.02 ьс Го) 0,005 0,01 0,015 0,02 
а) 6) 


Рис. 4. Силы резания, рассчитанные для: а — у = 10; 6 — у = 40° 


Полученные данные позволили установить, что силы уменьшаются с увеличением угла наклона 
инструмента от 10° до 40°. Кроме того, сила Ех наиболее чувствительна при изменении угла наклона 
поверхности. Отметим также, что силы резания, рассчитанные по вышеупомянутой методике, хорошо 
согласуются с полученными эмпирическими результатами работ [6, 17, 24, 26, 27, 30]. 

С ростом угла ‘у уменьшается предельная площадь погружения инструмента в заготовку (определяемая 
рабочим углом и активным числом зубьев <.). Следовательно, при чистовом  фрезеровании 
сфероцилиндрическим инструментом могут присутствовать «пульсирующие» силы, поскольку для угла 
наклона поверхности у > 0° число активных зубьев часто меньше единицы (< < 1). 

Чтобы установить <. проанализировали звуковые колебания в ходе резания. Авторы работы [31] 
предложили разделить спектр колебаний рабочих органов на диапазоны: 

— 20...300 Гц — низкочастотный, 

—300...1500 Гц — среднечастотный, 

— 1500 Гц и более — высокочастотный. 


Экспериментально полученный диапазон звуковых колебаний представлен на рис. 5. 


Рис. 5. Экспериментально полученный спектральный вид звуковой волны (сегмент равен 1 сек) 


К низкочастотному диапазону следует отнести колебания установочных элементов и приводов подачи, 
к средним и высоким — частоту инструмента, шпинделя и др. В данной работе с точки зрения 
акустической диагностики определены оптимальные режимы резания, поддерживающие устойчивую 
работу технологической системы. 

Предварительно определили частоту, соответствующую частоте резания (рис. 6). Применили фильтрацию к 
звуковым колебаниям при: 

— перемещении шпиндельного узла (У, 7), 
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— вращении (п), 
— перемещении стола (Х), 
— шуме в лаборатории. 


ЕЕ рии чет ЕЛЕ ТЕТ 


Рип 0 0.2 0.2 0.4 0:5 


Рис. 6. Спектральный вид звуковой волны после фильтрации (сегмент равен 1 сек) 


Стоит отметить некоторый разброс значений (рис. 6), связанный с огромным количеством факторов, 
которые сложно исключить при проведении натурного эксперимента. 

Середина выраженного диапазона частот, представленных на рис. 5, соответствует вращению 1500 об/мин с 
наклоном инструмента 40°. Акустическая диагностика подтвердила, что резание шло практически непрерывно, 
т. е. без ударных «пульсирующих» нагрузок. Последняя ситуация представлена на рис. 7. Здесь выделенное 


время оценки составляет 0,5 с, количество периодов — 25. 


Рис. 7. Фильтрованный фрагмент звуковой волны, которая основана на звуке резания при 7 = 40° (сегмент равен 1/2 сек) 


На рис. 7 представлена неоднозначно интерпретируемая звуковая волна. Однако ее довольно просто 
охарактеризовать, если углубиться в сущность механической обработки фрезерованием. При попутном 
фрезеровании и скруглении режущей кромки пластические деформации металла с остаточными 
напряжениями остаются на обработанной поверхности. Даже при должной подаче СОЖ в поверхностных 
слоях материала возникают температурные деформации и искажают кристаллическую решетку (остаточные 
напряжения Ш рода). Наряду с этим после удаления инструмента растянутые верхние слои металла 
приобретают остаточные напряжения сжатия. Если описывать механообработку со стороны кинематики 
процесса, стоит отметить фактор снятия стружки, а именно величину ©› — угол резания материала с 


наклоненным шпинделем в попутном направлении (рис. 8). 
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Рис. 8. Срезаемый при механообработке сегмент 2: а — вид сверху; б — проекция срезаемого сегмента [16] 


Уточним, что инструмент режет в зоне 2. Полученные данные подтверждаются результатами 
моделирования, но рассмотрены для идеального случая. 

Наряду с указанными факторами теоретические и эмпирические расхождения часто вызваны вибрациями, 
связанными с отклонением и износом инструмента, а также с геометрией зоны резания. На рис. 7 
прослеживается процесс подрезания материала заготовки ]-м зубом инструмента, который наряду с другими 
факторами определяет столь интересную форму полученной волны. 

В дополнение к вышесказанному на рис. 9 представлены зависимости, которые показывают, что при 
уменьшении амплитуды вибраций увеличивается угол наклона от 10° до 40°. 


А, мм 

в [ 
ооо ПОЙ Й ЗАЛ М, Да АЛ. 
ооди5 ГА [АаМАНАЙ АСАИАЮА ДАНИЯ АЕ 
о 


Рис. 9. Амплитуды колебаний в зависимости от угла наклона поверхности: сплошная линия — 10°, пунктирная — 40° 


М.Р. Гимадеев и В. В. Гусляков* проанализировали международную практику профильной оценки 
шероховатости поверхности. Авторы работы [32] выдвинули предложения по переходу на трехмерную оценку. 
В представленном исследовании для составления комплексной оценки измерения шероховатости проводились 
в ортогональной системе координат. 

Анализ полученных профилограмм (рис. 10) позволяет сделать заключение о регулярности микрорельефа 
поверхности после обработки сфероцилиндрической концевой фрезой. При этом нецелесообразно учитывать 
шероховатость инструмента, так как ее проявления незначительны. Некоторое дополнительное увеличение 
шероховатости наблюдается в направлении бокового шага а., что подтверждает увеличение упругих отжатий 
инструмента. Об этом свидетельствуют представленные профилограммы, где периоды профиля (расстояние 


впадин и выступов) соответствуют заданной величине подачи К = 0,2 мм/зуб. 
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Рис. 10. Профилограммы поверхностей после фрезерования с различными углами наклона: а — у: 7 = 10°, направление }; 
б — у= 10°, направление а.; в — 7 = 40°, направление }; г — у = 40°, направление а, 

Фотомикроскопические изображения обработанной поверхности имеют очевидную гексагональную 
форму (рис. 11). Качество обработанной поверхности в направлении бокового шага а, может снижаться из-за 
отжатия режущего инструмента и возрастания амплитуды колебаний технологической системы при 
фрезеровании. Это подтверждают результаты виброакустической диагностики. 
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0,5 мм[ ] 


Рис. 11. Схема фрезерования и фото обработанной поверхности при: а — 7 = 10; 6 —у=40 


Еще раз подчеркнем, что результаты экспериментов соотносятся с данными, полученными многими 
авторами, и согласуются с теориями упругости, формообразования, резания материалов и технологией 
машиностроения. 

Обсуждение и заключения. Изучено влияние угла наклона обрабатываемой поверхности на силы резания, 
вибрации и параметры микрорельефа. Создана модель силы резания и вибрации, включающая кинематико- 
геометрические параметры. Исследование показало, что количественно и качественно силы резания и 
колебания режущего инструмента зависят от угла наклона поверхности. 

Анализ экспериментальных данных позволил выявить при увеличении угла наклона до 30—40° тенденцию к 
снижению амплитудных параметров шероховатости на 10-20%. Качество обработанной поверхности 
снижается из-за упругого отжатия сфероцилиндрического инструмента и возрастания амплитуды колебаний 
технологической системы. 

Установлено, что силы растут в направлении а.(Х) и К(У) при 7 > 40°. Это можно объяснить распределением 
составляющих силы резания вдоль режущей кромки, которое зависит от наклона поверхности. 

Полученные авторами решения задач мониторинга и анализа параметров шероховатости за счет 
акустической диагностики позволяют существенно сократить объем экспериментальных исследований. Такой 
подход дает возможность значительно снизить затраты материальных ресурсов и времени на механическую 
обработку изделий. Результаты экспериментов представляют практический интерес и уточняют сведения о 
формировании микрорельефа, в том числе о параметрах шероховатости при механической обработке резанием. 
Необходима более детальная проработка предложенного метода акустической диагностики. 
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Определение динамических характеристик шестерённого насоса 
по методу вариации нагрузок с использованием специальных 
стендовых систем 

В. И. Санчугов ©, П.Д. Рекадз‹ @® 


Самарский национальный исследовательский университет имени академика С. П. Королева, Российская Федерация, 
г. Самара, ул. Московское шоссе, 34 
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Аннотация 

Введение. Пульсации рабочей среды в гидросистемах станков, системах топливоподачи авиационных 
двигателей и подачи жидкостного ракетного двигателя являются одним из главных факторов возникновения 
вибрации агрегатов, приводящей к усталостному разрушению корпусов элементов, трубопроводов и выходу из 
строя элементов насосов. Целью работы являлись реализация метода определения динамических характеристик 
объёмного насоса с использованием специальных стендовых систем и сравнение результатов расчёта с 
экспериментальными данными. Описаны этапы расчёта динамических характеристик объёмного насоса с 
учётом разработанных ранее специальных стендовых систем на примере шестерённого насоса внешнего 
зацепления производительностью 14 см/об. Показана реализация метода вариации нагрузок В. П. Шорина с 
использованием разработанных авторами специальных стендовых систем с заранее определёнными 
динамическими характеристиками. Излагаются основные этапы методики определения динамических 
характеристик шестеренного насоса. 

Материалы и методы. В работе применены методы спектрального анализа пульсирующего давления. 
Пульсации расхода жидкости на выходе из насоса определялись с помощью импедансного метода, метода 
вариации нагрузок и специальных стендовых систем. 

Результаты исследования. В работе реализована методика определения динамических характеристик 
шестеренного насоса в диапазоне оборотов приводного вала 500-2500 об/мин для четырёх гармонических 
составляющих спектра колебаний в широком диапазоне динамических нагрузок (от инерционного до 
емкостного характера). Представлен анализ стендовых систем, приводящих к расчёту динамических 
характеристик насоса с минимальной погрешностью на основании условия согласования динамической 
нагрузки и источника колебаний. Проведена проверка разработанного подхода к оценке динамических 
характеристик насоса посредством сравнения расчётных и экспериментальных данных пульсаций давления в 
стендовых системах с дросселем, полостью и протяжённым трубопроводом на выходе насоса. 

Обсуждение и заключения. Реализован метод определения динамических характеристик объёмного насоса с 
использованием разработанных авторами специальных стендовых систем. Результаты исследований 
показывают, что исследуемый шестерённый насос может рассматриваться как независимый источник 
колебаний расхода, для которого отклонение собственных динамических характеристик от средних значений не 
превышает 10 % по первой гармонической составляющей. 


Ключевые слова: объёмный насос, метод вариации нагрузок, стендовые системы, источник колебаний расхода, 
динамические характеристики насоса, независимый источник колебаний, пульсации расхода, импеданс. 
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Введение. Основные недостатки насосов заключаются в сопровождаемом их работу акустическом шуме, 
вибрации корпуса, присоединенных трубопроводов и возникающих вследствие этого разрушениях конструкции. 

Традиционный путь совершенствования конструкций насосов, используемый в машиностроении и 
заключающийся в сглаживании неравномерности течения жидкости, затруднен отсутствием средств прямого 
измерения расхода жидкости на выходе насоса. Это ограничивает любую экспериментальную оценку 
эффективности конструктивно-доводочных работ по насосам. Кроме того, это ограничивает аналитическое 
описание насоса как источника колебаний, необходимое разработчикам гидросистем для обеспечения 
требуемого динамического качества самих систем. 

По этой причине исследования, связанные с моделированием динамических процессов в гидромеханических 
системах [1,2] и применением новых подходов к решению задач динамики гидравлических систем [3, 4] на 
основе физических принципов, имеют актуальность. 
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Материалы и методы. Представлена реализация метода вариации нагрузок В.П. Шорина [4], 
заключающаяся в том, что определяются пульсации давления за насосом, взаимодействующим с несколькими 
стендовыми системами, динамические характеристики которых рассчитаны заранее. Далее обработанные 
амплитудно-частотные характеристики давления с использованием импеданса стендовых систем 
пересчитываются в собственные динамические характеристики насоса — пульсационную производительность 
и импеданс. Развитие подхода состоит в том, что предложены специальные стендовые системы, учитывающие 
соединительные штуцеры, переходники и внутренние каналы насоса, что позволяет точнее определять 
динамические характеристики насосов, а также осуществлять проверку применимости используемой модели 
насоса в виде эквивалентного источника колебаний. 

Результаты исследования 

Методика определения динамических характеристик шестеренного насоса 

В Самарском университете разработана расчетно-экспериментальная методика определения динамических 
характеристик шестеренного насоса (пульсационной производительности и импеданса), которая базируется на 
исследованиях Л. Тевенена, Э. Нортона и В. П. Шорина [5, 6]. Разработанная методика включает в себя 
следующие основные этапы: создание моделей и расчет динамических характеристик специальных стендовых 
систем; экспериментальное определение пульсаций давления на выходе насоса; расчет спектра возбуждаемых 
колебаний в различных стендовых системах; определение аппроксимирующих зависимостей для отдельных 
гармонических составляющих; расчёт собственных динамических характеристик насоса (пульсационной 
производительности и импеданса); оценка независимости (стабильности) динамических характеристик насоса 
от характеристик стендовых систем; проверка применимости используемой модели насоса в виде 
эквивалентного источника колебаний расхода (или давления). Далее приведен пример реализации 
предложенной методики с использованием шестерённого насоса внешнего зацепления удельной 
производительностью 14 см“/об. 

Модели и расчет динамических характеристик специальных стендовых систем 

Рассмотрены базовые варианты стендовых систем на выходе жидкости при испытаниях насосов. В 
сосредоточенных параметрах стендовая система может реализовывать активную расходно-зависимую, 
инерционную и упругую нагрузки. В распределённых параметрах стендовая система может быть реализована в 
виде протяжённого цилиндрического трубопровода, реализующего активную частотно-независимую 
неотражающую нагрузку. 

Модели насоса в виде эквивалентных источников колебаний расхода (ЭИКР) или давления (ЭИКД) 
определяются с использованием источника постоянного расхода О,, идеального источника колебаний давления 


р (или расхода 4) и внутреннего импеданса 7, [4]. На рис.1 представлены модели насоса в виде 


эквивалентного источника колебаний давления (и расхода) с присоединенной системой. 


стендовоя 
НаСОС | система насос | система 
а) 6) 


Рис. 1. Модели насоса в виде эквивалентного источника колебаний: а — давления; 6— расхода с присоединенной 


системой: О, 5 Б, = средний расход и давление; р„,4и— собственная пульсационная производительность насоса по 


давлению и расходу, 7, Иа — импеданс насоса и стендовой системы (адаптировано из [4]) 


и? 


При этом задачей выбора стендовой системы является обеспечение таких параметров системы, при которых 
реализуются наилучшие условия высокой точности измерений. Задачу выбора параметров присоединенной 
системы необходимо проводить при изменении импеданса в диапазоне от -© до +00. 

Принципиальные схемы стендовых систем представлены на рис. 2. 
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Рис. 2. Принципиальные схемы стендовых систем для определения собственных динамических характеристик насосов: 
а — сактивной дроссельной нагрузкой; 6 — с активной неотражающей нагрузкой; в — емкостной нагрузкой: 
1 — испытуемый насос; 2 — нагрузочный дроссель; 3 — расходомер; 4 — датчик давления; 5, 6 — соединительные 
штуцеры; 7 — магистраль; 8 — протяжённый трубопровод; 9 — полость 


Расчётные модели импедансов специальных стендовых систем: 
— с активной дроссельной нагрузкой (дроссель): 
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тр 
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— с активной распределенной нагрузкой (протяжённый трубопровод): 
И & ты _— 2 _ 


тр 


— с емкостной нагрузкой (полость): 


2 (3) 
те) 1-1 
С тех (О, Л ) = (ке? (0. ) -+ Кей ив (0, ) + Кей. (0. ) + р (0. )} + ЕЕ ы 
прС 
в (2 сы 
ве све „кеты ве, „(ль (6, 
В формулах обозначены: 7.» (0,) —Щ импеданс дросселя; Кей’ (©) —щ действительная часть импеданса 


трубопровода, включающего соединительную магистраль, штуцеры, внутренние каналы в насосе и дросселе; 
Кей ‚Кей _ ‚ Кей —Й действительная часть импеданса штуцеров и внутреннего канала в 
`лам о `турб о 'м.с. о 


насосе и полости при ламинарном, турбулентном течении жидкости за насосом и местных сопротивлениях 
соответственно; {, — частота К-й гармонической составляющей спектра колебаний за насосом; р — 


Плотность рабочей среды; 1; —Щ длина трубопровода, включающая длину магистрали 6, соединительных 
штуцеров 5и7и внутренних каналов насоса и дросселя (не показаны на рис. 7); ой — площадь поперечного 


сечения трубопровода; Г, — «инерционность» соединительных штуцеров и внутреннего канала в насосе и 
полости; С — упругость полости; К — номер гармонической составляющей спектра колебаний; | — мнимая 


единица (] = -1 ). 
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С использованием стендовой системы с дросселем на выходе насоса реализованы три инерционно-активные 
динамические нагрузки, в которых изменялось активное сопротивление изменением площади перекрытия 
дросселя Но: 

— для стендовой системы №1 площадь перекрытия Ео!=9,290х 10° м2, 

— для стендовой системы №2 площадь перекрытия Ео›=7,299х 10° м7, 

— для стендовой системы №3 площадь перекрытия Е‹з=2,247х 103 м?. 

С использованием стендовой системы с полостью на выходе насоса объёмом 2 л реализована инерционно- 
упругая динамическая нагрузка. 

С использованием стендовой системы с протяжённым трубопроводом длиной 106 м реализована активная 
частотно-независимая нагрузка. 

По моделям стендовых систем по формулам (1)—(3) рассчитаны зависимости импеданса (модуля и фазы) от 
частоты. В таблице 1 для 1-й гармоники приведены функции регрессии модуля импеданса стендовой системы 
7(Р и фаз В;(Е) для анализируемых стендовых систем. 


Таблица 1 
Функции регрессии динамических характеристик стендовых систем по 1-й гармонике 
№ Модуль импеданса 
стендовой Тип стендовой стендовой системы Фаза импеданса стендовой системы 
системы по системы кг В; (В, рад 
Я: 
рисунку 2 м‘с 
а— №1 1,412-107+ 2,201-107 7,066-10# + 0,199 
С дросселем на - 9 т 
а— №2 НЕОО-НАСОСЬ 1,610-10 + 4,258:10 6,906:10"Е+0,09039 
а— №3 4,980-107Ё + 1,503-10? 2,723-103#+ 0,251 
С полостью на 3,308-10“Ё — 1,416.10 + 4,234-10° Г +1,495.10? 
б- №4 9 атс ( - с) 
выходе насоса 1,873.10 1,430.10} +2,930.10 
С протяжённым 
в- №5 трубопроводом на 1,897-1010 0 
выходе насоса 


Экспериментальное определение пульсации давления на выходе насоса 

Исходной информацией для определения динамических характеристик насоса являлось измерение давления 
на выходе насоса. Пример осциллограммы и спектра пульсаций давления на выходе из насоса в стендовой 
системе с активной распределенной нагрузкой представлены на рис. 3. Режим работы насоса при этом: среднее 
давление РД, = 23,8 МПа ‚ частота вращения приводного вала п =100006 / мин. 
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Амплитуда пульсаций давления, МПа 
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Частота составляющих спектра колебаний, Гц 
6) 


Рис. 3. Давление на выходе насоса: а — осциллограмма; 6 — спектрограмма 


Амплитудно-частотная характеристика пульсаций давления на выходе насоса в диапазоне оборотов 
приводного вала насоса п=500-1250 об/мин представлена на рис. 4. Сплошной линией нанесены функции 
регрессии, построенные в среде Мсгобой Ехсе! при помощи инструмента «Линия тренда». 
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Рис. 4. Амплитудно-частотная характеристика пульсаций давления на выходе насоса по отдельным гармоническим 
составляющим (точками обозначены данные на: ®@ — первой, ® — второй, ® — третьей, ® — четвёртой гармониках, 
линиями обозначены функции регрессии) 


Анализ осциллограммы, спектрограммы и амплитудно-частотной характеристики пульсаций давления 
показал следующее: тип процесса полигармонический установившийся; число регистрируемых гармоник — 8, 
анализируемых' — 4; диапазон амплитуд колебаний давления 0,28—0,52 МПа; частотный диапазон приводного 
вала 500-1250 об/мин, анализируемых гармонических составляющих спектра 83-833 Гц. 

Для гармонического анализа периодических сигналов использовалось разложение в ряд Фурье [7]. Для 
разложения в ряд Фурье использованы участки осциллограммы давления, измеряемые относительно опорного 
сигнала на роторной частоте, т.е. за полный оборот ротора извлекался один участок осциллограммы. Опорный 
сигнал регистрировался при фиксированном положении шестерни, выставляемом при начале записи 
осциллограммы давления для пяти стендовых систем. Таким образом определялись амплитудно-частотные 


(АЧХ) А, ( @) и фазо-частотные (ФЧХ) ф, (<) характеристики. 


1 В работе принято, что крайняя высшая анализируемая гармоника по амплитуде не должна быть более 5 % от амплитуды первой 
гармоники. 
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Расчет спектра возбуждаемых колебаний в различных стендовых системах и определение 
аппроксимирующих зависимостей для отдельных гармонических составляющих колебаний давления 

На основе экспериментального определения давления в пяти стендовых системах: с дросселем (№1-3), с 
полостью (№4) и протяжённым трубопроводом (№5) на выходе насоса рассчитаны спектры возбуждаемых 
колебаний. В таблице 2 приведены функции регрессии амплитуд А,(Р) и начальных фаз пульсаций давления для 
1-й гармоники. 


Таблица 2 
Функции регрессии амплитуд и начальных фаз пульсаций давления 
№ стендовой . Амплитуда пульсаций давления Начальная фаза р 
Тип стендовой ы Частотный 
системы по кг пульсаций давления 
системы Ар (в, —— диапазон 
рисунку 2 мс ф;: (В, рад 
а— №1 1,003-107% —0,543 
С дросселем на Е 
а— №2 выходе насоса 1,243.10 Е —0,280 
а— №3 6,129.10 —5,941 
С 
6 №4 о 1,343-10°°Р — 4,298-104#+ 3,913-107 -5,373 83—417 
выходе насоса 
С протяжённым 
в №5 трубопроводом на 2,694-107% —0,658 
выходе насоса 


Расчёт динамических характеристик насоса 


Для расчета динамических характеристик насоса (импеданса и пульсационной производительности 
источника) необходимо провести эксперимент как минимум с двумя динамическими нагрузками”. При этом 
импеданс насоса по уточнённой авторами модели эквивалентного источника колебаний расхода с 


использованием аналитических зависимостей р, (/ ) и 7, (1) из таблиц 1 и 2 рассчитывается по формуле: 


1 
=2, а. (4) 

и И 

ры 

р. 2, 

Г 2 2 

где 7, — импеданс насоса (7, =|7, е*) — модуль импеданс насоса (|2, =\/Кейи +Пий и ), Вей „1т7, — 
действительная и мнимая часть импеданс насоса, ф — аргумент импеданса насоса, 7,,/, — входной 
импеданс динамических нагрузок № 1-2 (= е^, 2, =|2.|е*), 8,8, — аргументы имнедансов 
динамических нагрузок № 1-2, р,р, —щ пульсации давления за насосом в стендовой системе № 1-2 ( 
р = Ане" ‚Р› =А зе ), Ат, Ар» — амплитуды пульсаций расхода за насосом в стендовой системе № 1-2, 


ф,ф, — начальные фазы пульсаций давления динамических нагрузок № 1-2. 

Затем с использованием аналитических зависимостей из таблиц 1,2 определяется пульсационная 
производительность насоса по уточнённой авторами модели эквивалентного источника колебаний (переменная 
составляющая объёмного расхода а, и давления р, ): 


Ее 1-2 
о (5) 
“ 2, | и И 
Р Р> 2 
где 4,‚р, — переменная составляющая объёмного расхода (а, = Ане) и давления (р, = Аше” ”_) насоса, 
А„.Ам — амплитуда пульсаций расхода и давления насоса, ф,ф, — аргумент переменной составляющей 


объёмного расхода и давления насоса. 
Из рассчитанных характеристик по модели ЭИКР (а, ( Г) и ЭИКД (р, ( У )) по формуле (5) только одна 


будет отличаться стабильностью, т.к. эквивалентный источник колебаний имеет независимую динамическую 
характеристику только по одному из параметров (пульсирующий расход или давление). 


2 И ы 
Для достаточного объёма выборки необходимо использование хотя бы четырёх динамических нагрузок 2,, формирующих шесть 


неповторяющихся (в том числе 1 = ] ) комбинаций: бо, бл ба» Из, бд, Иа. 
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Расчет импеданса 7, и переменной составляющей объёмного расхода 4, насоса по формулам (4)-—(5) 


проведён для пяти динамических нагрузок, т.е. с использованием десяти различных комбинаций этих нагрузок: 
«1_2», «1_З», «2_З», «14», «2_4», «З_4», «1_5», «2_5», «З_5», «4_5» (см. таблицу 1). 

Для дальнейшего расчёта и анализа из десяти кривых амплитуд пульсаций расхода получена характеристика 
амплитуд пульсаций расхода насоса А„. Значения амплитуд А„ (на каждой частоте) рассчитывались по 


формуле из источника [8]: 


А„=А„ АА, 
где А„ — средняя амплитуда пульсирующего расхода, ДА„, — доверительный интервал амплитуды пульсаций 
19 
расхода насоса (ДА =-——), о — среднеквадратическое отклонение амплитуды пульсирующего расхода от 


Ут 
среднего значения, и — число комбинаций динамических нагрузок (т = 10), { — коэффициент Стъюдента ( 


1=2,262 , рассчитан для доверительной вероятности 0,95). 

Доверительный интервал амплитуд ДА„ формируется из погрешностей (инструментальной и методической) 
и степени зависимости динамической характеристики насоса А, (} ) от стендовой системы. Инструментальная 
погрешность состоит из погрешности определения давления и фаз, методическая погрешность состоит из 
погрешности неучтённого влияния динамических свойств приводной механической системы. 

Фазы пульсаций расхода ф, и их отклонения рассчитывались аналогичным образом. При расчёте фаз 
комбинации «1_2» и «1! _5» исключены по причине вырождения переменной составляющей объёмного расхода 
насоса 4, при использовании близких значений характеристик на одной частоте: ф =фир =р,. 

Результаты расчёта собственной динамической характеристики насоса в виде пульсационной 
производительности и импеданса представлены на рис. 5. Фазы для удобства представлена в градусах. 
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2 
= 156840 
м 
= 14Е+10 
8 1,2Е+10 
о 
[3 
= 10Е+0 
8 
8 8.08+09 
| 
8 6.0Е+09 
2 
= 408-409 
Е 
2 2,0Е+09 
Я 10Е+00 
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г) 
Рис. 5. Динамические характеристики насоса: а — амплитуда пульсаций расхода; 6 — фаза пульсаций расхода; 


в — модуль импеданса; г — фаза импеданса (сплошными линиями обозначены данные по: ча» — первой; 
чннинию — ВТОРОЙ; чине —— Третьей; чаю —— четвёртой гармониках) 


Амплитуда пульсаций расхода на 1-й гармонике имеет монотонно возрастающий характер. С 
увеличением номера гармоники наклон кривой амплитуд пульсаций расхода снижается до практически 
постоянного значения к 4-Й гармонике. Значения амплитуд А, на высших гармониках на порядок меньше, 


чем на 1-й гармонике. 
Модуль импеданса насоса на 1-Й гармонике имеет монотонно возрастающий характер. 
Для оценки реализованных динамических нагрузок использован критерий согласования модуля импеданса 


|2, источника колебаний и входного модуля импеданса динамической нагрузки |2 [99, 10]. Модули 


ст.с 


импедансов должны быть сравнимы (|7,|=| 


ст.с 


), чтобы количественные значения динамической 


характеристики насоса определялись с минимальной погрешностью. Для сравнения использовано отношение 


и и 


7 й При =1 динамическая нагрузка считается согласованной с источником. 


ст.с. ст.с. 


Ввиду различия форм кривых модулей импеданса насоса и динамических нагрузок в диапазоне частот 83— 


и 


417 Гц рационально использовать среднее отношение ‚› Которое позволит учитывать вклад большей 


Ст. С- 1 и Тнах 


части частотного диапазона. Т.е. отношение по 1-й гармонике составит: 
нах 

— для активной дроссельной нагрузки №1 — 1,16; 

— для активной дроссельной нагрузки №2 — 0,8; 

— для активной дроссельной нагрузки №3 — 2,65; 


— для емкостной нагрузки №4 — 4,64; 
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— для активной распределенной нагрузки №5 — 0,35. 

Т. е. динамические нагрузки стендовых систем №1,2 являются наиболее согласованными с источником 
колебаний. Таким образом, динамическая характеристика насоса, рассчитанная с использованием 
динамических нагрузок стендовых систем №1,2 («1 3», «23», «1 4», «24», «|5», «2 5») определена с 
минимальной погрешностью. 


Оценка стабильности динамических характеристик насоса от характеристик стендовых систем 


Стабильность динамических характеристик Ам и ф, ри, 


и ф,) насоса оценивалась с использованием 


инструментов математической статистики [8]. Так стабильность характеристики Ам оценена через расчёт 


отношения доверительного интервала [АА на частоте Тк среднему значению характеристики | на этой 
3; 
частоте Полученное относительное отклонение А, не должно превышать заданное отклонение Азот : 
ДА 
— Чи 
д. т А — Азот ? 
4иМе 
где Ашме среднее медианное значение амплитуды пульсаций расхода насоса; Аа задаваемое 


относительное отклонение амплитуды пульсаций расхода насоса от среднего значения. 


При относительном отклонении Д, —>0 (когда [АА 


—0) стабильность характеристики А„ максимальна. 


Исходя из опыта других исследователей, предпочтительно выбирать величину А, в пределах 5-30 % [11-17]. 


допа 


Для одночисловой оценки динамических характеристик насоса рационально использовать величину 
относительного отклонения. Расчёт такой величины возможен по оценке среднего медианного значения 


ДА, 
относительных отклонений в диапазоне частот каждой гармоники ( 1... Ли»): Ламе = = . Величину 
А. Лт --“Лтах 
АА х 
Ам рационально использовать, если — < 0,3 ‚ т.е. когда доверительный интервал ДА» не превышает 


44 У, 
30 % от средней амплитуды пульсирующего А„. 


Результаты оценки стабильности динамических характеристик насоса от характеристик стендовых систем“ 
составляют: по 1-й гармонике — 2%, по 2-й — 8%, по 3-й — 15%, по 4-й — 48%. Результаты расчёта 
динамических характеристик насоса показали, что по четырём гармоникам в диапазоне частот 500-2500 об/мин 
исследуемый шестерённый насос целесообразно рассматривать по модели эквивалентного источника 
колебаний расхода. 

Проверка применимости используемой модели насоса в виде эквивалентного источника колебаний расхода 

Для проверки разработанного подхода оценки динамических характеристик насоса было выполнено 
расчётное определение пульсаций давления в стендовых системах с дросселем, полостью и протяжённым 
трубопроводом на выходе насоса. С использованием найденных динамических характеристик насоса 
рассчитаны амплитуды пульсаций давления в анализируемых стендовых системах, которые сравнивались с 


экспериментальными. Ниже представлены результаты расчёта амплитуд А, ( 1) в стендовой системе с 
протяжённым трубопроводом. Амплитуды пульсаций давления р, рассчитаны в соответствии с моделью ЭИКР 
по формуле: 
4,:7,-7, (6) 
7,47, 
Результаты расчёта амплитуд пульсаций давления за насосом приведены на рис. 6. Пунктирной линией 
обозначены амплитуды, полученные экспериментально, а сплошной линией — рассчитанные по формуле (6). 


р; = 


ь Стабильность других собственных характеристик насоса оценивается аналогично. 
ь Необходимо учитывать методическую погрешность подхода и инструментальную погрешность определения давления. 
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Рис. 6. Амплитудно-частотная характеристика давления для стендовой системы с протяжённым трубопроводом длиной 
[=106 м (сплошными линиями обозначены расчётные данные по: ча» — первой; чение — второй; 
чниииинию — ТрстьсйЙ; чин — четвёртой гармониках; пунктирными линиями обозначены экспериментальные данные) 


Оценка сходимости результатов проверочного расчёта проведена по формуле отклонения разницы 
экспериментальных и расчётных значений от экспериментальных данных. Результаты проверочного расчёта для 
стендовых систем с дросселем, полостью и протяжённым трубопроводом на выходе насоса показали, что по 1-й 
гармонике относительное отклонение от экспериментальных данных для рассмотренных стендовых систем не 
превышает 10%, минимальное относительное отклонение характерно стендовой системе с протяжённым 
трубопроводом на выходе насоса — 2%. Амплитуды пульсаций давления 2-й-4-й гармоник находятся в пределах 
коридора, рассчитанного по формуле А, =А,-А, (где А, — средняя амплитуда пульсирующего давления). 

Обсуждение и заключения. Реализована методика расчёта динамических характеристик насоса в виде 
эквивалентного источника колебаний расхода жидкости на примере шестерённого насоса производительностью 
14 см*/об с использованием разработанных авторами специальных стендовых систем и уточнённой модели 
импеданса и пульсационной производительности по расходу. Показано, что исследуемый шестерённый насос 
внешнего зацепления рационально рассматривать по модели эквивалентного источника колебаний расхода. 
Проведена проверка предлагаемого подхода в виде расчёта амплитуд пульсаций давления в стендовых 
системах и сравнения с экспериментальными значениями. Отклонение расчётных амплитуд пульсаций 
давления от экспериментальных значений не превысило 10 %. 

Стоит отметить, что за рамками обзора осталось влияние на динамические характеристики шестерённого 
насоса приводной механической системы и износа деталей шестерённого насоса. На основе проведённой 
работы можно сформулировать направления развития предлагаемого подхода: 

— оценка динамических характеристик шестерённого насоса в течение всего срока службы; 

— разработка четырёхполюсной модели шестерённого насоса как эквивалентного источника колебаний, 
учитывающей приводную механическую и выходную гидравлическую системы; 

— разработка модели шестерённого насоса как эквивалентного источника колебаний, учитывающей его 
конструктивные особенности, процесс вытеснения, обратный гидроудар в межзубовую камеру, работу 
запертого объёма и перетечки из области нагнетания в область всасывания. 

При этом для подтверждения работоспособности предлагаемого подхода на весь класс машин объёмного 
вытеснения необходимо расширение номенклатуры испытуемых объёмных насосов. 
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Преобразование энергии гидравлического сопротивления системы 
в электричество 
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Аннотация 

Введение. Рассматривается гравитационный водовод с регулирующей задвижкой и регулятором давления в 
трубе самотечного водоснабжения в условиях, когда не нарушен гидравлический режим водопровода. 
Применительно к такой системе исследуются проблемы преобразования энергии местного искусственного 
гидравлического сопротивления в электричество. 

Материалы и методы. Изучена литература, в которой освещается возможность использования микротурбин 
для выработки электроэнергии в системах водоснабжения. Даны фактические значения, зафиксированные 
устройством непрерывной регистрации давления (логгер) в течение 12 часов, а также расход воды турбиной на 
данном участке узла (средний показатель в течение года) и перепад давления. Отмечено, что применение малых 
гидроэлектростанций в системах подачи воды значительно снижает стоимость их эксплуатации. Приведены 
показатели расхода воды в гидротурбинном агрегате в течение года с разбивкой по месяцам. Рассчитана 
максимальная мощность у входа турбины. Описаны принципы выбора типа гидротурбины. Указаны средние 
значения КПД для разных элементов гидротурбинной установки, усредненные параметры мощности МГЭС и 
соответствующие им показатели среднемесячной выработки электроэнергии. 

Результаты исследования. Оснащение исследуемых узлов специально разработанными турбинами позволит 
получить электрическую энергию, преобразуя искусственно созданную местным сопротивлением и 
погашенную механическую энергию. Подходы, описанные в данной статье, можно использовать при замене 
агрегатов регулирования давления системы водопроводной сети Еревана.  Спрогнозированы 
производительность МГЭС, сроки ее строительства и эксплуатации — 2 года и 30 лет соответственно. 
Предварительно рассчитаны затраты на сооружение и обслуживание. Приведены предполагаемые данные о 
доходах, расходах и чистой прибыли. Указано, что потенциально она будет близка к 6 млн драмов в год. 
Анализ данных позволил сделать вывод, что внутренняя норма доходности окажется на уровне 10,4 %, а срок 
окупаемости — 9 лет. 

Обсуждение и заключения. Рекомендуется заменить регулятор гидротурбиной с таким же гидравлическим 
сопротивлением и функцией автоматической регулировки расхода воды в системе. Выводы подтверждаются 
энергетическими и экономическими показателями гидротурбинной установки, расположенной на участке 
магистрального водовода «Арзакан — Ереван». 


Ключевые слова: гидравлическое сопротивление, энергия, гидротурбина, корректирующее устройство, 
водопровод, индикатор энергоэффективности, возобновляемые источники энергии, управление давлением, 
экспертные системы. 
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Введение. Традиционные невозобновляемые первичные источники энергии, ресурсы которых весьма 
ограничены, продолжают играть доминирующую роль в экономике Армении. Активное использование 
ископаемых источников энергии обусловило такие экологические проблемы, как загрязнение воздуха и почвы, 
дефицит воды и деградация экосистем. Следует учитывать рост численности населения, непрерывное бытовое 
и промышленное потребление энергии. Все это позволяет утверждать, что дальнейшее количественное 
(экстенсивное) развитие энергетической отрасли нецелесообразно. Использование возобновляемых источников 
энергии для пополнения общего энергетического баланса страны является императивом времени. Такой подход 
позволяет достичь устойчивого развития, сохраняя при этом окружающую среду. Водные и энергетические 
ресурсы необходимы для жизнедеятельности человека и подвержены регулярному экономическому, 
технологическому, демографическому и социальному воздействию [1]. По некоторым оценкам, 2—3 % мирового 
потребления электроэнергии приходится на системы водоснабжения (СВ) [2], и в 80-90% случаев 
используются насосные агрегаты. Расходы на их обслуживание — в числе основных операционных затрат 
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СВ [2]. Концепция рационального и эффективного использования водных ресурсов и электроэнергии имеет 
стратегическое значение для устойчивого развития и смягчения последствий климатических изменений. 
Рациональному использованию воды и электроэнергии препятствуют слабая инфраструктура и устаревшие 
эксплуатационные процедуры. Это особенно актуально для развивающихся стран. Более того, в соответствии с 
Целями развития тысячелетия' возникает потребность в альтернативных решениях: ставится задача сократить 
вдвое число людей, не имеющих постоянного доступа к чистой воде и основным услугам в области санитарии. 
В качестве возобновляемых источников энергии предлагается ввести в эксплуатацию малые 
гидроэлектростанции (МГЭС) на участках сети водоснабжения с местным гидравлическим сопротивлением при 
условии полуоткрытых задвижек или регуляторов давления. Это экономически выгодно для производства 
электроэнергии и водоснабжения и прекрасно иллюстрирует рациональное использование водных ресурсов. В 
прошлые десятилетия гидроэлектроэнергетика была одним из главных источников получения энергии из воды. 
Сегодня водное хозяйство рассматривается как прямой потребитель электроэнергии, и это влияет на 
распределение водных ресурсов”. Активное потребление энергии напрямую влияет на состояние водных 
ресурсов в мире, что в свою очередь является причиной климатических изменений. Рациональное 
использование энергии должно способствовать обеспечению устойчивого развития за счет грамотного 
распределения энергетических ресурсов на всех этапах преобразования. Для комплексной оценки проблемы 
приводится пример сооружения МГЭС на водопроводе «Арзакан — Ереван». 

Материалы и методы. Водопроводная ветка «Арзакан — Ереван» снабжает водой административные 
районы Еревана (Арабкир, Малатия-Себастия, Ачапняк, Эребуни) и прилегающие населенные пункты (Зовуни, 
Канакераван, Нор-Ачин, Нор- Гехи, Егвард и др.). Водопровод берет начало от источника Арзакан, а на 
территории административной общины Гетамедж разделяется на два водопровода с разными режимами 
давления. Рассматривается вопрос строительства МГЭС на одной из веток этого водопровода, питающего 
северо-западные и западные районы Еревана. Компания, работающая на разделительной ветке водопровода, 
установила полуоткрытую задвижку для регулировки расхода воды в системе. Так обеспечивается напор в 
водопроводе для населенных пунктов, расположенных выше основной магистрали (прилегают к городу Нор- 
Ачин). При полуоткрытой задвижке создается местное гидравлическое сопротивление около 13 м вод. ст. (до 
этого давление в водопроводной трубе — около 33 м вод. ст.). 

Предлагается вместо полуоткрытой задвижки установить МГЭС, не нарушая нормальный гидравлический 
режим работы водопровода. На рис. | представлен план местности в программе Соозе Еаг®, а в таблице 1 
сведены фактические значения, зафиксированные устройством непрерывной регистрации давления (логгер) в 
течение 12 часов. 


Рис. 1. План местности в Соо]е Еа 


' Цели развития тысячелетия / Организация объединенных наций : [сайт]. — ОВГ: В рз://у\\у.ип.ого/АеуеюртепЧеза/ли/тиИепиит- 
4еу@ортеп(-соа!5.В (дата обращения: 05.06.2022). 

Мао УизВЬ бВепе Нопе. Уапе Еааапе. Тфе ие со% оЁ соа. Опдегиапате СБштаз Епегоу Тапабсаре. 2008. ПВГ: 
Би рз://\ууи\ ипаег$апасЬтаепегоу.ого/отеепреасе-сШтафе-сЬа-зи$ашае-епегоу-ргосгату\и\Е-сЫпа-таш-аиогз-тао-уи$-зВеп?-Воп?- 
уапз-Рлайапз-2010-Шше-мие-соз-оЁ-соа1-ЕШ-герогИ (ассеззеа: 05.06.2022). 
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Таблица 1 
Данные регистраторов давления, установленных на исследуемом участке 
Зарегистрированное давление (м вод. ст.) 
Дата исследования Время — — 
До полуоткрытой задвижки После полуоткрытой задвижки 
15.09.2020 20:00 32,8 13,0 
15.09.2020 21:00 33,2 12,6 
15.09.2020 22:00 33,3 13,0 
15.09.2020 23:00 33,3 13,0 
16.09.2020 00:00 33,5 12,8 
16.09.2020 01:00 33,8 12,4 
16.09.2020 02:00 33,9 13,2 
16.09.2020 03:00 33,9 13,1 
16.09.2020 04:00 33,9 13,6 
16.09.2020 05:00 33,8 13,3 
16.09.2020 06:00 33,8 13,2 
16.09.2020 07:00 33,7 12,6 
16.09.2020 08:00 33,4 13,0 
Цель предложения —Й обеспечить преобразование искусственно созданного в локальном блоке 
сопротивления гидравлической энергии в электрическую. На месте строительства малой гидроэлектростанции 
(на ответвлении магистрального водопровода «Арзакан — Ереван») созданы благоприятные условия для 


комплексного использования водоэнергетического потенциала. Расход турбины на данном участке узла — 
1,15 м/с (средний показатель в течение года). Регулирующее устройство на трубе обеспечивает заданный 
расход. В результате образуется перепад давления от 20,3 м вод. ст. до 21,4 м вод. ст. (таблица 1). Следует 
отметить целесообразность практического использования так называемой «скрытой» энергии, образующейся на 
участках магистральных водопроводов в Республике Армения. Гидроэнергетический потенциал систем 
водоснабжения известен давно, однако не изучен должным образом. В литературе подробно описаны 
возможности использования микротурбин для выработки электроэнергии в системах водоснабжения [3—5]. 
Самотечные системы в районах с высокими уклонами рельефа формируют значительные давления в 
водопроводных и распределительных сетях, создавая предпосылки для выработки гидроэлектроэнергии. 
Турбины, установленные на водораспределительных сетях, можно использовать и как системы контроля 
давления. Для этого задействуют редукционные клапаны РКУ (от англ. ргеззиге гепеЁ уа]уе — клапан сброса 
давления), которые управляют потерями и утечками воды [6]. Они рассеивают энергию, таким образом снижая 
давление в системе. Что касается водяных турбин, они могут преобразовывать избыточное давление в 
электроэнергию”. Основные преимущества рекуперации гидравлической энергии в водопроводах, согласно 
Ф. Виейре и Х. М. Рамосу [3, 7], — повышение энергоэффективности системы за счет использования местных 
источников и снижение зависимости от внешней (сетевой) энергии. Кроме того, рекуперация гидравлической 
энергии способствует общему снижению эксплуатационных расходов. Ф. Виейра и Х. М. Рамос [3, 7] в своих 
исследованиях подчеркивают, что применение малых гидроэлектростанций в системах подачи воды 
значительно снижает стоимость их эксплуатации. Это объясняется тем, что предлагаемое решение позволяет 
заменить регулирующий клапан на равную по гидравлическому сопротивлению гидротурбину. В таблице 2 
приведены расходы воды в системе (указаны усредненные по месяцу показатели). 


Таблица 2 

Расход воды в гидротурбинном агрегате 

Месяц Расход, м/сек 
Январь 1,20 Июль 1,30 
Февраль 1,10 Август 1,30 
Март 1,10 Сентябрь 1,05 
Апрель 1,00 Октябрь 1,05 
Май 1,10 Ноябрь 1,10 
Июнь 1,30 Декабрь 1,20 


3 Тре МШеппит Оеуеортепе Соа!5 Вером 2011. Опйеа МаНопз. Ме\и УоК, 2011. ОВТ: 
Брз:/Ау\и\м ип. огелиШепиитеоа15/р9Е/(2011 Е)%20МРО%20Верог%202011 Воок%20Г.В РА (ассеззеа: 05.06.2022). 
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Учтем, что максимальный напор — Нштах = 21,4 м, минимальный — Нши = 20,3 м. 
р 3 

Среднегодовой расход, согласно табл. 2: О = 1,15 м’/с. 

При расчетном напоре, равном Н = 20 м, мощность потока у входа турбины [8]: 


М = раОН = 1000 - 9,81 : 1,15 20 =225632 Вт =225 кВт. 

Принципы выбора гидротурбины. При малом напоре и большем расходе (для МГЭС) следует выбрать 
турбину с большим коэффициентом быстроходности. Из реактивных турбин подойдет осевая. Из активных 
турбин малой мощности для пропуска относительно большого расхода следовало бы выбрать турбину Банки. 
Однако в данном случае недопустимо задействовать активную турбину, так как при ее наполнении у выхода 
формируется избыточное давление, и турбина Банки работает как очень плохая реактивная турбина с очень 
низкими энергетическими показателями. 

Рекомендуется установить гидротурбинный агрегат у задвижки на обходной линии. Входные и выходные 
задвижки обеспечат бесперебойную работу магистрали даже в случае отключения станции. Функционирование 
МГЭС будут контролировать датчики давления, установленные у входа и выхода турбины. Так обеспечивается 
требуемое давление на магистрали“. При установке оборудования на магистрали питьевого водоснабжения 
следует учитывать особые требования [8—10]. В частности, нужно исключить загрязнение воды смазочными 
маслами и иными материалами”. °. В целом же распределяемая вода должна соответствовать количественным и 
качественным нормам питьевого, бытового и промышленного потребления [11]. 

Энергетические показатели МГЭС. Для разных элементов гидротурбинной установки применены 
= 0,86 ‚ п.„ =0,94. 

Мощность на валу турбины: М, = п, М = 0,86 - 225 = 193,5 кВт, а на выходе генератора: Мен = Пхен №, = 0,94 
193,5 = 184 кВт. 

В таблице 3 представлены усредненные параметры мощности МГЭС и соответствующие им показатели 


следующие средние значения КПД: \ 


т 


среднемесячной выработки электроэнергии. 


Таблица 3 
Среднемесячная мощность и выработка МГЭС 

Месяц Мощность, кВт Выработка, кВт`час 
Январь 190 141360 
Февраль 174 116930 
Март 174 129456 
Апрель 158 113760 
Май 174 113456 
Июнь 206 148320 
Июль 206 153264 
Август 206 153264 
Сентябрь 166 119520 
Октябрь 166 123504 
Ноябрь 174 125280 
Декабрь 190 141360 

Итого 1573 474 


Следует отметить, что использование реактивной турбины МГЭС при большой отрицательной высоте 
отсасывания (й; = —1Зм), на взгляд авторов, в мировой практике предлагается впервые. В данном случае, с 


одной стороны, исключено возникновение кавитации [9, 12], что, безусловно, расценивается положительно. С 


* Маркарян А. Я., Токмаджян В. О. Регулирование производительности насосов с целью предотвращения кавитационных явлений // Вода: 
экология и технология : сб. докладов 7-го междунар. конгресса. М., 2006. Ч. 1. С. 566. 

°Запйагакап Капоппег уеу Моптег №2-Ш-А2-1, Кьтеш иг: гата{аКагаттап Кепгопа уаз Батакагоег! ! уогакш пегкауа’уоев ыетуемк 
рабаппег, Уогак1 Б5Кознигуип, Уегеуап, 20021, е] 11 [Санитарные правила и нормы №2-Ш-А2-1, Питьевая вода. Гигиенические требования к 
качеству воды централизованных систем водоснабжения, Контроль качества. Ереван, 2002. Стр. 11. На арм. яз.]. 

° Ваму, 1. А. \’аегеу: Епегеу апа У’аег ЕЁНс1епсу ш Митс:ра! Уаег бирру апа \У’аме\маег ТтеантепЕ Соз{-Е#есНуе бау!пе$ оЁ \аег апа 
Епегоу. ТВе АШапсе {о бауе Епегоу. \/азпе‘юп, 2007. 44 р. 
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другой стороны, неизвестно как влияет на КПД большое давление на выходе из турбины, поэтому в расчетах 
авторов данной статьи применен несколько заниженный показатель. 

Результаты исследования. Произведенную МГЭС электроэнергию (таблица 3) планируется поставлять в 
общую энергосистему Республики Армения на условиях гарантированной покупки по ценам, утвержденным 
Комиссией по регулированию общественных услуг Армении (в настоящее время действует тариф 
10,579 драмов за кВт‘час). 

Согласно данным таблицы 3, в течение года МГЭС произведет не менее 1,5 млн кВт‘час. По 
предварительным расчетам, срок строительства МГЭС — 24 месяца, а срок эксплуатации — 30 лет. Стоимость 


строительства и проектирования МГЭС составит 60,4 млн драмов (таблица 4). 


Таблица 4 
Затраты на строительство МГЭС 
Работы Стоимость, млн драмов РА 
Строительство / обслуживание трубопровода 2.2 
Строительство здания ГЭС 14,5 
Стоимость установки гидроагрегатов (гидротурбины, генераторы) 18,0 
Трансформаторная подстанция 10,7 
Строительство ЛЭП высокого напряжения 3,0 
Проектно-расчетные работы 7,0 
Иные затраты 5,0 
Итог 60,4 
Прогнозируемые доходы, эксплуатационные расходы и прибыль представлены в таблице 5. 
Таблица 5 
Прогноз годовых финансовых показателей 
Показатель Сумма, млн драмов РА 
Валовой доход 15,9 
Заработная плата, операционные расходы, текущий ремонт и ликвидация аварий 7,0 
Пошлины и прочие обязательные платежи 0,2 
Амортизационные траты (с учетом 30-летней эксплуатации) 2,0 
Иные затраты 0,2 
Прибыль до уплаты налога на прибыль 7,4 
Налог на прибыль 1,5 
Чистая прибыль 5,9 


Данные о годовых расходах, представленные в таблице 5, основаны на действующих нормативно-правовых 
актах Республики Армения. Согласно этим документам, налог на добавленную стоимость и налог на прибыль 
составляет 20 %7, а налог на имущество — 0,6 % (от остаточной стоимости). На основе данных о строительстве 
и текущих операционных расходах МГЭС в таблице б приведены основные показатели финансовой 
эффективности. 


Таблица 6 
Показатели финансовой эффективности 


Показатель Ед. изм. Значение 
Внутренняя норма доходности % 10,4 
Срок окупаемости лет 9,0 


7 Науауат! Напгарешуап Нагкаут огепзейк”, НО-165-М, упаипуа в 2016 РуаКат покетбен 4-т [Налоговый Кодекс Республики Армения, 
ЗА-165-Н. Принят 4 октября 2016 года. На арм. яз.]. 
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Обсуждение и заключения. Из-за гористого рельефа на большей территории Республики Армения для 
регулировки давления в водопроводной сети часто создают местное сопротивление. Как правило, речь идет о 
полуоткрытых клапанах и установке оборудования для регулировки давления. Оснащение узлов сети 
специально разработанными турбинами позволит преобразовать искусственно созданную местным 
сопротивлением и погашенную механическую энергию в электрическую. Следует отметить, что только в 
системе водопроводной сети Еревана установлено более 300 полуоткрытых задвижек и оборудования для 
регулировки давления. Многим из них требуется замена. На ремонтируемых участках можно применить 


решения, предложенные в данной статье. 
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Аннотация 

Введение. Движение носителя и внешние факторы (влияние атмосферы, температуры и давления) существенно 
ухудшают качество изображения смотрящих оптико-электронных систем (ОЭС) и точность позиционирования 
излучающих ОЭС. Рассматриваются проблемы повышения качества изображения и вероятность удержания 
изображения объекта наблюдения (ОН) на оптической оси следящей ОЭС. 

Материалы и методы. Создание системы автоматического управления оптико-электронного прибора (САУ 
ОЭП) предполагает решение многокритериальной задачи оптимизации с учетом ряда противоречивых технико- 
экономических (ТЭ) требований. Определение допустимых динамических погрешностей (ДДП) стабилизации 
изображения является ключевым вопросом в разработке бортовых оптико-электронных приборов (БОЭП). 
Уравнения Лагранжа П рода и смешанный метод Жильбера позволили получить математическую модель ОУ 
ОЭП. Затем выполнили декомпозицию двухсвязной САУ с нелинейными перекрестными связями в ОУ. 
Представлена функциональная схема модели формирования изображения БОЭП. Перечислены параметры 
матричного фотоприемника и требования, предъявляемые к динамической погрешности САУ ОЭП с учетом 
допустимой ФПМ ОЭП. Визуализированы функции передачи модуляции, линейного, гармонического и 
вибрационного сдвига изображения, соответствующие допустимым и достигнутым ДДП. В среде Ма@саа 
созданы логарифмические частотные характеристики. Двухсвязная система управления ОЭП при заданных 
параметрах ОУ для рассматриваемого движения представлена как два независимых канала управления по 
азимуту и углу места. 

Результаты исследования. Приведены процессы управления БОЭП в режимах стабилизации и слежения. Для 
исследования динамики пространственного управления БОЭП в соответствии с методикой САУ разработана 
компьютерная имитационная модель (КИМ) ЦСАУ ОЭП. Она реализована в среде МаЙаь и состоит из КИМ 
ОУ, приводов,  пропорционально-интегрально-дифференцирующих (ПИД) регуляторов с учетом 
нелинейностей, центрального вычислительного устройства (ЦВУ), программного устройства наведения, КИМ- 
носителя, реализующего уравнения движения. Описаны гармонические колебания носителя. Определены 
погрешности слежения и стабилизации в режиме слежения при дополнительно введенном управляющем 
воздействии в виде постоянной скорости. Исследована динамика пространственного управления БОЭП. 
Визуализированы компьютерная имитационная модель цифровой системы автоматического управления 
оптико-электронным прибором, результаты моделирования ЦСАУ ОЭП без учета движения борта и процессы 
управления ОЭП с учетом движения. 

Обсуждение и заключения. Для исследованных случаев рассчитана точность стабилизации. Установлено, что 
она в десятки раз превосходит заявленные ранее показатели, и это в десятки раз снижает требования к 
сходимости лазерного пучка и мощности излучения лазера при разработке оптического тракта 
рассматриваемого изделия. Предложенную КИМ можно использовать при разработке БОЭП. В этом случае 
применение представленной методики и КИМ поможет снизить трудозатраты и минимизировать ошибки. 
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Введение. При проектировании современных бортовых оптико-электронных приборов (ОЭП) и комплексов 
широко используются методы компьютерного моделирования оптико-электронных систем (ОЭС). Над этими 
задачами работали такие авторы, как Ю. Г. Якушенков, В. В. Тарасов, И. П. Торшина, В. П. Иванов, 
В. А. Балоев, В. А. Овсянников, В. Л. Филиппов. Компьютерное моделирование позволяет решать задачи 
рационального выбора структуры, параметров, элементной базы ОЭС, обеспечивающих требуемые показатели 
эффективности при заданных ограничениях без дорогостоящих натурных исследований и испытаний. 

Движение носителя и внешние факторы (температура, давление) существенно ухудшают качество 
изображения для смотрящих ОЭС и точность позиционирования для излучающих ОЭС. Сокращается 
вероятность удержания изображения объекта наблюдения (ОН) на оптической оси следящей ОЭС. Для 
излучающих ОЭС снижается точность позиционирования и вероятность выполнения задач наблюдения. Для 
удержания изображения ОН в поле зрения или в секторе облучения необходимо захватить довольно большой 
телесный угол обзора. Значит, нужно увеличивать габариты оптических систем и мощностей силовой 
электроники. Чтобы этого избежать, применяют управляемые излучающие ОЭС. Указанные выше особенности 
ОЭП ограничивают разработки САУ ОЭП. Эту проблему рассматривали В. А. Стрежнев, В. М. Матросов, 
А. С. Земляков, Н.Н. Маливанов, Е. И. Сомов, А. И. Маликов, В. А. Кренев, А. И. Карпов, Д. А. Молин, 
А. В. Михалицын. При проектировании управляемых ОЭП, работающих в режимах наведения и слежения, 
установлен следующий факт: время наведения и точность стабилизации оптической оси, а также динамика 
подсистем ОЭП существенно влияют на качество изображения. 

Таким образом, актуальна задача минимизации времени и увеличения точности наведения и удержания 
изображения ОН в поле зрения управляемого ОЭП. В этой связи необходимы: 

— разработка математических моделей ОЭП как объектов управления; 

— синтез алгоритмов управления; 

— создание компьютерных имитационных моделей (КИМ); 

— исследование систем управления, учитывающих динамику движения управляемого ОЭП; 

— определение параметров, влияющих на динамические свойства и качество изображения. 

Цель работы — улучшение точностных характеристик и повышение качества изображения ОЭП, 
работающих в режимах наведения, стабилизации и слежения, за счет рационального выбора их параметров при 
синтезе и моделировании. 

Материалы и методы. Разработка САУ ОЭП начинается с решения многокритериальной задачи 
оптимизации, учитывающей противоречивые технико-экономические (ТЭ) требования. В развитие работы [1] 
предложена модифицированная методика проектирования САУ ОЭП. За критерий качества САУ примем 
совокупность динамических характеристик каналов управления, удовлетворяющих условиям: 

Тоэп(Ун) > Тозп (Ун), (62) < 620" (и>), Док; < Дак", би, < Дбк”, (0 

М < МА", [Аф| > АА", [АЕ | > АГАоп. 

Здесь Т4°"(\„), Т.эп(ун) — допустимая и реализованная функция передачи модуляции (ФПМ) ОЭП на частоте 
Найквиста, 0(6); 6°" — амплитудная частотная характеристика перекрестных связей и ее допустимое значение 
декомпозиции; Да", Дак", Да, Да, — допустимые и установившиеся значения динамических погрешностей 
САУ по углу и угловой скорости при действии возмущений в условиях, близких к реальной эксплуатации САУ; 
к=1,2,... — номер канала управления, обеспечивающего качество изображения; Мхо = 1,05—1,25 — 
показатель колебательности, Аф; ДЕ, — запасы устойчивости по фазе и по амплитуде. 

Ключевой вопрос в разработке ОЭП — определение допустимых динамических погрешностей (ДДП) 


стабилизации изображения. 
Схема формирования изображения бортового оптико-электронного прибора (БОЭП) представлена на рис. 1. 
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Модель носителя 
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устройства 


МФП 


отображения 


Рис. 1. Функциональная схема модели формирования изображения БОЭП: ОН — объект наблюдения; ОС — оптическая 
система; МФП — матричный фотоприемник; УПУ — усилительно-преобразовательное устройство; ССИ — система 
стабилизации изображения 


ФПМ БОЭП должна а условию, обеспечивающему допустимое качество изображения [2]: 
Тоэс(у) = ТС) Тоб (Ти Туиу У)Тсси() > Тоэс ©, (2) 


2 (бес (м а ехр |- | 


16 т(2т) |" 

Здесь Тоэс(у) — ФПМ ОЭП; № — пространственная частота; Тоэс (у) — допустимая ФПМ ОЭП; Т.„(у) — ФИМ 
атмосферы; Тоб(у) — ФПМ объектива; Т»„(У) — ФИМ фотоприемника; Т„у(у) — ФИМ преобразования 
оптической информации [1]; Тсси() — ФПМ сдвига изображения (динамической погрешности ССИ), 
зависящая от вида динамического смещения изображения: линейного (Л) — х(Е) = УЕ, гармонического (Г) — 
х(Е) = а5т(®, случайного (СЛ) и расфокусировки. 

Допустимая ФИМ системы стабилизации изображения включает системы слежения, стабилизации, 
виброзащиты (СВ) и автоматической фокусировки (САФ) [3] 


р М 
Тсси(М) = Та(М)Те(М)Тв(М)Ть (М) > Теси(М) = м 
ТС») = тс аАалу), Те(у) = Ло (2пАаогу), Тв (У) = ехр[-2 (пАаву)? |, = 


2]. (А 21 (№). Меру еру 
Тьфу) = та Ду = 28п (“9”) (1- ти) 


1. (№) = (МТТ ФТ. 
Здесь Тл(М) — погрешность линейного смещения изображения; Т-(М№М) — амплитуда погрешности 
синусоидального колебания; Тв(№) — среднее значение погрешности амплитуды вибраций; ]› — функция 
Бесселя первого рода нулевого порядка; (4) — функция Бесселя первого рода первого порядка; Аз — 
погрешность фокусировки (мм); с — среднее значение волновой аберрации в долях длины волны; №», — 
средняя длина волны спектрального диапазона. 

Будем исходить из условия, что каждая подсистема САУ должна вносить одинаковую долю изменения 
качества изображения ОЭП и Т,.(М№) = сопзё за время наблюдения. Тогда упростится определение допустимых 


ФПМ ССл, СВ, САФ (3): 
вые = и (= Л,Г,В,Ф). (4) 


Разложим функции Тл(М№), Т-(М№), Тв(М№), Ть(М) (3) в ряд и получим выражения, определяющие ДДП: 


Е 0,824 тои Е 0,335 ТАЯТ) 
Аа А <= (5) 


Ун Ун 


[0.5 тн) 1 
де ЗОВ, = 28, Ар А 


ПУн 
Здесь у, — частота Найквиста. 
Предельная пространственная частота, которую должен разрешать ОЭП в процессе наблюдения, 
определятся критерием Джонсона: 


Упр = ^^, Ун = 0,5Упр, (6) 
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где упр — угловая (штр/рад) предельная пространственная частота; Мд — число элементов разрешения (числа 
Джонсона); [, — расстояние до ОН; й„, — критический размер ОН. 

При разложении отбрасываются все члены ряда, кроме первых, поэтому полученные допуски следует 
уточнить путем градиентного спуска. 

Представленные выше рассуждения позволили решить задачу определения точности стабилизации 
изображения бортовой оптико-электронной системы (БОЭС) для режимов слежения и стабилизации. Были 
разработаны методика и программа расчета ФПМ и допусков (патент В02021660340'). 

Назовем требования, предъявляемые к динамической погрешности САУ ОЭП с учетом допустимой ФИМ 
ОЭП. Исходные данные для расчета: 

— диаметр входного зрачка — О = 90 мм, 

— фокусное расстояние ОС — [`= 100 мм, 

— длина волны — 4. = 4 мкм, 

— расстояние до объекта наблюдения — [ = 10 км, 

— метеорологическая видимость — 5М = 15 км, 

— высота полета — Н = 5 км, 

— количество элементов МФП — 256х256, 

— частота кадров — Ё = 400 Гц, 

— вероятность правильного обнаружения — Рон = 0,8, 

— вероятность ложной тревоги — Ри. = 10“, 

— отношение сигнал/шум усилителя — оптимального фильтра — т = 1,105, 

— угловой размер объекта наблюдения — у; = 5:10 ° рад, 

— угловое расстояние между двумя соседними целями — у, = 3,34: 10 3 рад. 

Параметры МФП: 

—4=0,25а — толщина фоточувствительного слоя, 

— =, = 4-10? — неэффективность переноса заряда, 

— 0, = 2 — число выборок на 1 элемент, 

— К, = 0,6 — коэффициент заполнения, 

—т= 5,4: О е— время считывания, 

—Тр=5: 10 '%с — постоянная времени преобразователя. 

На рис.2 приведены ФИМ линейного, гармонического и вибрационного сдвига изображения, 
соответствующие допустимым и достигнутым ДДП САУ БОЭП. 


г” 
0,8 
0,6 ) Тлин. (у) 
Теары. (у) 
0,4. — Тибр. () 
0,2 
0. 
0 100 200 300 400 500 600 у, Ирад 


Рис. 2. Функции передачи модуляции Ту (у) 


В соответствии с (3)—(5) и исходными данными для расчета вычислили: У„р = 600 рад ' канала наблюдения 
БОЭП в инфракрасной области с вероятностью обнаружения Р=0,8. Для ФИМ БОЭП и ДДП САУ: 
дал = 3,5 угл. мин, До: = 3,1 угл. мин, Дов = 3,0 угл. мин. 


Бир://уезиК-доп$а.ги 


' Программа вычисления функции передачи модуляции оптико-электронного прибора, базирующегося на БПЛА : патент 2021660340 Рос. 
Федерация / К. А. Бурдинов ; Университет Иннополис. № 2021619174; заявл. 11.06.2021 ; опубл. 24.06.2021, Бюл. №7, 1 с. 
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На основе уравнений Лагранжа П рода и смешанного метода Жильбера получили математическую модель 
ОУ ОЭП [4] наведения и слежения в виде трех информационных каналов, приводимых в движение 
электродвигателями типа ДБМ по осям Карданова подвеса. Представили ее в матричном виде”: 

А(4)4 + В(&, 4)9 + Е(а, 4) + И’(&, а) + 9(&,4) = Му — Мур. (7) 
Здесь а = (а ВТ, а, В — углы поворота ОУ по азимуту и углу места. 


_ (В +В(В) —Р(В) _ 0 Вов (Е, а, В) 

А = ( —р(В) С. 869 = (ьзбав 0 ) 
‚, _ (-2аВР (а, В) — В2Е(В) _ (№6 а) 
—— @Р (В) } ис = (ша: 


г„(а) ра — — Млдвл № Мутр.151(6) 
04-4 ичААЕ). М» = (м. —- мт) 
Вов (Е, в В) = {(2А>(В) — С») соз а — 2Е (В) зт а}, (Е) — 2Е (В), (Е) 
— {(2А>(В) — С>) та + 2Е(В) с0$ 4, (6), 


ии (6,4) = [© (ОЕ@, В) + ву, (О (В, + В(В)) — „(© (о, В} + [Е(а, В). (6) — Е(а, В)» (6) + 
+2, Во(оу, (Е) + 241 (а, Вдох, (Во,, (©) — (о, В)», (в, (6), 
из (6, 9) = [@», (ОЕ (®) соз(а) — С, зт(е0] — Фу„(ВР(В) + 6, (Е (В) 5й (о) + С» с05(60]} — 


1 
- 5 (а + с05(20)) + р (В) зт(2) ох, (Е) — 2Е(В)оу, (Е) + [Е(В)(1 — с05(20)) — Р(В) зт(2а)]в52, (Е) + 
+2[2А› (В) со5() — Е(В) зо)», (Е) у, (Е) — 2[Е(В) зт(2о) — Б(В) со5(2а) (Е), (Е) — 
—2[Е (В) со$(а) + 2А>(В) $т(а) в. (Е), (6)} 


—хс, 5т(а) + 2с, с05(@) ^. 
т, (9) = ту 0 + 
—хс, с05(@) — 2с, 5 (а) 


—(хе, со5(В) — Ус, 5т(В)) 5т(а) + 2с, с05(@) . 
+т> 0 р 
— (хе, со5(В) — Ус, 5т(В)) с05(а) — 2с, 5т(а) 


(хе, тт (В) — ус, со(В)) соз(®0\" 
тв (а) = т> хс, со5(В) — Ус, (В) 
(хе, 5т(В) + Ус, со5(В)) 5т(а) 


Следующий этап разработки — декомпозиция двухсвязной САУ [5] с нелинейными перекрестными связями 
в ОУ, каким и является исследуемый ОЭП. Линеаризованные уравнения движения относительно траектории 
имеют вид: 

а хХьх (Е) = бЕЕО - 2п; В х„, (Е), = ВЕ ЕО + п/2, 
(ало + @11)А& + (Бо + В: 1) + с. Аа + аз АВ + В. АВ + с12 АВ = КА, 
ал Да + Ь.1Аб + ао ДВ + БоДВ + с›2АВ = К.2Ди», (8) 
Ди; = А, (р)Да!, Да, = Дав, — Да’; Ди» = В›(р)ДЕ», Де» = ДВь, — ДВ". 

Построив переменные коэффициенты системы (8), определим их значения при { = #. в «опасных» точках, в 
которых параметры значительно изменяются или меняют знаки в диапазоне углов (В’=0 + л/2, а’ =0-2л). 
Тогда, сведя систему (8) (при #{=&) к системе с постоянными параметрами и используя правило Крамера, 
выразим ее в форме передаточных функций (ПФ) в точках #= 1, : 

Да = И’ (р, «Ди, — И’, >(ф, «Ди», (9) 
ДВ = И/>2(ф, &)Ди» — И’, 1 (ф, ди. р 

Полученные ПФ (9) с учетом (8) для каждого выбранного момента времени (1= 1.) можно представить в 
виде структурной схемы с прямыми перекрестными связями (ПС, рис. 3). Здесь и далее для простоты записи 
ПФ обозначим Ир) = Укр, ых 


? Математическая модель и синтез системы автоматического управления бортового оптико-электронного прибора / К. А. Бурдинов [и др.] // 
ХПИ Всерос. съезд по фундаментальным проблемам теоретической и прикладной механики : сб. тр. Уфа, 2019. С. 190-192. 
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Рис. 3. Структурная схема двухсвязной САУ 


Оценим ПС САУ по их частотным характеристикам (ЧХ). Для этого запишем ПФ разомкнутой системы с 
учетом 2-го замкнутого канала управления при ДВь, = 0 (рис. 4). 


1 >) 


Рис. 4. Структурная схема САУ, разомкнутая по одному каналу управления 
Построим годограф И (76) = В, (°)И (6) и трубки вокруг него с радиусами для каждого момента 
времени #. : 
&1к (6) = А! (©) Ухи (6) 5х (6) 6, (6). (10) 
Учитывая (10) и критерий Найквиста, можно судить об устойчивости САУ (1), (2) и о влиянии ПС на ее 
устойчивость в диапазоне углов сканирования (3). 
В результате для каждого (1, -го) определим время ПС: 


Иа2 (6) М/ка1 (©) лир 
бк: (©) = 6,12 (©) 6.21 (6) ИЙ.; (©) = и о) И), С] = 1,2). (11) 


Рассмотрим ПС САУ, удовлетворяющие условиям декомпозиции. 

Согласно (8), (11) по данным геометрии масс ОУ, для 11 выбранных точек по времени определили (1. = 0; 
0,89; 1,35; 1,58; 1,71; 2,3; 2,97; 3,8; 4,5; 4,75; 6,28 с). Построили а;(#), В(0, с5®. В среде Маса4 создали 
логарифмические частотные характеристики (ЛАХ) 6,(6) (К = 1,11) (рис. 5). 


Амплитуда, дБ 


102 10° 10? 10° Частота, рад/с 


Рис. 5. Логарифмические амплитудные характеристики 6,(0) 


На рис. 5 видно, что все множество ЛАХ не превышает —50 дБ. Из анализа ЛАХ 6,(&) имеем 
представленные ниже результаты. 


УоеО, = (10`*=10°) рад/с. 
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Определены: тах бк(6) = 0,0136 (-37,3 дБ) < 6*°" = 0,0335 (-29,5 дБ), АЁ = М" / (Ми + 1) = 0,535 (-5,43 дБ) 
при Мю= М»= 1,1; М" = 1,15. 

Таким образом, двухсвязную систему управления (8) ОЭП при заданных параметрах ОУ для 
рассматриваемого движения можно представить как два независимых канала управления по азимуту и углу 
места. 

Результаты исследования. На рис. 6-7 показаны процессы управления БОЭП в режимах стабилизации и 
слежения. Масштаб графиков рассогласования увеличен в 600 раз. Графики программного управления и 


выходной величины визуально совпадают. 


30 


2Даг°" 
г 
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Рис. 6. Процессы стабилизации БОЭП по азимуту 
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Рис. 7. Процессы слежения БОЭП по азимуту 


Опишем гармонические колебания носителя: А, = А›=2° (амплитуды), Т! = Т»=1с (периоды колебаний) и 
А, = А, = 12%, Т, = Т, = 6,28 с. Укажем моменты нагрузки: Мхь! = 2,2 Нм, Мы = 0,3 Нм. При таких гармонических 
колебаниях и моментах нагрузки погрешность стабилизации оси визирования не превышает: 

— по азимуту — 1,2 угл. мин, 

— по углу места — 1,6 угл. мин. 

В режиме слежения при дополнительно введенном управляющем воздействии в виде постоянной 
скорости 12'7^/с погрешность слежения и стабилизации не превышает 1,7 угл. мин по азимуту и 1,3 
угл. мин по углу места. 

Для исследования динамики пространственного управления БОЭП в соответствии с методикой САУ разработана 
КИМ цифровой системы автоматического управления (ЦСАУ) ОЭП конкретного назначения (рис. 8) [1]. 
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Рис 8. Компьютерная имитационная модель цифровой системы автоматического управления оптико-электронным прибором 


Решение реализовано в среде МаНаЬ ЗипаПикК. Его элементы: 

— КИМ ОУ (ОБес®, полученная согласно (7); 

— КИМ приводов (Бпуе); 

— КИМ ПИД-регуляторов (Веза‘ог) с учетом нелинейностей и ЦВУ; 

— КИМ программного устройства наведения; 

— КИМ носителя, реализующего уравнения движения носителя. 

На рис. 6 и 7 приведены результаты моделирования ЦСАУ ОЭП без учета движения борта. При этом: 

— погрешность наведения не превышает по азимуту 18 угл. мин, по углу места — 10 угл. мин; 

— погрешность слежения — 0,6 угл. мин по азимуту, 0,6 угл. мин по углу места; 

— погрешность стабилизации — 0,45 угл. мин. по азимуту, 0,6 угл. мин по углу места. 

На рис. 9 показаны процессы управления ОЭП с учетом его движения. Погрешность наведения по 
азимуту — 2 угл. мин, по углу места — 5 угл. мин. Погрешность слежения и стабилизации — 0,8 угл. мин по 
азимуту, 1,3 угл. мин по углу места. 
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Рис. 9. Результаты моделирования системы в режиме наведения по двум каналам (градусы): а — рассогласование системы; 
158 б — графики переходных процессов. Во втором случае визуально совпадают график программного управления 
(арВа Еагое, Бей {агре!) и выходной величины (а!рБа, Беа) 


Бурдинов К. А. и др. Исследование системы автоматического управления стабилизации изображения 


Обсуждение и заключения. Полученные результаты представлены ниже. 

1. Предложена методика оценки допусков на точность стабилизации изображения БОЭП в условиях 
управляемых и температурных воздействий (патент К 92021660340). Представлены аналитические оценки ДДП 
стабилизации изображения по (4), (5) для режимов слежения и стабилизации на основе ФИМ ОЭП. Получены 
требования к САУ, обеспечивающие решение задачи обнаружения ОН в инфракрасной области. 

2. Назван критерий приемлемости декомпозиции, описана методика декомпозиции двухсвязной САУ ОЭП с 
применением компьютерных технологий. Для конструктивных параметров ОУ реального прибора определены 
требования, позволяющие синтезировать САУ ОЭГ как два независимых канала управления. 

3. Предложены КИМ САУ ОЭП для исследования динамики изолированных каналов управления и 
пространственного управления. При этом учтены нелинейность и нестационарность ОУ конкретного назначения. 

4. Разработана механическая модель ОЭП, получены уравнения динамики пространственного движения (7) ОУ. 

5. Синтезированы параметры ПИД-регуляторов, обеспечивающие устойчивое наведение и стабилизацию 
ОЭП при всех заданных режимах. Добротность по скорости: К, = (3700 = 4000)с`*, К,› = (3000 = 
3300) с *. Корректирующие звенья: Тнл1 = Тнл2 = 0,05с, Туо1 = Т.22 = 5 Х 10$ с. Это обеспечивает более 
точную работу системы планирования полета и построение траектории [6-8]. 

6. Показано, что с помощью КИМ можно получить требуемые время и погрешность наведения, 
стабилизации и слежения. При этом скорость наведения — 0,6 с. Соответствующий показатель зарубежных 
аналогов — 1 с*. 

7. В [9] рассматривается система лазерного противодействия. Ее технические характеристики позволяют 
выполнять задачу при У = 1000 Вт/стерадиан и мощности лазера Р = 3,1 Вт. В [19] показано, что необходима 
мощность излучения И = 2700 Вт/стерадиан при мощности лазера Р = 200 мВт. 


Рассчитаем требуемую ТОЧНОСТЬ стабилизации для каждого из представленных случаев по формуле А = 


2х180х%60 Р 
—— агссо$ (1 в —> =) [угл.мин] и сведем в табл. 1. 
п 


Таблица 1 
Сравнение требуемой точности стабилизации 


Мощность излучения | Мощность лазера | Требуемая точность стабилизации 
(И/), Вт/стерадиан (Р), Вт (д), угл. мин 
Установка [11] 2700 0,2 33,4 
Установка [10] 1000 Е 216,1 
Изделие 1,3 


Таким образом, ТОЧНОСТЬ стабилизации предлагаемой системы в десятки раз превосходит требуемую В 
статьях [9, 10]. И это в десятки раз снижает требования к сходимости лазерного пучка и мощности излучения 
лазера при разработке оптического тракта рассматриваемого изделия. 


Список сокращений и условных обозначений 
ОЭС — оптико-электронная система, 

ОЭП — оптико-электронный прибор, 

ОН — объект наблюдения, 

КИМ — компьютерная имитационная модель, 
САУ — система автоматического управления, 
ФПМ — функция передачи модуляции, 

ДДП — допустимая динамическая погрешность, 
ССИ — система стабилизации изображения, 
МФП — матричный фотоприемник. 


3 риесеа шёгагей Соищениеазиге (РВСМ) / ЕМ$ОРЕПГА // етзореФа.ого : [сайт]. — ЕВГ: Врз://м\у м. етзоре41а.ого/епез/Атес(еа- 
шйагед-соииегтеазиге-Чист/ (дата обращения: 01.05.2022). 
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Аннотация 
Введение. Зарождающиеся усталостные повреждения в металлических пролетных строениях мостов создают 


определенные угрозы безопасности эксплуатации. Для их своевременного выявления и диагностики 
используются различные методы неразрушающего контроля. Современным и популярным решением на 
данный момент является метод инфракрасной (ИК) термографии. В силу особенности работы ИК-камер для 
получения точного результата требуется дополнительная обработка записей, полученных с этих камер. Цель 
данной работы — представить метод обработки термофильмов и показать возможности его применения в 
реальных условиях. 

Материалы и методы. Приводится описание способа обработки термографических фильмов, позволяющего 
выявить температурные аномалии, используя только информацию с камеры. Приведены результаты его 
применения на элементах действующих металлических пролетных строений мостов. 

Результаты исследования. Показано, что для существующих дефектов имеются температурные аномалии. 
Это означает, что дефекты продолжают развиваться, что подтвердилось последующими наблюдениями за их 
состоянием. Кроме этого был выявлен случай температурной аномалии в бездефектной внешней области. Это 
может быть признаком зарождающегося дефекта, который не мог быть продиагностирован другими методами. 
Если наличие этого дефекта будет подтверждено при повторных обследованиях, то представится возможность 
диагностирования скрытых дефектов, которые еще не вышли на поверхность, и/или обнаружения потенциально 
разрушающихся мест. 

Обсуждение и заключения. Показана эффективность ИК-термографии как метода бесконтактного неразрушающего 


контроля, а также его работоспособность на реальных объектах, находящихся под случайной нагрузкой. 


Ключевые слова: ИК-термография, неразрушающий контроль, усталостные трещины, металлические мосты, 


дефектность конструкций, ИК-камера. 
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Введение. Усталостные трещины в металлических конструкциях мостов, образующиеся под воздействием 
нагрузки от транспорта, ветровой нагрузки и т. п., представляют серьёзную проблему. Хотя эти трещины сами 
по себе влияют только на долговечность сооружения, при определенных условиях они могут дать толчок 
развитию хрупких трещин, что уже непосредственно влияет на безопасность использования конструкций [1]. 

Одним из современных методов неразрушающего контроля, набирающего популярность, является ИК- 
термография [2-12]. Принцип метода основывается на том факте, что любое тело, испытывающее механическое 
воздействие, выделяет количество тепла, пропорциональное интенсивности воздействия. Поэтому, при 
обследовании участка конструкции, содержащего дефекты, можно сразу же получить информацию о 
распределении напряжений по этому участку. Отличающееся температурное распределение будет являться 
индикатором протекающего усталостного повреждения [13, 14]. 

Основной сложностью при работе с ИК-камерами (тепловизорами) является получение результатов с 
приемлемой точностью. В силу особенностей строения детекторов наибольшая точность достигается при 
использовании камер с охлаждаемыми матрицами. Однако слишком высокая стоимость подобных 
аппаратов является серьёзным препятствием для их использования. Появление камер, использующих 
неохлаждаемые массивы микроболометров, улучшило ситуацию, поскольку эти камеры оказались заметно 


дешевле. Недостатком таких тепловизоров стала меньшая точность по сравнению с камерами с 
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охлаждаемыми матрицами. Поэтому для получения устойчивых результатов требуется дополнительная 
обработка получаемых термограмм. 

В данной работе приводится алгоритм корреляционной ИК-термографии с применением сигнала, 
полученного непосредственно с тепловизионной съемки', а также некоторые результаты его применения. 


Материалы и методы. Метод ИК-термографии основан на анализе распределения температуры в 
материале. При воздействии растягивающих или сжимающих напряжений на материал его температура будет 
изменяться пропорционально величине этих напряжений в соответствии со следующей формулой:? 


[02 
АТ =-—Т.Ад, (1) 
Р.С, 
где а — коэффициент теплового расширения; р —ЩЩ плотность вещества; С, Е: удельная теплоемкость при 


постоянном давлении; Г — абсолютная температура, К; Ас — изменение главных напряжений, МПа. 

При сохранении адиабатических условий, а также при развитии только упругих деформаций связь (1) 
является линейной и обратимой, при которой изменение температуры следует за изменением напряжения. 
Как только в материале начинают зарождаться повреждения, он перестает быть упругим и развивает 
пластические деформации. Энергия системы, подверженной таким деформациям, интенсивно переходит в 
тепловую, привнося таким образом основной вклад в тепловой эффект. Таким образом, тепловая картина 
становится отражением внутреннего состояния материала (степени разрушения) и, следовательно, 
показателем уровня напряжений. 

Определить относительную амплитуду ИК-сигнала можно с помощью метода наименьших квадратов”. 
Положим, что на конструкцию действует нагрузка Е, воздействие которой в некоторой области можем 
измерить и обозначим как Кд. Будем считать этот сигнал опорным. Пусть также имеется ИК-видеосъемка этой 
конструкции в тот же промежуток времени, которую можно представить как последовательность кадров 
(кадр — матрица значений размером М хК) во времени и обозначим у,(1) значения ([,]) элемента этой 
матрицы (пикселя) в момент времени #, где 1<71<М ‚, 1<}<К. 


Обозначим приближающую функцию для этого элемента как: 
У (()= а, +Ь, 1 (1). (2) 


где а; — смещение; Ь, — коэффициент влияния опорного сигнала {. 

Будем приближать с помощью функции (2) сигналы у;(1) из тепловизионного фильма. Покажем, как найти 
эти коэффициенты для некоторого у,(1, для остальных —Щ рассуждения будут аналогичны. Будем 
минимизировать отклонение нашего приближения от сигнала, полученного с камеры: 


А = >6,()-(® — ши, (3) 


п= 


где М — количество кадров в записи; п — номер кадра, соответствующий времени. Отсюда Б,; может быть 
найдено следующим образом: 


№ 


>) Хз) музее) -У лети) 


Ь. —__ п п=1 = п=1 п п=1 (4) 


баКагапи Т, МазВипига Т, КиБо $, её а!. Оеуеюртепе оЁ а ЗеМ-геЁегепсе Госк-ш ТвегтостарВу юг Ветое Мопезииснуе Тезние оЁ Еанеие 
СгасКк (15 Веро, Рипдатеша! З4у Озше \Уе4еа Зее! батр!ез). Тгапзасйопз оЁ фе Ларап бослебу оЁ МесБашса! Епошеег$. Зейез А. 
2006;72:1860-1867. В рз://401.ог2/10.1299/Кала.72.1860 

^ТГВотзоп У". (Гога Кеш). Оп ше дупапуса! феогу оЁ Неае. Тгапз. В. бос. ЕашпьагеН. 1853;20:261-283. 

3 Тезшак ЛВ, Воусе ВВ, Но\епугаег С. ТвегтоеЙазис Меазигетепе Опаег Вапаот Г.оа@ ше. ш: Ргос. ЗЕМ 5риие СопЁ, 1998. Р. 504-507. 
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Повторяя (3) и (4) для всех (1, /) получим матрицу В = |} того же размера, что и исходный кадр съемки. 


Таким образом получим приближение каждого сигнала у;(В) ИК-съемки по времени посредством опорного 
сигнала КЕ). Значения матрицы В показывают относительную интенсивность изменения температуры в 
некоторой области по сравнению с интенсивностью изменения опорного сигнала. Опорный сигнал КА) не 
обязательно может быть получен со сторонней системы, например, с тензодатчика. Описанный подход, 
именуемый автокорреляционной ИК-термографией, может быть применен и с опорным сигналом, полученным 
с той же самой ИК-записи“ [15, 16]. 

Результаты исследования. Реализация предложенного алгоритма была выполнена в виде Руфоп скриптов. 
С помощью этих скриптов производилась обработка термографических фильмов, записанных на металлических 
пролетных строениях действующих автомобильных мостов, в конструкциях которых были ранее 
диагностированы усталостные трещины. Проверка работоспособности алгоритма была проведена на стендовых 
испытаниях, имитирующих поведение реальных конструкций с известной информацией о дефекте. Съемка 
проводилась с помощью инфракрасной камеры с неохлаждаемым микроболометром ЕшКе Т!1400, имеющей 
тепловую чувствительность менее 0,05 °С и частоту съёмки 9 Гц. Записи сделаны в момент воздействия 
автомобильной нагрузки на мост. В качестве опорных сигналов для каждой записи использовалась 
бездефектная зона размером 15х15 пикселей с того же конструктивного элемента недалеко от существующей 
и/или предполагаемой зоны разрушения. Результаты обработки наглядно представлены в виде изображений, 
полученных посредством значений матрицы В. Каждому значению матрицы соответствует один пиксель на 
изображении, преобразование числа в цвет осуществляется с помощью функции палитры из определенных в 
пакете таро т. 

На рис. | изображен фрагмент стальной балки пролетного строения. В месте соединения горизонтального 
листа со стенкой обнаружена трещина (рис. 1 в). На рис. 1 а показан один из кадров ИК-фильма этой трещины. 
На рис.16 представлен результат обработки, на котором отчетливо видно пятно, свидетельствующее о 
разогреве, происходящем в обозначенной зоне. Эта зона соответствует области вокруг вершины трещины, что 
свидетельствует о продолжающемся развитии дефекта. Это же было подтверждено повторными 
обследованиями, которые установили увеличение трещины на более чем 30 % за 4 года. 


Рис. 1. Результаты инфракрасной съёмки: 
а — кадр фильма; б — изображение, построенное по значениям матрицы В; в — фотография трещины 


4 СаНеш 0, Модиепо О, браепо!о Г.. А поуе! па! ргосеззте тео4 юг ТЗА аррНсанопз Меазигетеп! $с1епсе ап4 Тесвпооу. 2005;16:2251. 
10.1088/0957-0233/16/11/017 
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На рис. 2, 3 изображены трещины, зародившиеся в сварном шве, затем раздваивающиеся и уходящие в ребро 
жесткости и лист настила. Съемка проводилась во время испытаний моста. Для контроля напряжений в этих же 


местах были установлены тензометрические датчики, результаты с которых подтвердили развитие этих трещин. 


а) 6) 
Рис. 2. Диагностированная усталостная трещина: 
а — изображение, построенное по значениям матрицы В; 6 — фотография 


Контрольные метки для 
совмещения изображений 


в) 
Рис. 3. Диагностированная усталостная трещина: 
а — кадр фильма; б — изображение, построенное по значениям матрицы В; в — фотография трещины 


На рис. 4 изображены совмещенные показания тензодатчиков и температуры для одного и того же участка 
конструкции. Температурные показатели были получены усреднением в области обследования и последующим 
сглаживанием по времени. Тензометрические данные также были сглажены. Как видно, что изменение 


температуры полностью повторяет изменение напряжений с точностью до чувствительности прибора. 
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Нормированные значения, ед. 
> 
л 


Изменение напряжений 
Изменение температуры | 


| | | 
5 10 15 20 Время, с 


Рис. 4. Графики изменения напряжений, совмещенные с изменениями температуры (изменения температуры приняты с 
обратным знаком). Величины по оси ординат нормированы относительно единицы 


На рис. 5 изображен элемент конструкции, не имеющий ВИДИМЫХ повреждений, которые также не были 
однозначно определены другими методами. Однако, в результате обработки пакета термограмм было 
обнаружено, что в месте соединения двух элементов имеется тепловая аномалия. Это свидетельствует о том, 
что в этом месте вероятно идет саморазогрев материала, и под поверхностью происходит зарождение дефекта, 
который в скором времени должен будет выйти на поверхность. Таким образом, данный участок конструкции 
требует дополнительного внимания при будущих обследованиях. 


`, 


а) 6) 
Рис. 5. Результаты инфракрасной съёмки: 
а — сварной шов без видимых усталостных повреждений; 
б — изображение носка сварного шва, построенное по значениям матрицы В 


Обсуждение и заключения. Проведенные работы на пролетных строениях показали эффективность 
автокорреляционной ИК-термографии как метода бесконтактной диагностики. Показана работоспособность 
метода на реальных объектах, находящихся под случайной нагрузкой. Имело место соответствие изменения 
напряжений в конструкции с изменением регистрируемой температуры. Проведена диагностика некоторых 
известных дефектов, на которых было отмечено продолжающееся развитие повреждений, что в последствии 
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подтвердилось повторными обследованиями этих элементов конструкций. Кроме этого, выявлен случай 
саморазогрева во внешне бездефектной зоне, что указывает на вероятный процесс зарождения дефекта, 
который на момент съемки не мог быть однозначно диагностирован другими методами. Это свидетельствует о 
том, что данный участок требует дополнительного наблюдения. При условии обнаружения дефекта на этом 
участке в процессе последующих наблюдений можно говорить о возможности использования ИК-термографии 
для выявления скрытых дефектов и прогнозирования их появления. 
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Аннотация 
Введение. В статье представлен обзор современных моделей нейронных сетей для обработки естественного 169 
языка. Исследования обработки естественного языка представляют интерес в связи с тем, что в последнее 
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время возросла потребность в обработке больших объемов аудио- и текстовой информации, накопленной 
за последние десятилетия. Наиболее обсуждаемой в зарубежной литературе является особенности 
обработки разговорной речи. Цель работы представить современные модели нейронных сетей в области 
обработки устной речи. 

Материалы и методы. Прикладное исследование понимания устной речи является сложной и далеко идущей 
темой обработки естественного языка. Понимание на слух занимает центральное место в исследовании и 
представляет собой проблему. В этой статье предлагается метод понимания на слух, основанный на глубоком 
обучении. В статье кратко излагаются содержательные аспекты различных этапов создания нейронной сети по 
распознаванию речи, приводятся основные термины, связанные с этой теорией. Приводится краткая 
характеристика основных существующих на текущий момент нейронных сетей по обработке естественного 
языка. 

Результаты исследования. Проведен ретроспективный анализ зарубежных и отечественных литературных 
источников с описанием новых методов обработки устной речи, в которых использовались нейронные сети. 
Предоставлена информация о нейронных сетях, методах распознавания и синтеза речи. В работу включены 
результаты разноплановых экспериментальных работ последних лет. В статье подробно описаны основные 
подходы к обработке естественного языка и их изменения с течением времени, а также появление новых 
технологий. Рассмотрены основные проблемы, существующие в настоящее время в этой области. 

Обсуждение и заключения. Анализ основных аспектов систем распознавания речи показал, что в настоящее 
время не существует универсальной системы, которая была бы самообучающейся, шумоизоляционной, 
распознающей непрерывную речь, способной работать с большими словарями и в то же время имеющей низкий 
уровень ошибок. 


Ключевые слова: обработка естественного языка, устная речь, нейронные сети, автоматическая обработка 
естественного языка, семантическая согласованность. 
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Фае аге р!аппеа тю Ъе апа[у7е4, ап4 а согриз оЁ {ех{5 1$ сгеме4. Мех( %ер, Фе соЦеце таена| 15 гапзРегге4 ю ап 
ехрегЕ Ппхи15. Не ргезспЬез Ше гШез, сотрПез @сНопанез$, татКз$ ир {ех5 Юг Ше 14епйЯсайоп оЁ (агое( згтасаге$ т 
{(ех(5 Юг Ше Раег ‘1азК зошНоп. Апо@ег тео 1$ а15о изе4, 1 \сВ ап ехрегё Пп513ё плагК$ Ше {ехё шю тагое 
зегисгез ог саезоп7е$ {ех{5 ш сецаш с1а5зез, ап еп пасте 1еагито тео4$ ащотансаПу депуе зоте гШез ог 
то4е]з ог Виег зо сиггепе ргоетз. АЕ Ше еп4 оЁ Ше \уотК, Ше апаШу оЁ Фе те@о4$ 1$ а!\\ауз сБеске4. 

РЬ1о105155 зу зетайис$ оЁ Ше {ехЕ сопу14ете теап1пз$ оЁ ро!узетаййс ип$ ш сощехё етрВаз171ие Фа 
сощехЕ р1ауз а Рапдатепа! го]е 11 Ше \уога дейшиоп. ТБегеоге, е.5., Ше аиог5 Чзсоуег Фе сощехаа! пеап1п25 оЁ 
роузетапйс ипИз Шаё аге по! геол еге4 ш 1ех1соотарЫс зоигсез. ш Ше еайу $азез, зс1епН5{5 ргорозе4 тю 4туае апу 
зещепсе ш® а $е{ оЁ \’ог4з Ша сошШ4 Ъе ргосеззе4 шату14иаПу, \Ь1сВ \аз шисЬ еаз1ег ап ргосеззте а вое 
зещепсе. ТЬ$ арргоасН 15 зипПаг №0 Фе опе изе4 © {еасВ а пе\м 1апсиасе 1ю сБПагеп ап4 ада$. \Увеп уме НгзЕ 5 
1еагио а 1апоцаое, \е аге шио4исе4 © Из рак оЁ зреесВ. Ге! и5 сопз1Аег Еп?ИЗВ аз ап ехатр/е. [ Ваз 9 таш рат 
ОЕ зреесВ: пои, уетЬ, а@есйуе, а4уетЬ, ргопоип, агисе, ес. ТБезе раз оЁ зреесв Бе]р ю ипдег${апа ше РапсНоп оЁ 
еась \ота ш а зещепсе. Но\уеуег, 1 15 поЁ епоизВ 10 Кпо\/ Ше саесогу оГа \’ог4, езреслаПу Юг Шозе Ша тау Вауе 
тоге ап опе теапто. Зрес1ИсаПу, те \огА “Чеауез” сап Бе Ше уегб ‘Чо 1еауе” ш Ше Зт4 регзоп зпещаг ог ше 
рига] Фогта оЁ Фе поип ‘ЧеаР’, \сЬ зВошШ Бе соп$14егеа йот Фе рошЕ оЁ ]апопазе аз а зу$ет оЁ ииегге]ае4 апа 
шегерепаеп{ ипи5$. Те 14еа оЁ сопз1$епсу ш Ше 1ех1са| ап4 зетапис зрВеге оЁ 1апгиасе \’аз Нгз{ ехргеззе4 Бу 
М. М. РоКгоузКу, етрВаз171ие Фа{ “\\огз ап@ Феш теап1п25 40 поё Пуе а зерагае Ше Нот еасН оег, БиЁ аге 
соппесеа (ш опг $001), гехагез$ оЁ опг соп$с101$пезз, шо ЧШегепе отоир$, ап Фе Ъаз15 Юг Феш этопрх 15 
зипШагиу ог Апесе орроз1Яоп ш Феи Базе теаптз”. Рага@тлайс, зущаотайс, ап4 ер1@ютайс г@аНоп$ атопз 
1апопазе ип аге ипронап пап е${аНоп$ оЁ Ше зуетайс ап4 тезШаг паге оЁ 1апоцазе. Те гезеагсНегз пое фай 
У\ог4$ ещег Ше зушастайс ге]аНопз Базе оп Ше 1051са| сопасийу оЁ сопсер ап4, сопзедиепйу, Фет сотранб у 
УВ еасН оег [4]. 

\е пее4 ю ипаег$ап@ Ша от Ше рошЕё оЁ ме\ оЁ сотршег зс1епсе, зреесВ 15 поЕ зигасваге4 1шпЁогтайоп, Баё а 
зедиепсе оЁ сВагасегз. То епзиге фаЁ уо1се аа сап сопипие тю Бе изе4, Ше зреесВ гесосшиоп аррНсаНоп мапз$]айез 
{ шю {ехё. ТБе ассепть шага! шюпайоп$, ап етойопз аге атеау Бешх егазед ш Ше {ехЕ. \!МВеп даа аге 
папз]ае4 шю {ехё, Шеу аге гап$]ае4 \/ИВ 7егоз ап4 опез. 

'ТрегеРоге, сопрщегз пее4 а Баз1с ипаег\{ап те оЁ отаплтлаг {0 геЁег © # ш сазе оЁ сопРаз1оп. ТЬи$, фе пез Ююг 
Фе зкгасвге оЁ рЫгазез арреаге4. ТБеу аге а зе оЁ огаплтаг гез Бу \м51сВ а зещепсе 1$ сопзгасе4. ш Ел $8, 1 1$ 


Гогте4 \ийн фе Бер оЁ а пота! ап4 а уегб огоир. Сопз14ег Ше зещепсе, “Кае а Фе арр!е”. Неге, “Ка{е” 15 а поп 
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рВгазе ап4 “ае Фе арре” 15 а уег рвгазе. РЁегепЕ зещепсез аге югте4 изте ЧШегепЕ засигез. Аз Фе питЪег оЁ 
рВгазе згасфиге гШез шсгеазе$, а рагзе тее сап Бе сгеаже4 ю с1аз5Шу еасв \’ог4 т а рагис\аг зещепсе апа аглуе аё 11$ 


сепега! теапте (Но. 1). 


ЗЕМТЕМСЕ 


№ шип \/ога$ этоир 
Раегтттег$ 


УегЬ \ог4$ эгопр 


ргоре$5от а иотап 


Но. 1. Рагзе Тгее 


Тье зроКеп 1апопасе зубет пиезтайпе зреесв гесосп1оп ап4 зреесв зупез1$ 1$ Фе соге {есБпо]озу о{ Витап- 
сотрщег пщегасйоп, ап4 ога1 ипдег{ап@дс 15 Ше соге оЁ зроКеп 1апсиазе зузет [5]. 

АтЯйс1а| пеига1 пе\уотК$, стеае4 ш Ше Гогт о? сотршег то4е]$, соре зиссеззРиПу \ИН Ше 1а3К$ оЁ райегп 
тесорт оп. ТБеу аге ташаЫе, Шеу сап Ъе еазПу аЧаре4 (ап 11$ Ваз атеаду Бееп 4опе) ю зо!уе тапу 
ргасиса! ргоетз ге]ае4 10 зреесП гесост! оп, сопёго] оЁ уамои$ шас тез ап 4еу1сез, еуепЕ ргеЯ1сйоп, ес. 
ТБе Ъ100е5$1 Реабиге оЁ а 4еер пепга! пе \могК 1$ (гапише а 1агое атоциЕ оЁ даа, ап Феп ехгасипз спагацензИс 
шЕгтаноп. Свагасензис шогтайоп оашей Шгоиэй а пе\могК гасаге сап с1уе соо4 гези! 5 ш Ше зреесв 
сотргенеп$1оп {а$К$. 

Ап адуащасе оЁа пепга| пебуотК ш зреесЬ ргосеззте 1$ фаЁ Ше регсериоп сап регРоги @5сгитаие 1еагите 
Бегмееп зреесй ипИ$ ШФаЕ гергезепе фе оцёриЕ с1аззез оЁ Ше регсерноп. ТПе регсерёгоп поЁ ощу 1еагп$ ап орИп7ез 
Фе рагатеегз юг еасВ с1а5$ оп Фе даа Бе]опешс 10 (, Биё а[з0 И1е$ 10 ге]есЕ Фе даа Беопотшо 10 офег с1а5зез. ТБе 
регсериоп 1$ а 5гисйиге \Ий а ШеН 4естее о{ рагаПеЙзт, %сВ аПо\/з Ше изе оЁ рагаПе] Вага\уаге. ТВе Вг$Е пеига| 
пебмотК по4е]5 изе4 ш 5реесВ гесосшйоп зузепа$ \еге 4еуе]оре4 ог айс з1°па1$, ап поЁ ог Феи зедиепсе$ ог 
$12па1$ зиБ]есЕ 10 {етрога! уачабПиу. Гаег, тесштепе пеига] пеё\уотК$ ап@ сопуоийопа! пепга|! пеб\уогкз$ \ете 
ргорозеа [5]. 

Аз 1апоцасе то4е]з аге тгате4 оп 1агоег апа ]агоег {ех{5, Фе питбег оЁ итаце \огаз (Фе уосаб\Щагу) шсгеазез, 
ап4 сопНпиои$ зрасе етбедате Ве]рз ю аПеуае Ше сигзе оЁ Апптепзуопа!у ш 1апоиасе то4дейпе. Ап а[{егпае 
езсирноп 15$ Фаё а пепга| пеё арргохитайез$ Фе 1апоиазе РапсНоп. ТБе пеига| пеЁ агсоцесваге по Е Бе Ёе4 Гогмага ог 
тесигтепе, апа \ЮПе Фе огтег 15 зпирег, Ше 1авег 1$ тоге соттоп. Паргоуе4 аеогиил$, ро\егР] сотршетз, ап 
ап шсгеазе ш 41о1те4 даа Бауе Глее4 а геуошНоп ш тасЫте 1еагишо. Ап4 пе\ (есби1иез ш Фе 20105 гезаЦед ш 
“гар14 ппргоуетепЕ т {азКз” шс1ад те 1апечазе тапруабНоп, ш рагасщаг, папзКогтег — агсЬ\есвге Базе4 оп а 4еер 
1еагпте по4е!. уаз Нгзё шиоЧисед ш 2017 [6]. ТаЫе 1 ргезеп$ Ше таш (урез оЁ сштепИу ех15Нпо 1апоиасе 
по4е] о# пеига| пемогКз. 

ТаЫе 1 
Маш (урез оЁ пеига] пебмогК 1апоиасе то4е]$ 


Гапопаге то4е] Спагасет15с$ 


В1апесйопа! Епсодег Кергеземайоптз Нот Тгапзогтег$ 15 а пе\ тео оЁ ргегаиите 
ВЕКТ-Базе (2018) 1апоцазе. ВЕВТ 15 АШегепе Ъесаизе 1 15 дез1епте4 © геа@ ш Бош 4песйопз аё опсе. Озшя 
15 Мапесйопа| сара шу, ВЕКТ 1$ ргегаше оп {мо АШетепть БаЕ геажеа, МЕР ЧазКз: 
МазКеа Гапгиасе Мо4еПпх ап Мех( Зещепсе РгеФсНоп [7, 8]. 


Етбеа9115$ Кот Гапсцасе Моде] 1$ а мога етбеа4те тео4 Гог гергезепипе а зедиепсе 
ОЁ \ог4$ аз а соггезроп4те зедиепсе оф уестогз, Баё ип!Ке ВЕВТ, Ше \огА етБед4тз$ 


ЕСМ 
(201 ь ргодисед Бу е “Ъаз-оЁ-\ога5” то] 1$ а зпирИЁуше гергезетайоп. ЕГМо етБед4те$ аге 
сощех(-зепзиуе, ргодисш» ЧШегепе гергезешайопз ог \ог4$ Фаё зБаге Фе зате зреШпя 
БаЕВауе Чегет теаптпс$ (Ботопупл$) [9]. 
ОРТ 1$ а ТгапзЮгтег-БазеЯ агсЬцесйше ап@ Наше ргоседиге г паха! 1апсиазе 
СРТ (2018) ргосезз1те {а5Кз. Тгаште ЮПо\з а Гмо-%азе ргоседиге. Егзв, а 1апгиасе то4дейия об]еснуе 


13 изе4 оп Фе ищабее аа ю Пеагп Фе има! рагатеегз оЁ а пеига| пебуотк тоде|. 
биб5едиеп Ну, Фезе рагате(егз аге адар4 го а "агое! {а5К изшз Ше сотезропате 5ирегу15е4 
оесйуе [10]. 


Казакова М. А. и др. Анализ технологии обработки естественного языка: современные проблемы и подходы 


ЕЗРпе! тайу осизе$ оп еп4-ю-еп4 амютайс зреесВ гесозтНоп (АЗВ), ап4 адор8 улаеу- 
изе4 4упапис пеига! пебмотК юои$, Срашег ап РуТогсв, аз а тат 4еер 1еагиих епоше. 


Е$Рпее (2018) ЕЗРпе! ао оПо\з Ше Ка!41 АЗК оо зе Гог даа ргосезз1те, Ёеабаге ехгасйоп/огтаь, 
ап4 гес1рез 0 ргоу14е а сотр!ае зер Гог 5реесВ гесоршиоп ап оег зреесВ ргосеззше 
ехрегитепк [11]. 

Зазрег (2019) Моае! изез ошу 1. сопуош@оп$, Бас погта|7айопт, КеГ О, горо ап гез14иа1 


соппесНопз$ [12]. 


СРТ-2 (2019) 


ОРТ-2 апз]а{е$ {ехё, апз\ег5 дцезНоп$, зиплтаг17ез раззасез, ап сепегаез {ехё ошриё оп а 
1еуе] пав \уВПе зотейтез шЧ15Ипеи1$ВаЫе Нот а оЁ Витапз, сап Бесоте герениуе ог 
попзеп1са1 еп зепегаНпо 1юпх раззазез [13]. 


\МАУ2ГЕТТЕВ++ Г 15 ап ореп-зоигсе 4еер 1еагилие зреесВ тесостоп Нате\могК. \уау2]ецег+-+ 1$ ууиШеп 
(2019) епйге!у ш С++, ап4 изез фе АтгауЕпе {епзог ИБгагу Гог тахипит еЁй<епсу [14]. 
\МАУЗУЕС (2019) \ау2уес, 13 а ВОВ при пебмуотК ФаЕ {аКез ри ай 1о аз шриё ап4 сотршез а 
сепега] гергезегианоп Ша сап Бе 1приЕ © а зреесВ гесосш@оп зузет [15]. 
ТЬезе аге сгоз5-Низиа! ]1апгиасе то4е!з (ХГМ$): опе ипзирегу15е4 Фа ошу геЦез оп 
хим (2019) топоНиеиа| Ааа, ап опе зирегу1зе ФфаЕ Теуегасез рагаПе] даёа У а пе\у сгозз$-Ппеиа1 
Лапсцаее то4е! обесйуе. [ оМаштз зае-оРФе-агЕ гези5 оп сго$$-Ппеиа| <аз$Исайоп, 
ипзирегу1зе4 ап4 зирегу15е4 тасбше гап$айоп [16]. 
ХЕ Ме! изез а сепега|7е аиогестеззлуе геташите тео Фа епаез ]еагите ЫАпесНопа| 
ХЕМе сошех(з ШгоизВ пах11117105 Фе ехресе4 ПКкеПвоо4 оуег аП региаиайопз оЁ Фе ГасоптаНоп 
(2019) ог4ег ап4 аиюгесотезуе огианоп. ХЕМе{ ицеотаез 14еа5 Нот Тгапзогтег-ХГ,, 1е зае- 
о#-Фе-агЕ ашотеотезауе по4е], по гешаииае [17]. 
Ты$ парететайоп 1$ Фе зате аз ВЕКТ Моде! и а Нпу етБедЯтз»5 Е\\еаК аз \еП аз а 
ВоВЕЕТа (2019) зешр ог КоВЕВТа ыы по4ез. КоВЕВТа Баз Ше а С а ВЕЕТ Би изез 
а Буе-еуе! ВРЕ аз а юКешиег (зате аз СРТ-2) ап4 аррНез а ЧШегепе ргейаиите зсВете 
[18]. 
ЕТЕСТВА Еаеп у Геагишо Епсодег Та Ассигае!у СЛаз Нез ТоКеп Вер!асетеп 1$ а пем рге- 
(2020) 1еапше тефо@ Ша ошре{огт$ Аеу@ортепЕ езИтайоп Уфой шстеазше Фе 
сотршайопа] сое [19]. 
ТС зу%ет айп$ аё ти9-писгорВопе ишЯ-зреаКег зреесН гесост1Ноп ап4 Фамтайоп. ТВе 
зумет иН7ез зоН-асиуйу Базе оп Се Зоигсе Зераганоп (@$5$) Нош-епЯ ап а 
ТС зучет (2020) сотЫпаНоп оЁ ауапсе асоизис то4депие ‘1есбиаиез, шсааше (©55-Базе4 иашше ааа 


аистешайоп, ишШа-5614е ап4 ти-5еат зе-айепНоп ]ауегз, $аНзИсз 1ауег ап4 зресга1 
[20]. 


ОРТ-3 (2020) 


ОтНке офег по4е!5 сгеме4 тю здуе зресНс 1апиазе ргоетлз, Шеш АРГ сап зо уе “апу 
ргоеп$ ш ЕпеИзН”. ТВе а1еогибт \отК$ оп Фе рипс1р[е оЁ ащюсотр!еНоп: уои ещег Фе 
Бесшите оЁ Фе (ехё, ап4 Фе ргоэтат сепегаез Фе п1о5Е ПКеу сопипиайоп оЁ 1 [21]. 


АТВЕВТ (2020) 


АГВЕКТ шсогрогаез мо рагатаеег гедисНоп {есбаиез Фа ИЕ Фе та]ог об$асез ш 
зсаЙпе ргетатеЯ то4е!5. Тье ЙЯгзЕ опе 15 а Гасюпте4 етбед4те рагатаептайоп. Ву 
зрШипе а 1агое уосабщагу етбе41т» тайлх шю о зтаП плаблсез, 1 зерагаез Ше $17е оЁ 
Фе Ш 99еп 1ауетз гот Ше $17е оЁ уосаб\агу етбеЧЧ те. ТВе зесопа {есбтате 15 сгоз$-|ауег 
рагатеег звампе. ТЬ$ (есбилаче ргеуеп$ Ше рагатеег гот это\/те ми Фе дер оЁ Фе 
пемогк [22]. 


ВЕКТ-у\ут-ехё, 


Ргегашеа ВЕВТ хи Увое Уога МазКше дие ю Ше сотр ехйиу оЁ СЫшезе статлтаг 
зегисвте ап4 фе зетапис Фуегзиу, а ВЕКТ (\\т-ех@ \аз ргорозеа Базе оп Ше вое 


2021) Сыштезе \уогА таз ше, мН1сВ пугайез Фе ага\ФасК$ оЁ тазК1те рагна| \Мога Р!есе {юКеп$ 
ш ргенашеа ВЕКТ [23]. 
Ты$ 15 Ра \ау$ Гапоиазе Моде! 540-51Шоп рагатеег, 4епзе десо4ег. Ощу Тгапзогтег 
РаГМ (2022) то4е шашеа \и Фе Рай \ау$ зу$ет епае4 $ тю еЁйслепЙу тат а зте тоде] асго$$ 


шшар!е ТРО у4 Род$ [24]. 


А$ сап Бе зееп Нот ТаЫе 1, Фе Нг$Е гапзРогтег то4е], изте а МАнесйопа| сарабПиу, аПо\е4 мо ЧШетепе Би 
те|ае4 тазКз ог Фе МГР ю Бе зиЧеа БеГогевапа: зипу]аНпе а тазКе4 ]1апоцасе апа рге@сипе Ше пехЕ зещепсе. 
В14песйопа! Епсо4дег Вергезещаноп$ гот ТгапзРогтег$ соп$15{ оЁ (мо $ерз: Фе Нг5Е %ер 15 ргебаниие \Беге Фе 
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Чаба ещег Ше ]ауегз оЁ гапзФогтег, ап4 Ше гезш оЁ 1$ $ер аге уесюгз Юг \’огаз. ТВе зесопа %ер 15 Нпе ишише. ТВе 
ргейгатите %ер соп315{5 ОЁ мо $ерз: Фе тазКед ГМ апа Мех! Зещепсе Рге@сНоп (МЗР) [7, 8]. ВЕКТ 1$ поЕ \иош 
Нам, Ме по5Е оБу10и$ опе 15 Ше 1еагише теёфо4 — Фе пеига| пебуотК и1ез Ю сиез5 еасН \ога зерагаеу, св 
теапз ШФаё { 10зез зоте роз$1Ые соппесНоп$ Беб\мееп \’ога$ дите Ше 1еагишз ргосез5. Апо@ег опе 1$ фаЁ Ше пепга] 
пеё\уотК 15 таше4 оп тазКе4 юКепз, ап4 еп изе4 Гап4атетаПу а1етепе {а$К$, тоге сотшр]ех опез. 

Етфеаа!11$$ Нот Гапоцасзе Моде! 15 а деер сощехшаПтеЯ \уогЧ гергеземайоп а то4е!з Бош сошрМех 
спагасет1$Яс$ оЁ \отА изазе (е.о., зущах ап4 зетапйс$), ап Во\и 115 изазе уапез асгоз$ Ппе1$Яс сощех (1.е., 
то4е! роузету), зисН аз “БапК” ш “пуег БапК” апа “Ъапк Ба1апсе”. ТБезе у’ога уесюгз аге Теагпе4 РапсНоп$ оЁ ше 
егпа| $а{ез оРа Аеер Атесйопа] 1апоиасе то4е1 (Б.М), ушсв 1$ ргебате4 оп а 1агое {ехЕ согриз. ТВеу сап Бе 
еазПу а44е4 © ех1зИпх по4е|$ ап4 з1отйЯсап Иу ппргоуе Ше зе оЁ Ше ап асгозз а Бгоа4 гапое оЁ сваПепоте» МЕР 
ргоЫепл$, пса те ацезНоп апз\уегие, {еха| етаШтет, апа зептеп: апа]уз1$ [9]. 

То аЙемлае Фе ргоМет, зиегте гот Фе @15сгерапсу Бе \мееп ше ргебаште ап4 Впе-шпто асе Бесаизе Ше 
тазкшо юКеп [МАЗК] пеуег арреаг$ оп Ше Нпе-бшиз асе, ХГМе! \аз ргорозеа, \ысВ 1$ Базе оп Тгапзогтег- 
ХГ. То асШеуе И $ ©оа1, а поуе Г\мо-5еат зе{-айепйоп тесвап1$т, апд опе ю сВапее Ше ацюепсод тс ]1апоцасе 
то4е! шо ап ашогезтезяуе опе, \шсВ 1$ зиаПаг ю Ше ша41#опа! зайзИса! 1апоцазе по4е]з, \еге ргорозе4 [17]. 
КоВЕКТа, ЭТС Зумет, СОРТ то4е!$ \еге изе4 ш диИе а 1агое патЪег оЁ зуу\епаз. Ап Феу зВо\е4 ргейу гоо4 
тези 5. ТрБезе по4е!5 зисоезе4 ШаЁ ауегасше аП юКеп гергезещаНноп$ сопзу5$\{епЙу ш@исед Бебег зещепсе 
тергезеайотз$ Шап ито Фе юКеп етбе4Ате; сот ште Ше етбеЯ 4112$ оЁ Ше Бойот 1]ауег ап е юр Тауег 
ошрегогтеЯ Фе изе оЁ Ше юр Г\№о 1ауегз; ап@ погтаНие земщепсе етбед4$$ УИ а У\Юцешах аогт 
соп$1епЙу Боо$е4 Фе регЁогтапсе [18, 20, 21]. 

Те пехЕ %ер, ргобаЫу, \Ш Ъе ю заду Ше оуегзатрПпх ап4 ипдегзатр| ие оЁ {ехша| Ааа 10 паргоуе Ше оуегаП 
еп@(у гесоршиоп еЙеси. 

Вези5. ТВе апа[уз15 оЁ Фе Шегагу зопгсез 4езсгошс пе\’ ше!фо4$ оф ргосеззшс ога| зреесв, св ргоу14ез 
шЕгтаноп абои( пепга! пебуогК$, тефо@$ Рог Ше гасвиге ап зупе$1$ оЁ зреесН, та4е # роз51е ю даес{ 
фе ГоПо\ по: 

1. АЦ Фе то4е!$ ргезеще ш Фе геме\м гедише ]агое сотрийпе ромег ю зо уе пага! 1апоиасе ргосез$ те 
ргоЫетаз. [ 15 сошршайопаПу тоге ехрепзуе 4ие ю 1$ 1агоег згисвате. 

2. Мопе оЁ Ше сштепИу ех15Нпо {есбпо|о51ез епае зо!у1шх Фе ЁШ гапее оЁ ‘азК$ Гог гесо2и171пе сопНпаоц$, 
еРесйуе зреесп. 

3. Мо пабига! 1апэиасе ргосезз1ще по4е!$ аге 4ез1оте4 ю Бап е а \л4е уапе(у оЁ ЕпоП$Н Ч1аес апа 1410115. 

05си5$10п апа Сопс$10п5$. Уотсе аз315{ап(5 гергодисе ап4 гешогсе аП зегеобурез а1гогИвилз. ТВеу, аз а пише, 
тергодисе позе $егеогурез ШаЁ ех13Е по\/ ш зослейу. УВаЕ 4оез 11$ астеуетепи геаПу теап? Г теапз Фа Фе ус1се 
азз1збап{ 15 по \огзе (ог тауБе еуеп Бейег) ап ап ауегасе регзоп аё гесо211711е Ше зреесВ оЁ а регзоп \И а 
запаага Мог Атепсап ассеп. Ви ап АЁ1сап Атепсап зреаК$ го ап а3315бапф, Феп Ше ассигасу ул агор ю або 
80 %. ТЬ1$ 13 а Басе 41Шегепсе. Могеоуег, \Беп сопуегипх уо1се ю {ехЬ Ше зрес сз оЁ \упйпе, уШсВ сап Бе 
ппроцапе Гог зреаКегз, аге спагашее4 то Бе 1054. 

Уотсе а5$15(ап{5 40 по {аКе шю ассоциЕ Ше зреесВ ап4 изег Ваб1($ оЁ Пе е!4ейу ап4 реоре у зресла! пеед$. 

Апа Беге, 1 15$ поЁ еуеп а\мауз ве сотрехИу оЁ тесоет оп. ТПеге 15, е.с., засП а сопд1юп аз Аузана — а 
Теагаге оЁ Ше Рапсбоштя оЁ Фе соппеснопз Бебмееп Ше зреесВ аррагаба$ ап Фе пегуоиз зузет, \УмсН сап саизе 
ААси]нез ш ргопочист® шаг\1Аа| зоип4$ ог, ш сепега], т зреесв. 

А]50, 4ме ю Баг4маге ШшиИаНопз, апу сагелАсе \Ш тезай ш юо шапу то4е| рагатеег$ ап4 ипзиссез$ В! 
ехесиНоп. ТВе \гау тю зо]уе Ше ргоМет оЁ пшир!е сусез оё Ч1а]оие гедшшез Ва ег гезеагсВ. 
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